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The  cover  shows  a  drawing  of  the  results  of  a  surface  oil  flow  test 
conducted  in  the  USAF  Academy's  2  ft.  x  3  ft.  subsonic  wind  tunnel. 

The  model  tested  is  the  Northrop  top~mounted  inlet  vertical  assisted 
takeoff  and  landing  (VATOL)  configuration.  The  model  was  mounted  at 
a  45~degree  angle  of  attack  in  a  freestream  of  100  ft/sec.  Additional 
results  for  this  configuration  are  given  in  Major  Jonas's  article  in 
this  issue  of  the  Aeronautics  Digest. 
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PREFACE 

This  report  Is  the  eleventh  issue  of  the  Air  Force  Academy 
Aeronautics  Digest*.  Our  policy  is  to  print  articles  which  represent 
recent  scholarly  work  by  students  and  faculty  of  the  Department  of 
Aeronautics,  members  of  other  departments  of  the  Academy  and  the  Frank 
J.  Seiler  Research  Laboratory,  researchers  directly  or  indirectly 
involved  with  USAFA'^sponsored  projects,  and  authors  in  fields  of 
interest  to  the  USAFA. 

In  addition  to  complete  papers,  the  Digest  Includes,  when 
appropriate,  abstracts  of  lengthier  reports  and  articles  published  in 
other  formats.  The  editors  will  consider  for  publication  contributions 
in  the  general  field  of  Aeronautics,  including: 


•Aeronautical  Engineering 
Aerodynami cs 
Flight  Mechanics 
Propulsion 
Structures 
Instrumentation 
•Fluid  Dy nara 1 cs 

•Thermodynamics  and  Heat  Transfer 
• B iomechanlcs 
•Engineering  Education 
•Aeronautical  History 

Papers  on  other  topics  will  be  considered  on  an  Individual  basis. 
Contributions  should  be  sent  to: 


Editor,  Aeronautics  Digest 
DFAN 

US  Air  Force  Academy 
Colorado  Springs,  CO  80840 


The  Aeronautics  Digest  is  edited  at  present  by  MaJ  Jay  E.  DeJongh, 
PhD;  William  H.  Heiser,  PhD;  and  Maj  Robert  M.  Hogge,  PhD,  who  provided 
the  final  editorial  review.  Our  thanks  also  to  Contract  Technical 
Services,  Inc.  for  editing  and  illustration  services. 


"Previous  Issues  of  the  Digest  can  be  ordered  from  the  Defense  Technical 
Information  Center  (DTIC),  Cameron  Station,  Alexandria,  VA  22324. 
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EXTERNAL  FLOW-FIELD  MEASUREMENTS  ON  A  TOP-MOUNTED,  INLET-VATOL 
CONFIGURATION  AT  HIGH  ANGLES  OF  ATTACK 

Frederick  M.  Jonas* 

Abstract 

This  paper  presents  the  results  of  a  research  effort  conducted 
for  tiie  NASA  Ames  Resear'cfi  Center  in  which  a  seven-hole  pr’ assure  probe 
was  used  to  map  an  unknown  flow  field.  Also  included  in  this  effort 
are  results  that  define  the  probe's  measurement  capabilities.  Using 
a  seven-hole  pressure  probe,  external  flow-field  measurements  were 
made  on  a  1 0 1 h  -  sea  1  e  ,  top-mounted  inlet  model  of  the  Northrop  VATOL 
(vertical  assisted  taxeoff  and  landing)  configuration.  The  tests  were 
conducted  in  the  USAF  Academy's  2  ft.  x  3  ft.  Subsonic  Wind  Tunnel. 

The  model  was  mounted  at  a  fixed  it5-degree  angle  of  attack  and  the 
test-section  velocity  was  nominally  100  ft/sec  for  all  subject  tests. 
External  flow-field  measurements  were  made  over  most  of  the  upper 
surface  of  the  wing  and  fuselage.  The  results  show  the  presence  of 
two  counter-rotating  vortices  produced  by  the  wings'  leading-edge 
extensions.  These  vortices  help  keep  the  flow  over  the  wing  and 
fuselage  upper  surfaces  attached.  Vortex  disintegration  appeared  to 
occur  Just  aft  of  the  leading  edge  of  the  inlet  compression  ramp. 

I .  Introduct ion 

This  research  effort  was  conducted  to  accomplish  two  objectives: 
first,  to  measure  the  external  flow  field  in  and  around  the 
top-mounted  inlets  on  a  1 /i)0th  scale  model  of  the  Northrop  VATOL  at  a 
fixed  115-degree  angle  of  attack  using  the  seven-hole  pressure  probe 
and  second,  to  determine  the  measurement  capabilities  of  the 
seven-hole  pressure  probe  in  high  and  low  shear  flow  environments. 

All  tests  were  conducted  in  the  USAF  Academy's  Aeronautics 
Laboratory's  2  ft.  x  3  ft.  Subsonic  Wind  Tunnel  at  a  nominal  test 
section  velocity  of  100  fps. 

The  Northrop  VATOL  configuration  is  a  top-mounted  inlet  fighter 
aircraft  with  a  delta-wing  planform  and  wing  leading-edge  extensions 
(LEX).  Since  this  configuration  was  designed  to  operate  at  high 
angles  of  attack,  an  examination  of  the  interactions  of  the  LEX 
vortices  with  the  top-mounted  Inlets  during  these  conditions  of  flight 
is  valuable.  The  seven-hole  pressure  probe  was  used  to  map  the 
external  flow  field  on  the  upper  surface  of  the  wing  and  fuselage, 

•Major,  USAF,  Associate  Professor  of  Aeronautics,  DFAN 
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Starting  just  forward  of  the  wing  LEX  and  extending  over  the  entire 
body.  Both  qualitative  and  quantitative  Information  on  these 
interactions  was  obtained.  Crossflow  velocities  and  local  total 
pressures  are  presented  in  graphical  form  for  the  data  planes 
investigated  (Figures  12  and  14).  Surface  oil  flow  results  for  this 
test  condition  are  also  presented  (Figure  15). 

In  addition,  as  part  of  this  research  effort,  the  measurement 
capabilities  of  the  seven“hole  pressure  probe  in  high  and  low  shear 
flows  were  determined.  It  should  be  noted  that  this  test  was  part  of 
an  ongoing  research  effort  at  the  USAF  Academy  to  establish  the 
capabilities  and  limitations  of  the  seven-hole  pressure  probe  in 
mapping  unknown  flow  fields.  The  results  of  our  research  show  this 
probe  to  be  a  valuable  and  reliable  instrument  in  making  these  kinds 
of  investigations.  For  this  test,  seven-hole  pressure  probe 
measurements  were  compared  to  hot-wire  anemometer  results  in  the  wake 
flow  behind  an  airfoil  (low  shear  flow).  Next,  seven-hole  pressure 
probe  measurements  of  local  total  pressure  at  a  point  in  the  vortex 
wake  (high  shear  flow)  created  by  the  wing  LEX  of  the  Northrop  VATOL 
model  were  compared  to  measurements  of  local  total  pressure  obtained 
using  a  total  pressure  probe.  For  this  last  effort  the  total  pressure 
probe  was  aligned  first  with  the  freestream  flow  and  then  with  the 
local  flow. 

The  results  of  these  Investigations  appear  in  Section  V.  A 
review  of  the  basic  flow  characteristics,  in  terms  of  velocity  and 
pressure  in  viscous  wake  phenomena,  appears  in  the  Theory  section 
below.  This  section  serves  as  a  review  of  what  one  should  expect  to 
see  when  measuring  these  wake  phenomena.  In  the  following  section, 
the  apparatus  and  procedure  used  in  the  conduct  of  this  investigation 
are  described.  The  paper  closes  with  a  description  and  summary  of  the 
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results  of  the  investigation. 

1 1 .  Theory 

In  the  measurement  of  wake  phenomena,  one  has  certain 
expectations  regarding  pressures  and  velocities.  For  example,  in  the 
measurement  of  an  airfoil  wake  in  a  steady  incompressible  viscous 
flow,  one  expects  the  velocity  profile  through  the  wake  to  appear  as 
shown  in  F i gur e  1  . 


y 


Figure  1.  Velocity  Profile  Through  an  Airfoil  Wake 

Depending  on  how  far  behind  the  airfoil's  trailing  edge  measurements 
are  made,  one  may  or  may  not  see  any  variation  in  the  local  static 
pressures  through  the  wake.  In  fact,  beyond  a  distance  of  less  than 
one  chord  length  behind  the  trailing  edge,  changes  in  the  local  static 
pressure  through  the  wake  are  negligible.  On  the  other  hand,  the 
local  total  pressure  deficit  through  the  wake  is  expected  to  have  the 
same  shape  as  the  velocity  profile  regardless  of  distance  behind  the 
trailing  edge.  The  local  total  pressure  deficit,  like  the  velocity 
profile,  would  then  diffuse  or  decay  slowly  as  one  moves  farther  away 
from  the  airfoil's  trailing  edge,  since  the  local  total  pressure  at 
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any  point  is  a  measure  of  both  local  velocity  and  static  pre;')3ure.  At 
any  point  in  a  steady  incompressible  flow 


+ 


(  I 


The  variation  of  P,,  through  the  wake  is,  of  course,  caused  by  viscous 
effects . 

We  would  expect  to  see  similar  trends  in  a  simple,  two-dimensional 
viscous  vortex.  These  trends  are  easily  verified  by  using  a  model:  a 
vortex  filament  that  dissolves  with  time  because  of  fluid  friction 
(Ref.  1 ) . 

Using  a  two-dimensional,  t i me -dependen t  viscous  vortex  to  model 
the  characteristics  of  a  three-dimensional,  v i scous - vor t ex  filament 
in  steady  flow  may  at  first  seem  confusing,  since  we  will  be  measuring 
the  latter.  The  decay  of  a  three-dimensional  viscous  vortex  with 
downstream  distance  in  a  steady  flow,  however,  can  be  related  to  the 
decay  of  a  two-dimensional  vortex  with  time  according  to  the  relation 
that  distance  equals  velocity  times  time,  or 


t  =  dA„  or  d  =  t 

40  00 


Thus,  in  the  discussion  that  follows,  for  the  two-dimensional  vortex, 
time  can  be  replaced  with  distance  along  the  three-dimensional  vortex 
filament.  The  change  in  the  properties  with  time  described  for  the 
two-dimensional  vortex  thus  corresponds  to  the  change  in  properties  for 
the  vortex  filament  as  the  vortex  proceeds  downstream  along  its 
length.  The  variation  in  the  circumferential  velocity  component  with 
respect  to  time  is  described  by 


U,  (r.t) 


2irr 


( 


where  T,  is  the  circulation  around  an  infinitely  large  circle  (Ref. 
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1).  The  change  in  velocity  with  respect  to  time  is  shown  in  Figure  2 


Figure  2.  Variation  of  Velocity  with  Radius,  Viscous  Vortex  Filament 


Note  that  as  the  center  of  the  vortex  is  approached,  the  velocity 
approaches  0,  and  as  the  radius  becomes  very  large,  the  velocity  again 
approaches  0.  From  momentum  considerations  one  can  show  that 


J, 

T 


(A 


which  describes 
f i lament .  Note 


If  =  P  “  =  7  I  ^ 


the  variation  of  local  static  pressure  through  the 
that  for  Eqn.  (5) 


(5 


iP 

3r 


>  0  V  r  >  0 


(h 


and 


-2^  -►0  as  r  -►  ® 
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or,  that  the  static  pressure  increases  from  some  in  i  n  i  in  u  v  i  ’ 


moves  away  from  the  vortex  center  and  that  the  static  pr,>:; 


off  far  away  from  the  vortex  center.  Finally,  if  one  as  suin-’ t  F  .  1  it 
any  point  in  this  incompressible  viscous  flow  the  local  total  ;>  r  i- s  s  ij  r 
is  the  sum  of  the  local  static  and  dynamic  pressures,  or 

P_  =  P|  + 

Ul 

then  the  change  in  the  local  total  pressure  with  radius  is  d(> scribed 


°'^o.  3P,  3U 

=  — 

Substituting  Eqn,  (5)  into  this  expression  for  '•P/'ir  and 
differentiating  Eqn.  (3)  to  obtain  1 U,,/'.  r  results  in  the  following 
equation,  describing  the  change  in  local  total  pressure  through  the 


vortex  . 


LfiV  J.e 


Again,  note  that 


>  0  V  r  >  0 


'  0  as  r  -*■  »  and  as  r  -*■  0 


or  that,  as  befor'e  with  the  local  static  pressure,  the  local  total 
pressure  increases  from  some  minimum  value  as  one  moves  away  from  the 
center  of  the  vortex  and  then  levels  off.  From  these  descriptions  we 
would  expect  local  total  and  static  pressures  to  vary  through  the 


vortex  filament,  as  shown  in  Figure  3.  at  any  time,  t. 


Figure  3*  Variation  of  Pressure  with  Radius,  Viscous  Vortex  Filament 


A  three-dimensional  vortex  in  steady  incompressible  viscous  flow, 
although  more  complicated,  would  be  expected  to  exhibit  the  same  local 
variations  in  velocity  and  static  and  total  pressures  as  the 
two-dimensional  vortex  filament  would,  as  previously  discussed.  The 
pressure  boundary  conditions  are  represented  by  the  freestream 
inviscid  values  of  P_  and  .  At  no  time  would  one  expect  to  see  the 

oc- 

local  total  pressure  in  the  vortex  exceed  the  freestream  values  (Ref. 
2).  Intuition  tells  one  that  in  any  viscous  phenomenon,  the  effects 
of  viscous  interactions  or  fluid  friction  can  be  expected  to  dissipate 
energy,  causing  a  loss  in  local  total  pressure  through  the  wake. 

These  changes  in  pressures  through  viscous  wakes  are  why  pressure 
probes  can  be  used  to  measure  wake  phenomena  (Ref.  2)  and  map  their 
locations. 

III.  Apparatus 

A .  Wind  Tunnel 

The  Subsonic  Wind  Tunnel  in  the  Aeronautics  Laboratory  of  the 
USAF  Academy  is  a  con t i nuous -f low ,  c losed -c 1 r cu i t  wind  tunnel 
operating  at  atmospheric  pressures.  The  test  section  is  fixed  in 


geometry,  nominally  2  ft.  high  by  3  ft.  wide  and  70  in.  in  length; 


i;SAFA-TR-a4-3 


its  sidewalls  diverge  approximately  1/2  degree  with  respect  to  the 
tunnel  centerline.  The  Mach  number  range  is  0.04  M  •  0 .  3^  >  which 

corresponds  to  a  speea  of  50  to  400  fps.  The  corresponding  Reynolds 
number  per  foot  is  0.2  (10^)  to  1.3  (10*’)  and  the  dynamic  pressure 
ranges  from  1.8  psf  to  130  psf.  A  schematic  of  the  wind  tunnel  is 
shown  in  Figure  4. 

B.  Model 

A  1/40th  scale  model  of  the  Northrop  VATOL  cor. figuration 
(Refs.  3  and  4)  was  used  in  this  investigation.  Flow  into  the 
flow-through  inlets  of  the  model  was  adjusted  to  match  conditions 
typical  of  0.9  Mach  cruise,  Ao/Ac  -  0.66,  for  all  tests.  The  model 
was  strut-mounted  at  a  fixed  angle  of  attack  of  45  degrees  with 
respect  to  the  aircraft  horizontal  reference  line  (HRL).  All  tests 
were  conducted  at  a  nominal  test-section  velocity  of  100  fps.  A 
schematic  of  the  model  is  shown  in  Figure  5,  while  a  typical 
installation  is  shown  in  Figure  6. 

C.  Pressure  Probes 

A  0.109  in  00  seven-hole  pressure  probe  (Ref.  5)  was  used  to 
map  the  external  flow  about  the  model.  The  probe  was  aligned  with  the 
freestreara  flow  direction  for  all  of  the  tests,  as  shown  In  Figure  6. 
In  subsonic  compressible  flows,  relative  flow  angles  of  up  to  70 
degrees  with  respect  to  the  probe  centerline  are  measurable  to  within 
±1.6  degrees,  while  the  ratio  of  local  total  pressure  to  local  dynamic 
pressure  is  measurable  to  within  ±5  percent  at  95  percent  certainty 
(Ref.  5).  The  seven-hole  pressure  probe  data  was  collected,  reduced, 
and  displayed  by  means  of  the  seven-hole  probe  data  acqu'sition  system 
(Ref.  5),  which  uses  the  POP  11/45  computer  system.  The  probe 
traverse,  axis  system,  and  nomenclature  are  shown  in  Figure  7. 
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Two  separate  total  pressure  probes  were  used  to  measure  local 
total  pressures  at  specific  points  above  the  model  for  later 
comparison  with  the  seven-hole  pressure  probe  measurements.  The 
first  total  pressure  probe,  shown  in  Figure  8,  is  a  0.125  in  OD 
hemi spher i cal -head  pitot  tube.  This  probe  was  aligned  with  the 
freestream  flow  for  all  tests.  The  second  total  pressure  probe, 
consisting  of  a  0.109  in  OD  and  0.086  in  ID  tube,  was  bent  at  a  23.8+ 
0.1  degree  angle  to  correspond  to  the  local  flow  angle  (measured  using 
the  seven-hole  pressure  probe)  at  a  specific  point  above  the  model. 

The  tunnel  installation  at  this  point  is  shown  in  Figure  9.  Pressure 
data  for  both  probes  was  collected  using  the  PDP  11/45  data 
acquisition  system. 


D.  Airfoil  Model 

An  NACA  0015  airfoil  model  was  used  to  obtain  wake  measurement 
comparisons  between  the  seven-hole  pressure  probe  and  the  hot-wire 
anemometer  in  the  subsonic  wind  tunnel.  The  airfoil  model  spanned  the 
wind  tunnel,  thus  representing  two-dimensional  flow.  The  model's  chord 
length  was  6  in.  and  spanned  the  wind  tunnel  from  the  top  to  the 
bottom  walls.  The  airfoil  was  set  near  0  degrees  angle  of  attack,  and 
tests  were  conducted  at  60  fps  and  100  fps. 


E.  Hot-Wire  Anemometer 

A  constant-temperature,  hot-wire  anemometry  system 
(Thermo-Systems  Model  1050)  belonging  to  the  Frank  J.  Seiler  Research 
Laboratory  (an  AF  Systems  Command  tenant  organization)  was  used  to 
make  airfoil  wake  velocity  measurements.  A  single-wire  sensor  probe 
was  employed  to  measure  freestream  velocity  components  through  the 
airfoil's  wake.  Linearized  anemometer  signals  were  processed  using 


signal  conditioning  circuits  supplied  with  the  system.  Data  wa 
recorded  using  a  Houston  Instrument  Model  2000  x-y  recorder. 


Figure  8.  Total  Pressure  Probe  Aligned  with  Freestream 
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IV .  Procedure 

A.  Flow  Surveys 

1.  Seven-Hole  Probe/Northrop  VATOL  Model 

The  seven-hole  pressure  probe  was  used  to  collect  flow 
field  data  at  18  data  planes  chosen  to  survey  flow  regions  in  and 
around  the  top-mounted  inlets.  The  probe  was  parallel  to  the 
freestream  velocity  ( a  =  0  degrees)  for  all  data  planes.  Each  data 
plane  was  defined  by  a  constant  value  of  Z  (Figures  7  and  10),  and 
data  was  collected  in  each  plane  normal  to  the  freestream  flow 
(defined  by  x-y  planes).  The  reference  point  for  all  data-grid  points 
(x-0,  y-0,  z=0)  was  the  centerline  top  of  the  inlet  lip,  as  shown  in 
Figure  6  and  10.  All  data  planes  consisted  of  200  grid  or  data  points 
for  a  total  of  3600  data  points.  The  nomenclature  pertaining  to  each 
data  plane  is  defined  as  follows: 

Z  M/P  xPxx 


location,  in.,  constant  value  of  Z  defining 
data  plane  location 

minus  or  plus,  corresponding  to  axis  system  with 
respect  to  reference  point 

Z-ax  1  s 


For  example,  ZM5P25  corresponds  to  a  data  plane  located  at  Z  =  -5.25 
in.  Data  planes  are  shown  in  Figure  10  and  are  listed  in  Table  I. 


2.  Pressure  Probes/Northrop  VATOL  Model 

For  the  measurement  of  local  total  pressures  both  the 
seven-hole  pressure  probe  and  total  pressure  probes  were  used  for  the 
purpose  of  comparison.  The  measurements  were  made  in  the  data  plane 
ZMP75.  Regions  of  local  total  pressure  greater  than  the  freestream 
total  pressure  were  first  identified  using  the  seven-hole  pressure 
probe.  These  regions  are  shown  in  Figure  11  and  coincide  with  the 


16 
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outboard  portions  and  top  centerline  of  the  vortices.  Data  grid 
points  and  local  flow  angles  in  these  regions  were  then  identified  for 
further  investigation  using  the  total  pressure  probes.  For  the  first 
test,  a  total  pressure  probe  was  aligned  with  the  freestream  flow  (as 
was  the  seven-hole  pressure  probe)  and  10  measurements  were  made  at 
three  identified  grid  points.  For  the  second  test,  a  total  pressure 
probe  was  then  aligned  with  the  local  flow  at  one  selected  grid  point 
(Figure  11)  to  measure  the  actual  local  total  pressure  associated  with 
the  total  velocity  vector.  A  total  of  69  data  points  was  recorded  at 
this  grid  point.  Note  that  each  data  point  represents  an  average  of 
50  samples  rather  than  an  Instantaneous  value. 

3.  Airfoil  Model  Test 

For  the  wake  Investigations  of  the  NACA  0015  airfoil  model, 
both  the  seven-hole  probe  and  hot-wire  anemometer  were  used.  The  data 
was  collected  3  chord  lengths  downstream  from  the  model.  Each  survey 
consisted  of  a  data  pass  through  the  airfoil  wake  in  a  direction 
normal  to  the  airfoil  wake  centerline.  Both  probes  were  parallel  to 
the  freestream  velocity,  and  data  was  collected  at  the  model  midspan 
to  minimize  wind-tunnel  wall  effects. 

B.  Surface  Oil  Flow 

Surface  oil  flows  were  conducted  on  the  Northrop  VATOL  model 
mounted  at  *15  degrees  angle  of  attack.  Test  section  velocity  was  100 


Figure  10. 


Data  Plane  Locations,  Grid  Sizes 


Data  Plane 

Grid  Size 
(x  by  y,  inches) 

L oc a  t  ion 

ZP6 

16x8 

1  inch  aft  of  nozzles 

ZP2 

16x8 

Junction  vertical  tail 

LE  and  f uselage 

ZP1 

1  Mx7 

— 

ZO 

8x14 

Reference  plane,  top 
of  inlet  lip 

ZMP25 

8x4) 

— 

ZMP5 

8  X 14 

— 

ZMP75 

8x44 

LE  of  lower  inlet  ramp 

ZMl 

8x44 

— 

ZM1  P5 

8x44 

— 

ZM2 

8x14 

Intersection  of  LEX 
and  wing  LE 

ZM2P5 

7x3.5 

— 

ZM3 

6x3 

— 

ZM3P25 

5x2.5 

— 

ZM3P5 

5x2,5 

— 

ZM3P75 

44x2 

— 

ZM«P25 

44x2 

— 

ZMHP75 

14x2 

Approximately  where  LE 
of  LEX  joins  fuselage 

ZM5P25 

2x1 

Approximately  where 
canopy  LE  Joins 
fuselage 
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Figure  11.  Regions  of  Positive  ,  Data  Plane  ZMP75 


V .  Discussion 

In  this  section,  the  results  of  this  research  effort  will  be 
presented  and  discussed  as  follows: 

A.  Seven-hole  pressure  probe  external  flow  field  measurements, 

1/40th-3cale  Northrop  VATOL  model  at  45  degrees  angle  of  attack, 
test  section  velocity  of  100  ft/sec.  Surface  oil  flow  results 
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will  also  be  presented. 


B.  Seven-hole  pressure  probe  and  hot-wire  anemometer  measurements  of 
an  airfoil  wake,  NACA  0015  airfoil,  at  a  test  section  velocity  of 
60  and  100  ft/sec,  nominal  angle  of  attack  near  0  degrees. 

C.  Investigations  of  local  total  pressure  in  a  vortex  wake  using  the 
seven-hole  pressure  probe,  pitot  probe  aligned  with  freestream 
flow,  and  pitot  probe  aligned  with  local  flow.  Testing  was 
conducted  on  Northrop  VATOL  model  mounted  at  45  degrees  angle  of 
attack,  test  section  velocity  of  100  ft/sec. 

In  reporting  the  results  of  the  pressure  probe  testing  the 


nondimens  ional  coefficient,  ,  will  be  presented  where 


^Totai  =  ;  ^ 

Oco  " 

=  local  total  pressure 

P^  =  freestream  total  pressure 
00 

P^  =  freestream  static  pressure 


( 13) 


A.  Seven-Hole  Pressure  Probe/Northrop  VATOL  Model 

The  results  of  the  external  flow  field  measurements  on  the 
Northrop  VATOL  model  mounted  at  45  degrees  angle  of  attack  at  a  test 
section  velocity  of  100  ft/sec  are  presented  in  Figures  iPa-i,  13f  and 
I4a-c.  Each  data  plane  consists  of  200  individual  test  or  grid 


points. 

Figures  12a-i  present  the  crossflow  velocities  and  then  contours 
first  for  negative  and  then  for  positive  values  of  .  The  data 

is  presented  by  data  plane,  starting  from  the  nose  of  the  model  and 
moving  downstream.  The  data  planes  are  viewed  looking  upstream  from 
behind  the  model.  As  seen  from  the  plots  of  crossflow  velocity,  there 
is  a  slight  side-wash  or  positive  sideslip  in  the  flow.  This  was  due 
to  the  model  mounting  in  the  wind  tunnel  and  not  to  a  flow  condition 
in  the  test  section,  -iS  was  later  verified.  In  general,  the  magnitude 
of  this  sideslip  was  negligible  and  did  not  appreciably  affect  the 
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qualitative  results.  Clearly  identifiable  in  the  cr oss f 1 ow- ve 1 oc  i  t y 
plots  is  the  circulatory  nature  of  the  flow  corresponding  to  the 
formation  of  the  f us  el  age / 1 ead i ng  edge  extension  (LEX)  vortices.  The 
vortices  are  counter-rotating  with  the  direction  of  rotation  being 
counterclockwise  for  the  right  vortex  and  clockwise  for  the  left 
vortex . 

As  previously  discussed  in  the  Theory  section,  the  vortex  wakes 
are  clearly  evident  in  the  contour  plots  of  Figure  12  for  the  negative 
values  of  C  .  The  contours  (constant  values  of  C  )  are 

generally  elliptical  in  shape.  The  values  of  proceed  from  0 

at  the  boundaries  to  a  minimum  negative  value  of  that  should 

correspond  to  the  vortex  core.  Note  that  in  these  figures  both  the 
minimum  and  maximum  values  of  for  each  data  plane  are 

presented,  as  well  as  their  location  (x,y).  The  vortex  cores  do  not 
clearly  appear  until  data  plane  ZM3P25  (Figure  12c),  just  aft  of  the 
canopy.  Their  appearance  is  indicated  by  the  shape  of  the  contours 
and  the  maximum  negative  value  of  .  The  vortices  then  steadily 

grow  in  size  (or  diffuse)  as  one  moves  downstream  towards  the  inlets. 
Multiple  low-energy  cells  first  appear  at  data  plane  ZM1 ,  as  shown  by 
the  negative  contour  values  of  .  This  data  plane  is  located 

Just  in  front  of  the  inlet  ramp.  The  presence  of  these  low-energy 
cells,  as  well  as  their  growth  and  diffusion,  is  more  apparent  in  data 
planes  ZMP5,  ZMP75,  ZO,  ZP1,  ZP2,  and  ZP6.  At  this  high  angle  of 
attack,  the  vortices  have  probably  burst  at  or  near  the  LEX  wing 
Juncture;  thus,  what  we  see  in  these  data  planes  is  the  interference 
and  blockage  effects  of  the  inlet.  These  effects  are  also  evident  in 
the  c r OSS f 1 ow- ve 1 o c i t y  plots  as  the  vortices  lose  their  definition  due 
to  this  interference  and  the  mixing  of  the  burst  vertex  cores  with  the 
local  freestream  air. 
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Figure  13.  Vortex  Core  Locations  from  Minimum  Values  of  C 
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The  trace  of  the  vortex  cores  corresponding  to  the  location  of 
the  maximum  negative  values  of  is  shown  in  Figure  13  for  data 

planes  ZM3P25  to  ZP2.  As  seen  in  this  figure,  the  vortex  cores  seem 
to  originate  from  the  point  at  which  the  leading  edge  of  the  LEX  joins 
the  fuselage.  The  decay  in  the  strength  of  the  vortices  as  one  moves 
downstream  is  evidenced  by  the  decreasing  absolute  magnitude  of  the 
value  of  Cjotal  identified  with  the  vortex  core.  These  values  are 
listed  in  Table  II  by  data  plane  with  the  location  of  the  particular 
vortex  associated  with  the  minimum  value  of  (L  being  the  left 

vortex  and  R  being  the  right  vortex)  given.  The  approximate  location 
of  the  vortex  cores,  in  inches,  is  also  presented  with  respect  to  the 
probe  axis  system  for  each  data  plane.  The  values  in  parentheses 
indicate  approximate  or  average  locations. 

Shown  in  Figure  1iia-c  are  axonometric  or  three-dimensional 
projections  of  the  data  for  each  data  plane  (note  that  the 

order  of  presentation  is  reversed).  All  projections  were  plotted 
using  the  same  scale  factor,  and  the  growth  (diffusion)  and  decay 
(dissipation)  of  the  vortices  is  clearly  visible.  Further  evidence  of 
mixing  is  the  relative  flattening  of  the  peaks  of  the  values 

(negative  being  up  or  out  of  the  plane  of  projection). 

Table  II 

VORTEX  CORE  LOCATIONS 


Data  Plane 

Maximum  Negative 
^ Total 

Locat ion  ( x 
Left  Vortex 

,  y  -  inches ) 

R i ght  Vortex 

ZM3P25 

-2 

.80R 

1  0 

50, 

1  2 

00 

1  2 

20, 

1  2 

00 

ZM3 

-2 

.  41  L 

1  0 

60, 

1  1 

70 

1  2 

40, 

1  1 

70 

ZM2P5 

-2 

.57L 

1  0 

20. 

1  1 

40 

1  2 

50, 

1  1 

30 

ZM2 

-2 

.  49L 

10 

'15. 

1  1 

40 

1  2 

60  , 

1  1 

40 

ZM1  P5 

-2 

.  28L 

1  0 

00, 

1  0 

90 

1  2 

30, 

1  1 

30 

ZM1 

-2 

.  20R 

( 1  0 

00  , 

1  0 

70) 

1  2 

70, 

1  0 

70 

ZMP75 

-  2 

.  1  2R 

(10 

00  , 

1  0 

70) 

1  2 

70, 

1  1 

00 

ZMP5 

-  2 

.07R 

( 1  0 

00, 

1  0 

70  ) 

1  2 

70, 

1  0 

70 

ZMP25 

-  2 

,06R 

( 1  0 

00  , 

1  0 

.70) 

(  1  2 

70, 

1  0 

70) 

Z 

-2 

.OOL 

(10 

00, 

1  0 

70) 

(  1  2 

60  , 

1  0 

00  ) 

ZPI 

-  1 

,78L 

(  9 

80  , 

8 

90) 

(  1  2 

80  , 

9 

60) 

ZP2 

-  1 

.841 

(  8 

70. 

9 

80) 

(  1  2 

80  , 

8 

90  ) 

ZP6 

-  1 

.  6  3R 

- 

-- 

- 

-- 
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Surface  oil  flow  results  at  the  same  test  conditions  are  shown 
Figure  15.  The  traces  of  the  counter-rotating  vortices  producing  an 
outwash  and  the  separation  line  near  the  leading  edge  of  the  LEX  are 
clearly  visible.  The  surface  oil  was  applied  rather  liberally  to 
determine  if  there  were  any  strong  reverse  flow  patterns  near  the 
trailing  edge  of  the  wing.  As  can  be  seen,  the  flow  has  fully 
separated  over  the  entire  aft  portion  of  the  wing  with  no  indication 
of  reverse  flow.  In  fact,  the  main  lifting  surface  shows  regions  of 
separated  flow  except  near  the  inlets.  The  flow  around  the  inlets 
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remains  attached  from  the  upper  surface  of  the  inlet  nacelle  to 
approximately  the  location  of  the  horizontal  tail.  The  flow  on  the 
nozzle  afterbody  has  separated.  Interestingly,  there  is  a  separation 
line  on  the  fuselage  nose  that  appears  to  correspond  to  the  separation 
vortex  forming  off  the  nose. 

Positive  values  of  Cjoj-gi  are  also  shown  in  Figures  12a-i.  The 
location  of  these  regions  of  positive  C Total  is  consistent  from  data 
plane  to  data  plane.  These  regions  of  positive  correspond  to 

the  regions  of  upwash  associated  with  each  vortex  and  the  downwash 
region  where  the  two  vortices  meet  above  the  fuselage.  The  strength  and 
location  of  these  regions  (same  data  planes  as  identified  with  the 
vortex  cores)  are  shown  in  Table  III. 

Table  III 

POSITIVE  locations 


Data  Plane 


ZM3P25 

ZM3 

ZM2P5 

ZM2 

ZM1  P5 

ZM1 

ZMP75 

ZMP5 

ZMP25 

Z 

ZP1 

ZP2 

ZP6 


Maximum 
^  Total 


Location  (x,y  -  inches) 


11.35,  13.39  centerline,  downwash 

11.35,  13.37  "  " 

11.35,  12.94 

13.52,  10.50  right  vortex,  upwash 
14.01,10.00  "  " 

14.91,  9.40  " 

11.29,  12.51  centerline,  downwash 
1 1  .29,  1 2.4  6  "  " 

11.29,  12.46  "  " 

11.35,  12.41  "  " 

11.35,  11.99  " 

11.35,  11.54 

12.24,11.33  "  " 


In  general,  it  seems  that  the  maximum  positive  C^o^al  values  first 
grow  and  then  remain  constant  until  after  passing  over  the  inlet  lip. 
A  value  of  0.45  'for  corresponds  to  a  local  total  pressure  0.2 

percent  higher  than  the  freestream  total  pressure  at  the  nominal  test 
conditions.  This  value  of  has  a  magnitude  greater  than  the 


expected  data  uncertainties  at  the  95  percent  confidence  level  (Ref. 
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5).  That  these  readings  exist  and  are  consistent  from  data  plane  to 
data  plane  indicates  that  they  either  exist  or  are  the  result  of  an 
idiosyncrasy  of  the  seven-hole  pressure  probe.  The  second  alternative 
appears  more  likely  at  this  time;  that  is,  the  positive 
readings  result  from  the  physical  interaction  and  limitations  of  the 
probe  in  regions  of  high  shear  flow.  Several  studies  of  the 
seven-hole  pressure  probe  are  being  conducted  at  present  to  quantify 
these  limitations.  The  results  of  two  of  these  efforts  --  a  correlation 
of  the  seven-hole  probe  measurements  with  hot-wire  anemometer 
measurements  in  an  airfoil  wake  or  low  shear  flow  and  an  attempt  to 
measure  the  actual  pressures  in  these  positive  regions  of  using 

total  pressure  probes  --  appear  in  Sections  B  and  C  below. 

B.  Airfoil-Wake  Measurements 

Measurements  of  the  velocity  profile  through  the  wake  of  an 
NACA  0015  airfoil  were  made  at  a  location  3  chord  lengths  aft  of  the 
airfoil's  trailing  edge.  The  airfoil  model  was  set  at  a  nominal  angle 
of  attack  near  0  degrees,  and  the  tests  were  conducted  at  test-section 
velocities  of  60  ft/sec  and  100  ft/sec.  Measurements  of  the  velocity 
profile  (velocity  deficit)  were  made  using  a  single-wire,  hot-wire 
anemometer  and  the  seven-hole  pressure  probe.  The  wake  surveys  were 
made  at  the  centerline  spanwlse  station  of  the  airfoil  model.  The 
seven-hole  pressure  probe  data  consisted  of  three  passes  through  the 
airfoil  wake  at  this  location.  Each  pass  through  the  wake  was  in  the 
opposite  direction  from  the  previous  pass.  Although  the  pressure  data 
is  not  presented,  there  was  no  evidence  of  positive  values  of 
through  the  wake,  and  the  static  pressures  were  very  nearly  equal  to 
the  freestream  static  pressure  at  this  data  plane  location.  This 
correlation  is  to  be  expected  from  the  theoretical  model  presented 


b’SAFA-TR-84-5 


earlier.  As  can  be  seen  from  Figure  16,  the  velocity  profiles  compare 
very  well,  especially  at  the  higher  velocity  (100  ft/sec),  as  one 
would  expect.  The  seven-hole  pressure  probe  data  is  within  3*^ 
percent  at  60  ft/sec  and  0.7  percent  at  100  ft/sec  of  the  peak 
velocity  deficit  measured  by  the  hot-wire  anemometer.  Both  values  are 
within  the  95  percent  confidence  level  of  uncertainty,  based  on  Mach 
number,  of  the  seven-hole  pressure  probe  data  acquisition  system  (Ref. 
5).  There  was  no  evidence  of  hysteresis  in  the  seven-hole  pressure 
probe  measurements  based  on  the  direction  of  the  probe' --  traverse 
through  the  wake. 

C.  Positive  Measurements 

Measurements  of  local  total  pressure,  and  thus  C  ,  were 

made  using  three  different  pressure  probes  at  the  data  plane  ZMP75  on 
the  Northrop  VATOL  model.  Regions  of  positive  local  total  pressure 
coefficient,  Cjotal  ’  were  identified  using  the  seven^hole  pressure 
probe,  as  shown  in  Figure  11.  At  specific  points  in  these  regions, 
two  total  pressure  probes  (one  aligned  with  the  freestream  flow  and 
one  aligned  with  the  local  flow)  were  inserted  to  measure  the  local 
total  pressure.  The  results  for  the  point  (11.80,  12,50)  in  the  data 
plane  ZMP75  are  summarized  in  Table  IV.  It  should  be  noted  that  the 
local  flow  angle  at  this  point  was  determined  from  the  seven-hole 
pressure  probe  data  to  be  m  =  -23. *11  degrees  with  negligible  sideslip. 
The  angle  of  the  bent  probe  was  23.8  degrees,  which  is  within  the  data 
uncertainty  of  the  seven-hole  pressure  probe  measurements. 


Seven -Hole 
Pressure  Probe 

Total  Press u 
A1  igned  with 

F r ees  tr earn 

re  Probe 

Aligned  with 
Local  Flow 

Number  of 
Readings  ( 50 
Data  Samples 

Per  Reading 

1 

1  0 

69 

Average  Value 
^  Total 

+0.197 

-0.096 

-0.013 

Standard 
Deviation,  a 

+0.003 

+0.0237 

+0.023 

Maximum  Value 

C Total 

— 

-0.003 

+0.090 

Minimum  Value 
^ Total 

— 

-0.067 

-0.055 

Percentage  of 
Positive  CTotal 
Readings 

... 

0% 

26% 

As  can  be  seen  from  the  data,  the  total  pressure  probe 
measurements  of  the  positive  value  of  C'pQj-al  this  point  do  not 

agree  with  the  measurements  made  by  the  seven-hole  pressure  probe. 

For  the  total  pressure  probe  aligned  with  the  freestream  flow,  there 
were  no  measurements  made  where  was  greater  than  0.  However, 

due  to  the  flow  angularity  at  this  point,  the  flow  has  probably 
separated  at  the  probe  tip  and  does  not  represent  a  measure  of  the 
true  local  total  pressure. 

For  the  total  pressure  probe  aligned  with  the  local  flow,  26 
percent  of  the  measurements  resulted  in  a  local  total  pressure  greater 
than  that  of  the  freestream  value.  These  measurements,  however,  are 
lower  than  those  obtained  using  the  seven-hole  probe.  It  should  be 
kept  in  mind,  however,  that  this  total  pressure  probe  was  bent  at  an 
angle  based  on  measurements  made  using  the  seven-hole  pressure  probe 
and  that  this  measurement  is  being  made  at  a  region  of  high  shear 
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flow.  If  the  shear  or  velocity  gradient  at  this  point  is  such  that 
holes  on  opposite  sides  of  the  seven-hole  pressure  probe  see 
significantly  different  velocities,  then  the  local  flow  angularity 
measurement  could  be  in  error  as  well.  If  this  measurement  of  local 
flow  angle  is  then  substantially  in  error,  the  flow  around  the  tip  of 
the  total  pressure  probe  supposedly  aligned  with  the  local  flow  could 
also  be  separated.  Separated  flow  around  the  tip  would  lead  to 
measurements  which  would  not  represent  the  true  local  total  pressure. 
Nevertheless,  the  average  value  for  based  on  these  measurements 

was  negative,  which  corresponds  to  the  theoretical  model  presented 
earlier.  The  maximum  positive  value  of  measured  with  this 

total  pressure  probe,  +0.040,  is  within  2.3  standard  deviations  of  the 
mean  value,  -0.013i  or  within  the  97.9  percent  uncertainty  band. 

Although  these  contradictions  do  not  exclude  the  possibility  that 
the  regions  of  positive  measured  by  the  seven-hole  pressure 

probe  do  exist,  they  seem  to  indicate  that  the  readings  may  be 
erroneous  in  regions  of  high  shear  flow.  It  must  be  recalled  that  at 
present  the  seven‘-hole  probe  is  calibrated  in  a  uniform  flow.  The 
regions  where  the  positive  values  of  have  been  measured  by  the 

seven-hole  pressure  probe  are  regions  of  high  shear  gradients  whose 
effect  may  not  be  properly  accounted  for  in  this  data  reduction  and 
whose  limits  need  to  be  defined.  These  results  do  not  imply,  however, 
that  the  seven-hole  pressure  probe  data  is  not  useful,  since  it  has 
more  than  demonstrated  its  value  in  mapping  hitherto  undefined 
external  flow  fields. 

V  I .  Summary 

The  external  flow  field  on  a  1/40th-scale  model  of  the  Northrop 


VATOL  configuration  was  mapped  using  the  seven-hole  pressure  probe 
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data  acquisition  system.  The  model  was  mounted  at  45  degrees  angle 
of  attack  relative  to  the  horizontal  reference  line,  and  the  test 
section  velocity  was  100  ft/sec.  Two  counter-rotating  vortices 
produced  by  the  leading  edge  extensions  were  mapped,  and  their  growth 
and  decay  were  observed.  Flow  interference  effects  due  to  the 
top'^mounted  inlets  were  identified  by  the  appearance  of  multiple 
pressure  contour  cells  and  a  low-energy  region  in  the  cores  of  the 
vortices  as  they  passed  over  the  inlets.  Flow  is  shown  to  be  attached 
from  the  upper  fuselage  to  the  inlet  locations  and  separated  over  the 
majority  of  the  upper  surfaces  of  the  aircraft  wings. 

Comparisons  of  seven-hole  pressure  probe  measurements  with 
hot-wire  anemometer  measurements  were  made  in  a  low  shear  flow  wake 
3  chord  lengths  aft  of  a  two-dimensional  airfoil  model.  The 
comparisons  showed  good  agreement,  with  the  seven-hole  pressure  probe 
measurements  of  local  velocity  being  within  3  •  *1  percent  at  a  test 
section  velocity  of  60  ft/sec  and  0.7  percent  at  a  test  section 
velocity  of  100  ft/sec  of  the  peak  velocity  deficit  measured  by  the 
hot-wire  anemometer.  The  overall  signature  of  the  airfoil  wake  was 
clearly  identified  by  both  measurement  systems. 

Finally,  the  measurements  of  local  total  pressures  made  by  the 
seven-hole  pressure  probe  were  compared  to  those  made  by  total 
pressure  probes.  Tests  were  conducted  using  the  1/40th-scale  Northrop 
VATOL  model  mounted  at  45  degrees  angle  of  attack  and  a  test  section 
velocity  of  100  ft/sec.  This  investigation  was  made  to  determine  the 
actual  values  of  local  total  pressures  in  flow  regions  where 
(according  to  seven-hole  probe  measurements)  these  values  are  greater 
than  the  freestream  total  pressure.  The  measurements  made  using  total 
pressure  probes  did  not  confirm  the  seven-hole  pressure  probe 
measurements.  To  Investigate  this  phenomenon  further,  research  into 
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seven-hole  pressure  probe  characteristics  in  regions  of  high  shear 
flow,  such  as  the  flow  typically  encountered  in  viscous  vortex  wakes 
is  recommended. 


English  Symbols 
A 

ALPHAT 

BETAT 

^Static 

'^Total 

d 

HRL 

M 

P 

r 

(  r  ,  'O 
t 

,U,,  ) 

(  u  ,  V  ,  w  ) 


V 

(  X  ,  y  ,  z  ) 
Greek  Symbols 


9  9 

r./ » wJ- 


i) 


Symbols 

9  T 

area,  inches-  or  feet- 
angle  of  attack,  degrees 
angle  of  sideslip,  degrees 

dimensionless  coefficient,  Pj  -  P„yP  -  P,„ 
dimensionless  coefficient,  P^,  -  P  /P^, 

?  eu 

distance,  inches  or  feet 
horizontal  reference  line 
Mach  number 
pressure,  psi  or  psf 
radius,  inches  or  feet 
polar  coordinates 
time,  seconds 

velocity  components,  polar  coordinates,  fps 
velocity  components,  Cartesian  coordinates, 
fps 

velocity,  fps 

Cartesian  coordinates,  inches  or  feet 

angle  of  attack,  degrees 

angle  of  sideslip,  degrees 
■> 

circulation,  ft^/sec 
k i nemat ic  viscosity,  ft"/see 
dens ity,  slug/ft- 
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t  s  t  a  n  d  a  r  d  d  e  V  i  a  t  i  o  n 

Subscripts 

t  localconditions 

o  total  or  stagnation  conditions 

“>  freestream  conditions 
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analytical  predictions  for  these  planforms. 

The  first  planform  configuration  (Figure  1)  is  ,)  tent  it  the 

midchord  so  that  the  trailing  half  simulates  a  trailing  'ia:  ;i-f!eoted 
downward  by  a  20-degree  angle.  The  angle  of  attack  of  thin  planform  is 
measured  with  respect  to  the  front  half  of  the  plate.  This  bent  plate 
configuration  is  interesting  because  of  the  effect  of  the  bend  on  the 
structure  of  the  tip  vortex.  The  discontinuity  in  slope  at  the  raidchord 
produces  a  secondary  vortex,  which  interreacts  with  the  main  vortex 
created  at  the  wing  tip.  In  fact,  when  we  observe  experiments  with  very 
low  Reynolds  numbers,  we  see.  that  this  secondary  vortex  core  becomes 
embedded  in  the  wake  wrapped  around  the  core  of  the  wing  tip's  main 
vortex.  This  secondary  vortex  core,  caused  by  the  midchord  bend, 
provides  a  pressure  characteristic,  or  signature,  in  the  main  vortex 
system  that  can  be  observed  at  downstream  locations.  The  rotation  rate 
of  the  main  vortex  can  then  be  determined  by  observing  the  change  in 
angular  position  of  this  pressure  signature  in  successive  locations 
(planes)  downstream  of  the  plate. 

The  second  planform  configuration  is  shown  in  Figure  2.  This 
highly-swept  delta  configuration  produces  a  leading  edge/tip  vortex 
system  similar  to  that  observed  around  the  strakes  of  modern  fighter 
aircraft.  This  vortex  system  builds  in  strength  when  moving  streamwise 
along  the  leading  edge.  Downstream,  it  wraps  the  trailing  edge  wake  up 
and  in  towards  the  centerline  of  the  model.  Since  modern  fighter 
configurations  use  this  strake  vortex  to  delay  flow  separation  over 
their  primary  wings,  it  is  of  critical  importance  to  know  the  location 


and  strength  of  this  type  of  1 i f t i ng -sur f ace  vortex  system. 
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1 1  .  Background 

Test  engineers  working  at  the  Aeronautic^  Laboratory  of  the  USAF 
Academy  have  developed  an  experimental  technique  that  can  measure 
pressure  and  velocity  information  from  steady  fields  even  with  highly 
irregular  flows.  The  test  technique  uses  a  multi-port  pressure  probe 
linked  to  a  computer  system  that  produces  real-time  flow  field  data  for 
any  geometric  location  in  which  the  pressure-probe  tip  can  be 
positioned.  This  particular  seven-hole  pressure  probe  was  first 
developed  to  measure  incompressible  flow  at  the  USAF  Academy  under  the 
sponsorship  of  the  NASA  Ames  Research  Laboratory  (Ref.  3).  The  probe's 
capabilities  were  later  expanded  to  include  supersonic  testing.  This 
probe  is  now  being  used  in  ongoing  research  to  document  the 
characteristics  of  1 i f t i ng -sur f ace  wakes,  particularly  those  created  by 
the  X-29  Forward-Swept  Wing  Flight  Demonstrator.  In  the  first  wake 
survey,  experimenters  used  the  seven-hole  probe  to  study  a  generic 
canard/swept-wing  configuration  (Ref.  5).  This  study  of  the  delta  and 
bent-plate  planforms  provides  background  information  on  1  i  f t i ng -sur f ace 
wakes  that  will  be  used  later  in  an  in-depth  investigation  of  the 
canard/wing  wake  interaction  of  a  1/10  scale  model  of  the  X-29 
Forward-Swept  Wing  Demonstrator  Aircraft. 

The  seven-hole  probe  used  for  these  wake  studies  is  shown  in  Figure 
3.  Because  the  probe  is  small  (.109  inches  in  diameter),  one  can 
position  it  inside  the  near  flow  fields  of  moderately-sized  wind  tunnel 
models  without  s i gn i f 1  can 1 1 y  disturbing  those  flow  fields.  The  position 
of  the  probe  is  maintained  by  an  automated  device  called  a  traverse, 
which  can  precisely  position  the  probe  anywhere  in  the  test  section. 
Since  this  positioning  traverse  is  controlled  by  the  wind  tunnel 
computer,  the  probe  location  can  be  loaded  into  the  computer  as  a  series 
of  coordinate  points. 
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The  probe  is  made  up  of  a  pressure  port  at  its  centerline 
surrounded  by  six  pressure  ports  on  its  conical  face.  The  seven 
pressure  ports  are  connected  to  electronic  pressure  transducers  by 
smal 1 -d i ame t er  tubing.  With  these  transducers  providing  electrical 
analogs  to  the  pressures,  a  computer  can  then  directly  supply  the 
digitized  values  of  raw  pressure  at  these  port  locations. 

Seven  sets  of  fourth-order  polynomial  equations  were  used  to  derive 
flow  field  pressures  and  velocities  from  the  raw  pressure  values 
measured  at  these  seven  pressure  ports.  These  equations,  along  with  the 
non-dime nsionalized  forms  of  their  coefficients,  are  presented  in  Ref. 

4.  Calibration  tests  at  known  flow  conditions  establish  the 
coefficients,  which  then  remain  stored  in  a  da t a -ac qu i s i t i on  computer 
for  use  with  unknown  flow  fields.  Using  the  raw  pressure  data  from  the 
seven-hole  probe,  these  equations,  with  their  calibrated  coefficients, 
provide  a  means  of  calculating  the  more  useful  values  of  total,  static, 
and  dynamic  pressures  and  velocity  vectors  at  any  selected  point  in  the 
flow.  The  pressure  and  velocity  values  for  each  data  point  are  recorded 
on  high-speed  disk  storage  for  later  use. 

The  data  points  for  the  tests  are  arranged  in  planes  having  common 
streamwise  coordinates.  This  procedure  makes  it  easier  both  to  acquire 
and  present  data.  The  planes  are  chosen  at  various  streamwise  locations 
as  measured  from  the  trailing  edge  of  each  model.  By  locating  the  wake 
characteristics  at  each  streamwise  position  in  the  flow,  one  can  observe 
and  interpret  their  development  as  the  flow  proceeds. 

With  a  complete  set  of  pressure  and  velocity  data  stored  for  each 
data  point,  one  can  quickly  generate  veloc i ty- vector  and  pressure- 
contour  plots  in  the  data  planes  as  testing  proceeds.  This  insures  that 
all  of  the  important  locations  in  the  flow  field  are  surveyed  in  each 
data  plane.  Later,  these  plots  can  be  used  to  locate  important  flow 
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characteristics  such  as  vorticity,  wakes,  and  downwash  for  final 
interpretation  of  the  lifting-surface  wakes. 

Figure  4  is  a  sample  cr oss^ vel oc  i  ty ,  flow-field  data  plot  generated 
with  data  taken  for  the  bent  plate  at  20  degrees  angle  of  attack.  This 
figure  shows  the  cross-velocity  vectors  at  the  data  points  surveyed  on 
the  data  plane  immediately  downstream  from  the  trailing  edge  of  the  bent 
plate.  These  vectors  represent  the  projection  on  the  data  planes  of  the 
total  velocity  vector  measured  at  each  probe  location.  The  origin  of 
each  vector  gives  the  geometric  location  of  the  data  point;  its  length 
gives  the  relative  magnitude  of  the  in-plane  velocity  component;  and  the 
arrow  head  indicates  the  flow  direction. 

In  Figure  “4  the  view  is  looking  upstream  at  the  port-side  wing. 

The  cr OSS- ve loc i ty  vectors  indicate  a  strong  clockwise  swirling  motion 
(typical  of  a  strearawise  vortex)  centered  about  1.6  inches  in  the 
negative  x  direction.  The  direction  of  the  rotation  in  the  cress  flow 
is  that  expected  about  the  tip  vortex  generated  by  positive  lift  on  the 
port  wing.  The  center  of  this  motion  will  then  be  the  core  filament  of 
the  vortex.  From  these  velocity  fields,  we  can  determine  the  locations 
and  relative  strengths  of  vortex  systems,  as  well  as  the  direction  of 
the  downwash  flow  field. 

Pressure  contour  plots  complement  those  containing  the  cross- 
velocity  vectors  and  can  be  most  useful  in  wake-location  calculations. 
Figure  5  shows  the  total-pressure  coefficient  contours  for  the  same  data 
points  viewed  in  the  same  direction  as  in  Figure  4.  This  total-pre;-, sure 
coefficient  is  defined  by  total  pressure  comparisons  between  local  and 
freestream  conditions. 
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The  difference  between  local  (P  )  and  freestream  (P  )  total 

Local  Tunnel 

pressure  is  nond  i  mens  1  ona  1 1  zed  by  the  freestream  dynamic  pressure.  Cj, 

o 

becomes  negative  as  irreversible  losses  in  total  pressure  occur  locally, 
because  of  viscosity  in  regions  having  large  velocity  gradients.  These 
regions  occur  near  the  core  of  vortex  filaments  and  in  the  boundary 
layers  of  lifting  surfaces.  In  downstream  flow,  this  loss  in  total 
pressure  provides  a  signature,  which  identifies  the  freestream  wake 
locations. 

In  Figure  5  contours  of  various  values  of  Cp  indicate  regions 

o 

where  losses  in  total  pressure  have  occurred.  Note  that  the  wake  of  the 
plate  is  apparent  in  the  pressure  values.  Downstream  from  the  tip  chord 
of  the  plate,  this  wake  wraps  around  the  tip-vortex  core.  Using  contour 
plots  like  those  in  Figure  5,  we  can  determine  the  wake  locations  and 
the  vortex  systems  (at  their  edges)  for  both  the  bent  plate  and  delta 
planforms  in  each  plane.  The  plots  of  their  wake  locations  appear  at 
the  end  of  this  report. 

III.  Test  Apparatus 

We  tested  two  planforms  in  the  subsonic,  con t i nuous- f 1 ow  wind 
tunnel  (described  in  Ref,  5)  in  the  Aeronautics  Laboratory  at  the  USAF 
Academy,  The  freestream  velocity  was  maintained  at  100  ft/sec.  The 
bent-plate  model  of  Figure  1  was  mounted  on  a  thin  vertical  strut  from 
the  tunnel  celling.  To  develop  the  following  data,  each  model  was  set 
at  two  angles  of  attack;  10  and  20  degrees  with  respect  to  the 
freestream  direction  (as  measured  from  the  forward  half  of  the  model). 
The  unit  normals  of  the  planes  defining  the  data  point  streamwise 
locations  were  always  parallel  to  the  freesteara  direction. 

The  delta  planform  model  was  also  run  at  10  and  20  degree 
angles  of  attack.  Again,  the  unit  normals  of  the  planes  used  to 
organize  the  data  points  were  parallel  to  the  freestream  direction. 
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I  V .  Results 

The  data  obtained  for  the  two  planforms  is  presented  graphically. 

A  cross-velocity  plot  appears  for  each  set  of  data  points,  followed  by 
a  second  plot  showing  the  lifting  sui’face  wake.  The  wake  is  plotted 
from  the  centerline  outward,  showing  the  way  it  becomes  wrapped  in  the 
tip  vorticity  until  its  pressure  signature  cannot  be  distinguished  from 
the  inner  field  of  the  tip  vorticity.  The  core  location  of  the  tip 
vortex  is  also  noted.  If  the  data  plane  intersects  the  wing  planform, 
this  too  is  noted.  The  steamwlse  development  of  the  wake  for  different 
wing  angles  of  attack  can  be  compared  for  both  planforms. 

A.  Bent-Plate  Results 

The  bent  plate  data  is  taken  from  the  centerline  towards  the 
port  side  of  the  plate.  Since  all  flow  configurations  are  symmetric 
about  a  vertical  plane  passing  through  the  planform  centerline,  the  flow 
field  of  only  one  side  of  the  planform  is  presented  to  allow  greater 
detail  in  the  plots.  The  coordinate  system  used  for  all  the  data  points 
has  its  origin  at  the  intersection  of  the  vertical  plane  of  symmetry  and 
the  trailing  edge  of  the  plate.  The  vertical  direction  (perpendicular 
to  the  freesteam  velocity  vector)  up  (the  positive  lift  side)  from  the 
plate  is  the  positive  y  coordinate.  The  lateral  direction 
(perpendicular  to  both  the  freestream  direction  and  the  y  direction)  is 
positive  towards  the  starboard  wingtlp  and  is  designated  the  positive 
X  direction.  All  the  plots  are  presented  in  such  a  way  that  the 
freesteam  velocity  vector  is  out  of  the  data  plane.  Thus  the  plots  are 
to  be  viewed  as  looking  upstream. 

The  data  planes  for  10  degrees  angle  of  attack  are  at  streamwise 
locations  relative  to  the  trailing  edge  origin  of  -2,  .2,  2,  4,  6,  9, 


and  12  inches.  The  data  planes  for  20  degrees  angle  of  attack  are  at 
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locations  of  -2,  .2,  2,  6,  and  11.5  inches.  The  data  is  presented  in 

progressive  steamwise  locations.  The  cross-velocity  and  wake-location 
plots  are  shown  for  each  plane. 

Figures  6  through  19  show  the  location  of  the  wake,  the  velocity 
field,  and  the  vortex  core  for  the  bent  plate  at  10  degrees  angle  of 
attack.  Figures  20  through  31  show  the  same  at  20  degrees  angle  of 
attack.  For  both  angles  of  attack  and  any  streamwise  location,  the  x 
location  for  the  vortex  core  remains  approximately  -1.5  inches.  There 
are,  however,  changes  in  the  y  location  of  this  vortex  core. 

The  tip  vortex  core  location  in  the  y  direction  is  influenced  by 
the  direction  and  strength  of  the  downwash  from  the  lifting  surface. 

The  downwash  field  —  and  therefore  the  movement  of  the  vortex  core  in 
the  downwash  direction  —  is  stronger  at  20  degrees  angle  of  attack  than 
at  10  degrees.  This  is  evident  vfhen  one  compares  the  locations  of  the 
tip  vortex  core  in  Figure  15  and  29  or  Figures  19  and  31. 

B.  Delta  Planform  Results 

The  results  of  the  low  aspect-ratio  delta  wing  planform  tests 
are  presented  in  Figures  32  through  55.  Again,  the  data  points  lie 
above  and  in  the  wake  of  the  port  side  of  the  wing.  The  same  coordinate 
system  is  used  for  the  delta  planform  as  was  used  for  the  bent  plate, 
with  the  origin  at  the  centerline  trailing  edge  location  of  the  delta 
planform. 

The  data  planes  for  10  degrees  angle  of  attack  are  at  streamwise 

locations  relative  to  the  trailing  edge  origin  of  -2,  .2,  2,  i4  ,  6,  9, 

and  12  inches.  The  data  planes  for  20  degrees  angle  of  attack  are  at 

I  streamwise  locations  of  -2,  .2,  2,  *1,  and  6  inches.  The  plots  of  the 

t 

,  wake  locations  and  cross -vel oc 1 ty  vectors  for  each  plane  of  data  points 


are  presented  in  progressive  streamwise  fashion  to  indicate  the  movement 
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of  the  wake  as  it  travels  downstream  from  the  delta  wing. 

The  delta  wing  results  at  10  degrees  angle  of  attack  are  contained 
in  Figures  32  through  45.  The  results  at  20  degrees  angle  of  attack 
appear  in  Figure  46  through  55.  The  highly-swept  leading  edge  develops 
a  strake-like  vortex  system  that  travels  in  both  the  spanwlse  and 
downwash  directions  as  progressively  downstream  data  planes  are 
surveyed.  The  downwash  velocity  vectors  can  be  seen  in  the  velocity 
vector  plot  for  each  data  plane.  Note  the  clockwise  rotation  of  the 
port~side,  leading-edge  vortex  system.  This  is  consistent  with  the 
positive  lift  condition  created  at  10  degrees  angle  of  attack.  Similar 
trends,  only  stronger,  are  found  at  20  degrees. 

The  tip  vortex  core  appears  to  seek  the  same  location  at  either  10 
or  20  degrees  angle  of  attack:  approximately  -2.3  inches  on  the  x  axis 
and  .2  inches  on  the  y  axis.  This  location  for  the  tip  vortex  core 
appears  to  have  been  achieved  for  both  angles  of  attack  in  the  data 
plane  6  inches  downstream  of  the  trailing  edge. 

The  shape  of  the  delta  wing's  downsteam  wake  differs  from  that  of 
the  bent  plate.  Downstream  of  the  delta  wing,  the  portion  of  the  wake 
that  becomes  wrapped  around  the  vortex  core  is  wrapped  in  a  looser 
spiral  than  that  found  behind  the  bent  plate.  This  is  especially 
noticeable  in  the  x  locations  of  the  outboard  wake  that  is  wrapping 
around  the  core. 


Borit  Plnte  Cross  Ve 


Bent  Plate  Cross  Veloc 
Vector  Plot,  ‘  =  10°  a 


L'SAFA-TR-84-5 


/  /  /  /  /  I  J  J I 
f  / 1 1 1  i}X^ 
/Min 

I  1 M  \ 

\  \  \  \  N  's, 

N  ^  V  N  •S 
^  N  V  V  ^ 

•  N  K  “s.  --s  ^  ^  ^  »— ^— 
.  V  •s  •».  ^  w— 


N  V\  V  i  n  \  "  '•I 
rr/v  /  i  M  '  ^ 

/  /  /^  /  /  I  I  ‘  M 

7  /  V  /  /  /  t  I  V 

^  /  /  /  '  t 

y  //  I  i  I 

^  y  //  I  i  I 

‘ ^  ^  i  i  I 

^  ^  ^  /  i  t 

^  ^  ^  y  /  / 

^  ^  ^  ^  y'  y  ^  ✓ 


o  w 

V-  c 

o  -  ? 

AJ  O 


Qj  O  O 


c 

o  t— 

Cl 

Cl 

tfl 

>  O' 

T  n  r  h  C’ 


WAKE— A 


Plate  Cross  Velocity  Figure  19.  Bent  Plate  Wake  location 

r  Plot,  I  =  10°  and  t  =  10°  and  12  Inches 


Plot,<=20°and  i=20°and2Inchcs  Upstream 


Velocity  Figure  27.  Bent  Plate 


nt  Plate  Wake  Location 
=  20°  and  11.5  Inches 


eltn  Wing  Cross  Velo 
ector  Plot,  '*=  10° 


c  he  s 


Wing  Cross  Ve 
Plot ,  I  =  10 


Delta  Wing  Cross  Velocity  Figure  49.  Delta  Wing  Wake  Location 

Vector  Plot,  ‘  =  20°  and  '  =  20°  and  .2  Inches 


, / / / / in^79v,i} 

, ,-,  , , , 


i  ‘ 

/  /  ''/It 


\  \  \  \'\ 

^  \  \  \  N  •s,»-v»' 

^  S  N  N 

\  s  V  N  V,  ^ 


y/  [  '  /  /  ! 

^  /'  ^  /  /  i  f 

^  ^  //  / 

^  ^  Z'  /  / 

Z*  z'  z 


r  /^  / 


Vector  Plot,  ‘  =  20 
2  Inches  Downstream 


52.  Delta  Wing  Cross  Velocity  Figure  53.  Delta  Wing  Wake  Loco 

Vector  Plot,  t  =  20°  and  ^  =  20°  and  4  Inches 

4  Inches  Downstream  Downstream 


L'  S  A  1'  A  -  r  R  -  «  4  -  5 


V .  ConcI usions 

The  low  aspect-ratio  delta  wing  demonstrates  wing  wake 
sensitivities  to  leading-edge  sweep.  The  leadihg-edge  sweep  of  the 
planform  s i gn i f i ca n t 1 y  affects  both  the  way  the  wake  becomes  wrapped 
around  the  lifting-surface,  tip-vortex  system  and  the  movement  of  the 
vortex  system  downstream.  With  large  leading-edge  sweep,  the  wake 
bulges  in  a  more  spanwise  fashion.  Also,  the  tip  vortex  moves  spanwise 
due  to  sweep  effects  but  tends  to  go  to  approximately  the  same 
downstream  location  relative  to  the  delta  wing  trailing  edge  provided 
the  angle  of  attack  is  not  too  large. 

The  bent-plate  configuration  addresses  camber  effects  on  wake 
characteristics  using  the  20-degree  bend  at  midchord.  The  secondary 
vortex  system  caused  by  the  discontinuity  In  slope  at  the  midchord  of 
the  bent  plate  could  not  be  discerned  in  the  pressure  data.  This  is  due 
to  the  turbulence  of  the  Interior  vortex  flow  field.  The  wake  strength 
tended  to  determine  the  downstream  vortex  core  location  within  the 
downstream  range  available  in  this  test.  Very  little  outward  drift  of 
the  vortex  core  was  noted  regardless  of  plate  angle  of  attack. 

The  actual  pressure  and  velocity  values  for  each  data  point  in 
these  test  cases  (recorded  on  magnetic  tape)  can  be  obtained  by 
contacting  the  Director  of  the  Aeronautics  Laboratory,  Department  of 
Aeronautics,  'JSAF  Academy. 
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SECTION  II 


Aeronautical  History 


U  S  AFA-  TK- rt-'4  -  T 


T!iK  KAKFY  FAYF  ;F  A  F  h  0  N  A  ! '  T  I  F  A 
i  r  J4‘;;rrr'(^y  I.  Triylor' 

K  (1  i  t  o  r-  *  s  N  o  t  (' 

T  h  <?  !'  ■ lowiM*;  artiolp  h  a  r.  F  pt-n  r  t-p  r  i  n  t  <■  cj  ft-om  the  Journal  of  the 
Royal  A  e  r  o  r.  a  u  t  i  c  -  4  1  F  o  c  i  e  t  y  Volume  70,  Number  661,  J  a  n  u  a  r  y  1  966,  with 
the  p e  r’ m  i  s  3  i  ( ) n  of  the  Royal  A  e r-  o n a  u  t  i  c a  I  Society.  The  author,  Sir  G  .  I  . 
Taylor,  was  a  prolific  researcher  who,  in  his  60  plus  years  as  a 
scienti;it  and  academician,  published  more  than  200  papers  on  a  variety 
of  subjects  including  aerodynamics,  hydrodynamics,  meteorology,  and 
mechanics  of  solids.  His  career  spanned  the  development  of  aeronautics 
from  the  first  days  of  heavier-than-air  aircraft  to  the  days  of  manned 
spaceflight.  Hi.-,  expertise  and  breadth  of  experience  combine  to  make 
his  recollections  of  the  early  days  of  aeronautics  a  very  interesting 
piece  of  reading.  Most  of  Taylor's  published  works  have  been  collected 
in  the  four-volume  set  The  Scientific  Papers  of  G.  I.  Taylor,  edited  by 
G.  K.  Batchelor  and  published  by  the  Cambridge  University  Press  between 
1958  and  1971.  Anyone  looking  through  these  volumes  can  only  marvel  at 
the  quality  and  quantity  of  research  produced  by  this  remarkable  man. 

It  is  always  rather  risky  to  ask  old  people  to  talk  about  what 
things  were  like  when  they  were  young.  Having  turned  on  the  tap  you  may 
not  see  how  to  turn  it  off  again.  Fortunately  there  is  a  circumstance 
which  keeps  this  trouble  under  control.  Old  people  remember  a  lot  of 
things  when  there  is  no  call  for  them  to  do  so  but  seldom  when  there  is. 
That  is  the  reason  why  I  have  written  down  most  of  what  I  can  say  for 
this  lecture . 

I  cannot  attempt  to  give  a  complete  or  an  objective  picture  of 
aeronautical  science  when  I  started  thinking  about  it.  The  best  I  can 
do  is  to  describe  a  few  of  the  things  which  happened  to  me  in  those  days 
to  try  to  convey  some  pictures  of  our  primitive  ideas. 

My  earliest  aerodynamic  experiences  were  all  connected  with  sails. 

I  built  a  boat  in  our  home  in  St  John’s  Wood  while  I  was  at  school.  It 
was  before  the  days  when  woodworking  firms  would  s\jpply  ready-made  kits 
so  I  had  to  do  all  the  design  and  construction  work  myself.  Fig.  1  was 
taken  while  the  sails  my  mother  had  made  were  being  fitted.  Those  of 
you  who  know  about  boats  will  see  it  is  old  fashioned,  tiut  it  was  quite 
fast.  After  building  the  hull  in  my  bedroom  on  the  second  floor  we  had 
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to  get  it  through  the  window.  It  was  a  sa.y.h  window  which  ecu  id  tak'cc 
out.  We  turned  the  boat  on  it.s  .side  becau.se  it  was  wider  than  the 
win<!ow  frame.  We  put  slings  round  it  and  attached  them  to  a  rope  wiiich 
we  slung  over  a  borrowed  builder's  pulley  at  the  end  of  a  joist 
projecting  from  my  parents'  bedroom  on  the  top  floor.  After  making  the 
lower  end  of  the  rope  fast  to  the  garden  roller  we  got  the  boat  out  of 
the  window  and  left  it  hanging  while  we  went  down  into  the  garden.  We 
then  detached  the  rope  from  the  garden  roller  and  my  brother,  my  cou.sin, 
and  I  started  to  let  it  out.  As  we  were  doing  this  the  boat  above  our 
heads  slipped  slightly  in  the  slings.  My  brother  let  go  to  run  upstairs 
to  refix  it  and  we  then  discovered  that  the  boat  was  a  little  heavier 
than  my  cousin  and  me.  Fortunately  we  both  held  on,  no  further  slipping 
occurred  and  ray  brother  caught  us  almost  before  we  had  left  the  ground. 
It  was  a  good  lesson  in  practical  mechanics. 


Figure  1.  The  boat  in  the  garden 
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Aftor  iHurichirif:  rhe  Scat  at  il  amme  r  sm  i  t  h  I  sailr-A  a.r,w:,  t  R.  o  a  I  h  *■  r:  •. 


and  a  c  r  (.1  s  s  to  o  c  h . 


o  a  c  h  a  i  d  e  o  f  t  h  i 


about  tho  ba!  aru 


ling  vt-ry  uncomfortably  with  o  r:  a  1 


d  o  I'i o  ,  T  d  i  f!  not  think  rt  u  c  h  i  n  tho  a 


a  ti  d  h  y  d  r  o  "  d  y  n  a  m  i  c  fore  o  s  b  1  w  e  e  n  n  a  i  I  r,  and 


hull.  It  waa  ('nougi!  to  know  th.  at  a  ails  push  one  along,  but  ‘bat  voyag*' 
taught  me  a  little  about  tlio  ways  of  seamen.  On  the  I'eturn  journey  up 
tho  Thames  I  met  a  ti-'esome  sea  and  a  head  wind  in  Sea  Heach.  Thcia; 
were  some  ‘■•0  or  6  0  sailing  barges  beating  up  at  the  time  and  I  ma  n  a  gr-d 
to  hook  onto  one  of  them  and  get  a  tow  for  the  boat  and  .a  passage  for 
me.  Somewhere  off  Woolwich  my  barge  was  involved  in  a  very  interest  int’ 
triple  collision  with  two  others,  each  barge  hitting  both  of  t  h ci  o  t  h e  r  . 
One  of  them  ultimately  sank.  The  verbal  exchanges  which  passed  between 
the  three  .skippers  and  their  mates  when  we  interlocked,  were  an 
eye-opener  to  me.  I  listened  with  admiration,  never  thinking  that  an 
analogous  experience  might  one  day  call  forth  an  unexpected,  although  of 
course  far  smaller  talent  in  me.  I  was  taking  my  19-ton  cutter  before  a 
strong  wind  into  the  sea  lock  at  the  end  of  the  Crinan  Canal  when  my 
crew  failed  to  belay  In  time  the  rope  which  should  stop  the  way.  I 
thought  we  were  going  to  ram  the  lockgates  and  in  the  heat  of  the  moment 
said  a  few  words  about  the  situation.  I  was  later  told  that  a  summer 

visitor  who  wa.s  standing  on  the  lock  side  watching  this  manoe  was 

seen  to  pick  up  her  little  P e k e  ,  and  stop  up  its  ears. 

The  triple  collision  was  also  educative  in  another  way.  All  tlv' 

three  skippers  gave  different  accounts  of  the  cause  of  the  accident. 

The  skipper  of  the  barge  called  the  "Jolly  Dogs"  was  certainly  in  the 
w r n g  ,  hut  he  put  the  blame  on  my  skipper  who,  fortunately  for  m f’  ,  was 
entirely  correct  in  his  actions.  A  pretty  ca.se  was  brewing  for  decision 
in  the  law  courts  and  the  counsel  for  my  skipper  asked  me  to  g '!  to  i i  i 
office  and  d  e  s,  c  r  i  b  e  what  happened.  I  w  a  .s  able  to  give  a  c  o  n  s  i  s  1  e  n  t. 
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picture  of  the  seQuence  of  events  and  air  reed  to  e;  i  vo  v  i  i  r.  tht 

court.  A  short  time  afterwards  the  lawyer  wrote  that  this  ret  ' 

necessary  because  when  he  told  the  opposing  lawyers  that  hr-  rai  a  truiv 
independent  witness  who  could  describe  the  collision,  tl;ey  hriew  it  was 
useless  to  go  on,  even  though  the  witness  was  only  a  s  c  h  o  o  1  s  v  ,  arc  tr'^v 
withdrew  their  allegations. 

My  first  introduction  to  fluid  mechanics  as  a  science  ca:r--  throuich 
meteorology.  After  .studying  the  subject  in  its  theoretical  asp':;:!  f  r'  a 
time,  I  was  appointed  meteorologist  to  the  old  wooden  whaling  si;  ip 
Scotia  which  was  chartered  in  1913  by  the  government  and  t  ii  e  shipping 
companies  to  look  for  icebergs  in  the  North  Atlantic  and  report  their 
positions  by  means  of  the  newly  invented  m, ethod  of  wirele.ss  telegraphy. 
You  may  remember  that  the  T 1  '  a  n i c  went  down  in  1912  through  striking  an 

iceberg.  The  Scot  i  a  was  a  wooden  barque  of  238  tons.  Fig.  2  sh,  ows  her 

quarterdeck.  I  took  this  picture  one  sunny  morning  when  we  were  stuck 

in  ice  firm  enough  for  me  to  stand  on.  The  Scotia  had  been  u  .s  e  d  b  y  t  h  e 

explorer  D  r  u  c  o  in  the  Antarctic  and  was  equipped  with  an  auxiliary  steam 
engine  of  almost  unbelievable  inefficiency  which  would  d ^ i v e  her  up  to 
4  1/2  knots  in  a  calm,  alLhougi'.  would  make  about,  n  knots  under  sail 

in  a  gale.  She  had  two  whale  boat.a  which  were;  seldom,  if  ev“r,  used 
on  our  expedition.  They  were  the  original  ho.ats  u.sed  whi'i;  w  ;i  .a i  •' g  was 
carried  Gut  by  rowing  up  to  a  wh.ale  and  .otiokirig  a  ii  a  r' no  u' n  int'i  it. 

They  were  very  much  in  my  way  wh«n  T  flew  kit'"-:'  ■' .a  r  t' "y  i  n  g  'r  e  t  <' o  r  !"■  j  c  g  i  e  ; 

1  n  s  t  r  u  m e  n  t  .0  from  h  *> , ■  q  u  a  r  t  !-■  r  d  i-  n  y  .  A  n  c  t  h  e  r  o  b  u  *■  i  ■  n  "w  .a  o.  t  "ri  e  , ■ 

boom.  I  noticed  th.at  ti".  o  nkiryo'  r-  no-.-en  \  '  h--  r  ]■■  ■  s- s.  ,i  ;  i  •  i 

him  why.  H.o  .said  the  bO':rri  war,  q-iito  rotte-i.  u  t  wh'"-!  i  .  ■  ■ 

kept  it  in  p  o  n  ■  t  i  o  n  t;  e  war  h  <:■  r  r  i  f  i  ■  ■  i  n  1  ;■  a  i  i  ,  ”  W  1.  c  _  .• 

my  spanker  room'.'  t!  n  t  h  i  n  k  a  t  "  Afti-r 

2  u  n  d  e  e  while  tri'-  e^^'W  rot  .-r  r  .-.■a o 
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way.  Most  of  our  time  was  spent  on  the  banks  of  Newfoundland  where  the 
i ceberg- 1 aden  cold  water  comes  south  from  Bafflns  Bay.  Although  the  air 
at  sea  level  was  near  freezing  point  the  air  a  few  hundred  feet  up  had 
frequently  come  straight  off  the  heated  land  of  the  American  continent 
and  was  so  hot  that  my  kites  felt  warm  to  the  touch  when  I  reeled  them 
in.  It  was  thinking  about  the  way  heat  transferred  to  and  from  the 
earth's  surface  that  roused  my  interest  in  turbulence. 


Figure  2.  Quarter  deck  of  the  Scotia 

While  I  was  writing  up  the  results  of  the  Scotia  expedition,  war 
broke  out  on  ilth  August  1914.  The  only  way  I  could  think  of  in  which 
my  experience  might  be  of  value  was  as  a  meteorologist  so  I  went  to  the 
War  Office  and  asked  whether  I  could  give  weather  forecasts  for  military 
operations  to  the  expeditionary  force  then  being  organised  if 
telegraphic  facilities  could  be  made  available.  The  officer  whom  I  saw 
there  did  not  seem  to  doubt  that  I  could  give  accurate  forecasts  “ 


which  he  might  very  well  have  done 


but  said  that  they  would  be  of  no 
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picture  of  the  sequence  of  events  and  agreed  to  give  evidence  in  the 
court.  A  short  time  afterwards  the  lawyer  wrote  that  this  would  not  be 
necessary  because  when  he  told  the  opposing  lawyers  that  he  had  a  truly 
independent  witness  who  could  describe  the  collision,  they  knew  it  was 
useless  to  go  on,  even  though  the  witness  was  only  a  schoolboy,  and  they 
withdrew  their  allegations. 

My  first  introduction  to  fluid  mechanics  as  a  science  came  through 
meteorology.  After  studying  the  subject  in  its  theoretical  aspect  for  a 
time,  I  was  appointed  meteorologist  to  the  old  wooden  whaling  ship 
Scotia  which  was  chartered  in  1913  by  the  government  and  the  shipping 
companies  to  look  for  icebergs  in  the  North  Atlantic  and  report  their 
positions  by  means  of  the  newly  invented  method  of  wireless  telegraphy. 

You  may  remember  that  the  Titanic  went  down  in  1912  through  striking  an 

iceberg.  The  Scotia  was  a  wooden  barque  of  238  tons.  Fig.  2  shows  her 

quarterdeck.  I  took  this  picture  one  sunny  morning  when  we  were  stuck 

in  ice  firm  enough  for  me  to  stand  on.  The  Scotia  had  been  used  by  the 

explorer  Bruce  in  the  Antarctic  and  was  equipped  with  an  auxiliary  steam 
engine  of  almost  unbelievable  inefficiency  which  would  drive  her  up  to 
4  1/2  knots  in  a  calm,  although  she  would  make  about  9  knots  under  sail 
in  a  gale.  She  had  two  whale  boats  which  were  seldom,  if  ever,  used 
on  our  expedition.  They  were  the  original  boats  used  when  whaling  was 
carried  out  by  rowing  up  to  a  whale  and  sticking  a  harpoon  into  it. 

They  were  very  much  in  my  way  when  I  flew  kites  carrying  meteorological 
instruments  from  h:;r  quarterdeck.  Another  obstruction  was  the  spanker 
boom.  I  noticed  that  the  skipper  never  set  the  spanker  sail  and  I  asked 
him  why.  He  said  the  boom  was  quite  rotten,  but  when  I  asked  him  why  he 
kept  it  in  position  he  was  horrified  and  said,  "What,  go  to  sea  without 
my  spanker  boom?  Unthinkable!"  After  spending  24  hours  at  anchor  off 


Dundee  while  the  crew  got  sober  enough  to  work  the  ship,  we  got  under 
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value  in  military  operations  even  if  accurate.  "Soldiers  don't  go  into 
battle  under  umbrellas  “*  they  go  whether  it  is  raining  or  not."  I 
expect  he  was  politely  telling  me  that  the  army  had  all  the  science  it 
needed,  but  when  six  months  later  our  people  found  that  the  French  army 
derived  great  benefit  from  their  meteorological  service  our 
meteorological  office  was  asked  to  organise  one;  and  one  of  its  staff, 
Colonel  Gold,  did  so  very  successf ul 1 y . 

Returning  to  the  War  Office,  the  officer  asked  me  to  leave  my  name 
and  address  in  case  anything  turned  up  for  which  a  scientist  could  be 
used.  I  went  home,  I  must  confess,  rather  disappointed,  for  I  had  no 
desire  to  become  a  dispatch  rider.  That  evening  a  messenger  came  from 
the  War  Office  saying  that  if  I  went  down  to  the  Royal  Aircraft  Factory 
at  Farnborough  I  would  be  employed  there  as  a  scientist.  I  went  there 
on  the  5th  or  6th  August.  What  had  happened  in  the  War  Office  was  that 
Major  Sefton  Brancker,  then  I  think  head  of  the  flying  wing  of  the  Army 
had  happened  to  go  into  the  room  of  the  man  with  whom  I  had  just  had  an 
interview.  This  officer  had  told  Brancker  of  my  visit  and  asked  him 
whether  he  could  find  any  use  for  academic  scientists  in  his  branch. 
Brancker  said  he  could  use  as  many  as  he  could  lay  his  hands  on.  He 
took  my  name  and  sent  his  messenger  at  once. 

On  arriving  at  Farnborough  O'Gorman,  who  was  then  the  director. 
Introduced  me  to  a  fascinating  character  called  Teddy  Busk  who  asked  me 
to  join  his  mess  in  Arnold  House.  This  mess  was  run  chiefly  by  members 
of  the  Naval  airship  section  among  whom  was  Cave-Brown-Cave  who 
afterwards  became  Professor  of  Engineering  at  Southampton  University. 

At  first  only  Busk  and  I  and  F.  M.  Green  were  not  in  the  Navy,  but 
later  we  were  joined  by  other  civilians,  F.  W.  Aston,  the  inventor  of 
the  mass  spectroscope  and  the  principal  discoverer  of  isotopes,  and 
Lindemann,  afterwards  Lord  Cherwell.  The  mess  subsequently  moved  to  a 
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house  called  "Chudleigh":  I  cannot  remember  who  were  members  at  the  time 
of  the  move,  but  Farren,  Glauert,  R.  H.  Fowler,  George  Thomson,  E.  D. 
Adrian  and  Melvill  Jones  were  certainly  "Chudleighites." 

When  I  arrived  at  Farnborough  I  was  put  to  work  in  H  department 
onto  which  the  odd  jobs  that  did  not  fit  in  with  any  other  branch  were 
decanted.  It  was  housed  In  a  building  with  a  tin  roof  which  looked  out 
over  the  airfield.  The  aeroplanes  used  to  come  in  close  to  the  roof 
making  a  fearful  noise  and  quite  often  crashed  just  outside  our  windows. 

Some  of  the  things  we  were  set  to  do  now  seem  futile  in  the  light 
of  hindsight,  and  some  even  then  seemed  rather  wild,  but  of  course  we 
were  under  orders  from  the  Army  and  had  to  do  what  we  were  told. 
O’Gorraan,  however,  gave  us  a  pretty  free  hand  and  we  were  able  to  turn 
down  quite  a  number  of  things  which  could  not  possibly  have  worked.  One 
job  which  was  put  on  us,  although  we  never  thought  it  would  be  useful, 
was  the  designing  of  darts  to  be  dropped  on  troops  from  aircraft. 

Melvill  Jones  and  I  were  the  principal  people  concerned  and  both  he  and 
I  have  described  the  sequence  of  experiments  Involved,  so  I  will  only 
give  a  brief  account  now.  The  French  had  used  pointed  steel 
"flechettes"  but  we  found  that  they  whirled  when  dropped,  even  though 
the  centre  of  pressure  when  falling  obliquely  was  well  behind  the  centre 
of  gravity.  In  any  case  it  seemed  to  be  a  matter  of  prestige  that  we 
must  have  English  darts  so  we  designed  and  tested  some  with  the 
necessary  aerodynamic  properties.  Our  last  test  was  to  find  their 
distribution  on  the  ground  when  dropped  from  an  aeroplane  out  of  a 
container,  and  we  got  a  pilot  to  drop  them  onto  an  unused  piece  of 
ground.  To  measure  their  distribution  we  searched  the  ground  and 
whenever  we  found  the  rear  end  of  a  dart  projecting  from  it  we  pushed 
a  square  of  paper  over  it  to  mark  its  position.  We  had  just  completed 
this  work  and  were  preparing  to  photograph  the  field  when  a  cavalry 
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officer  came  up  and  asked  what  we  were  doing.  When  we  expliined  tf-.at 
the  darts  which  covered  the  field  had  been  dropped  from  an  aeroplane  he 
looked  at  them,  seeing  each  piercing  a  white  square  and  said,  in  a 
surprised  tone,  "If  I  had  not  seen  this  I  would  never  have  believed  that 
it  was  possible  to  make  such  good  shooting  from  the  air."  Having 
finished  that  job,  the  darts  were  never  used,  not  apparently  for  the 
reason  we  had  originally  raised  as  an  objection,  namely  that  missiles 
descending  vertically  would  have  a  much  smaller  chance  of  hitting  a 
standing  man  than  the  same  weight  of  bomb  going  off  on  contact  with  the 
ground  and  throwing  its  fragments  horizontally.  The  reason  we  were 
given  was  that  darts  were  regarded  as  inhuman  and  could  not  be  used  by 
gentlemen . 

All  the  things  we  did  were  not  so  absurd  as  our  work  on  darts  and 
by  way  of  contrast  I  should  like  to  trace  one  set  of  ideas  which 
originated  at  Farnborough  at  that  time.  One  of  the  troubles  which  we 
were  asked  to  look  at  was  the  weakening  of  propeller  shafts  for 
transmitting  power  owing  to  key-ways  cut  in  them.  About  this  time  a 
young  man  called  A.  A.  Griffith  arrived  at  the  factory.  He  was 
obviously  very  able  and  he  came  and  worked  with  me  on  the  subject.  We 
devised  what  is  now  called  an  analogue  computer  for  obtaining 
mathematical  solutions  of  the  stress  distribution  near  Internal  corners 
in  key ways.  Thus  we  formulated  rough  rules  for  deciding  on  the  amount 
of  rounding  necessary  to  prevent  too  big  a  concentration  of  stress  at 
the  internal  corners,  both  in  torsion  and  in  bendirig.  This  work  led  to 
our  thinking  about  what  it  is  that  determines  the  strength  of  materials. 
We  knew  that  consideration  of  intermolecular  forces  leads  to  an 
expectation  that  solids  should  be  much  stronger  than  they  are  and  we 
came  to  the  conclusion,  as  I  think  others  had  already  done,  that  all 
materials  have  cracks  or  faults  in  them  that  entail  far  greater  local 
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stresses  than  the  mean  applied  stress.  This  qualitative  idea  had  little 
value  unless  one  could  form  a  quantitative  model  for  relating  definable 
faults  to  the  reduction  in  mean  stress  which  they  cause.  It  was  a 
brilliant  idea  of  Griffith's  that  made  this  possible.  He  pointed  out 
that  it  is  not  necessary  to  know  the  exact  distribution  of  stress  near 
the  sharp  ends  of  a  crack,  but  only  the  strain  energy  which  is  spread 
over  a  much  larger  area,  and  the  surface  energy  or  surface  tension  of 
the  material.  The  condition  that  the  crack  can  extend  can  be  found  by 
equating  the  loss  of  elastic  energy  during  this  extension  to  the  surface 
energy  gained.  Although  this  idea  proved  very  fruitful  and  was  verified 
by  some  beautiful  experiments  which  Griffith  did  at  Farnborough  with 
brittle  materials,  glass  and  quartz,  it  explained  only  the  reduction  in 
strength  due  to  a  crack.  The  longer  the  cracks  the  weaker  the 
material.  Griffith  and  I  and  Lockspeiser,  who  was  also  interested,  knew 
that  in  fact  most  metals  tend  to  get  stronger  when  they  are  strained  so 
that  the  Griffith  crack  theory  could  not  be  used  to  give  a  theoretical 
picture  of  the  strength  of  metals.  We  made  some  unsuccessful  attempts 
to  modify  the  Griffith  theory  to  include  such  cases.  I  tried  to 
describe  mathematically  the  situation  where  the  stress  distribution 
round  one  Griffith  crack  might  Inhibit  the  growth  of  another,  and  thus 
appear  to  strengthem  the  metal,  but  mathematical  difficulties  defeated 
this  attempt.  It  was  not  until  19  years  later,  in  193*1,  that  I  formed 
the  conception  that  when  the  Griffith  crack  is  a  shear  crack  which  does 
not  pull  the  two  newly  formed  faces  apart  but  allows  them  to  slide  over 
one  another,  the  molecular  force  will  still  operate  across  it.  Under 
these  circumstances  the  faces  may  Join  up  again,  reforming  th*.  crystal 
structure  and  leaving  a  highly  stressed  region  near  both  its  ends.  The 
stress  concentrations  at  the  two  ends  become  what  the  mathematicians 
call  singularities  and  the  stress  distributions  due  to  each  of  these  can 
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be  summed  as  though  they  were  independent,  irrespective  of  how  they  were 
produced.  This  idea  cleared  up  the  mathematical  difficulties,  and  the 
dislocation  theory  of  the  strength  of  metals  started  as  an  offshoot  of 
or,  in  the  modern  phrase,  a  fallout  from,  the  association  of  Griffith 
and  me  in  our  efforts  to  think  about  stress  concentrat i one  in  crankshaft 
keyways  of  aircraft. 

Much  of  our  work  at  Farnborough  consisted  in  designing  instruments 
for  use  in  aircraft  and  it  soon  became  obvious  to  me  that  many  of  our 
difficulties  arose  because  we  could  not  form  a  clear  picture  of  the 
experiences  that  pilots  tried  to  explain  to  us.  In  many  cases  they  did 
not  really  understand  what  they  tried  to  describe  and  often 
unconsciously  described  as  facts  things  that  were  really  only  their  own 
theories  about  the  causes  of  their  experiences.  It  seemed  to  me  that 
I  ought  to  fly  myself  if  I  was  going  to  be  any  use  as  a  deviser  of 
gadgets  for  use  in  the  air.  After  some  discussion,  O'Gorman  came  to 
agree  with  me,  but  the  idea  that  a  scientist  should  fly  appeared 
unorthodox  and  indeed  unnecessary  to  the  Service  chiefs  at  that  time, 
and  they  turned  it  down.  Instead  O'Gorman  released  me  from  his  staff 
and  in  May  1915  I  joined  the  Royal  Flying  Corps  as  a  pupil  in  the 
ordinary  way,  leaving  my  scientific  proclivities  in  abeyance.  I  was 
sent  to  the  flying  school  on  the  old  motor  racing  track  at  Brooklands 
and  took  a  course  on  the  old  1913~'nodel  Maurice  Farman  then  called  a 
Longhorn  or  a  "Rumpety"  for  reasons  which  I  never  fathomed.  It  was  a 
pusher  biplane  with  an  elevator  in  front  and  a  tailplane  supported  on 
long  struts  behind.  Its  maximum  speed  was,  I  think,  about  60  miles  an 
hour  and  it  would  stall  about  37  mph.  One  tried  to  land  it  at  about  42. 
At  that  time  aircraft  were  so  scarce  that  pupils  were  not  allowed  to  fly 
if  there  was  enough  wind  to  put  a  match  out.  We  went  on  to  the  field 
in  the  calm  of  the  early  morning  and  again  in  the  evening.  If  there  was 
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any  wind  we  had  the  day  or  night  off.  It  was  an  easy  life  i)ut  to  keep 
in  it  one  had  to  learn  certain  rules.  In  my  case  tiie  most  important  was 
never  to  touch  the  controls  of  the  dual  control  machine  on  which  I  was 
being  instructed.  The  reason  for  this  was  that  the  instructor  to  whom 
I  had  been  assigned  was  a  nervous  character  who  had  been  been  injured, 
or  at  any  rate  frightened,  in  a  crash  caused  by  a  mistake  of  one  of  his 
pupils,  and  he  decided  this  must  not  happen  again,  so  if  one  of  his 
pupils  applied  any  appreciable  force  to  the  control  or  joy  stick,  as  it 
was  called,  he  reported  that  this  man  was  heavy  handed  and  would  never 
make  a  good  pilot.  At  this  time  many  more  people  wanted  to  learn  to  fly 
than  there  were  machines  to  teach  them  in,  so  the  danger  of  being  thrown 
out  of  the  course  was  a  very'real  one.  The  consequence  was  that  the 
first  time  this  instructor's  pupils  went  solo  was  the  first  time  they 
had  taken  control  and  consequently  they  were  apt  to  do  funny  things  on 
that  occasion.  I  was  fortunate  because  having  been  nearly  a  year  at 
Farnborough  I  knew  what  the  controls  were  expected  to  do.  My  first 
landing  was  too  far  along  the  field  and  ended  in  the  sewage  farm  just 
outside  i t . 

One  of  the  things  that  troubled  the  pupils  was  that  at  the  low 
speeds  our  machines  would  fly,  a  cross-wind  made  one  seem  to  be 
side-slipping  badly  if  one  looked  at  the  ground.  People  were  apt  to 
correct  this  by  manoeuvres  which  gave  rise  to  a  dangerous  real 
side-slip.  Although  I  fear  I  was  a  bad  pilot  when  It  came  to  judging 
distances  I  was  not  bad  In  the  air  and  could  do  the  loops  and  steep 
banking  that  people  did  In  those  days.  I  was  able  to  do  fairly  violent 
manoeuvres  close  to  the  ground  because,  unlike  most  pilots  at  that  time, 
I  had  complete  faith  In  my  instruments  and  did  not  look  at  the  ground 
till  it  was  necessary  to  land  on  it.  A  rumor  even  went  round  that  on 
one  occasion  this  faith  was  misplaced  because  while  I  was  turning,  as  I 
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thought  about  2000  feet  up,  the  tail  of  my  eye  caught  sight  of  a  tree 
quite  close  to.  It  turned  out  that  the  instrument  I  was  relying  on  was 
the  rev  counter,  not  the  altimeter. 

After  getting  my  wings,  as  the  completion  of  the  training  course 
was  called,  in  July  1915  I  did  a  short  course  on  what  was  then  called 
reconnaissance.  This  consisted  in  flying  over  a  farm  which  was  supposed 
to  be  the  target  of  artillery  fire.  People  on  the  ground  let  off  puffs 
of  smoke  in  its  neighborhood  and  we  were  instructed  to  tap  out  in  morse 
code  the  bearing  and  distance  of  the  puff  from  the  target.  At  that  time 
wireless  telephony  had  not  been  invented  and  the  transmitting  apparatus 
had  a  long  wire  aerial  with  a  weight  at  the  end  which  the  pilot  had  to 
let  down  when  he  wanted  to  transmit.  I  failed  at  the  first  attempt 
because  I  forgot  to  wind  it  up  again  before  I  landed  and  it  was  torn  off 
near  the  edge  of  the  airfield. 

When  I  had  passed  the  necessary  tests  O'Gorman  asked  me  to  return 
to  Farnborough  and  rejoin  his  group  of  experimental  officers.  Although 
I  had  not  been  given  a  flying  training  as  a  scientist  I  think  that  the 
fact  that  nothing  very  dreadful  occurred  when  I  was,  may  have 
contributed  to  the  decision  of  the  authorities  about  a  year  later  to 
accept  other  scientists  like  Lindemann,  Farren,  G.  P.  Thomson,  Melvill 
Jones  and  others  as  suitable  candidates  for  this  experience. 

One  of  the  first  things  I  did  on  my  return  to  Farnborough  was  to 
design  and  then  use  apparatus  for  making  pressure  measurements  over  a 
section  of  a  wing  in  flight.  To  understand  the  urge  to  do  this  work  one 
must  think  of  the  state  of  the  art  at  that  time.  Although  Lanchester's 
ideas  about  trailing  vortices  were  known  to  us,  Prandtl's  independent 
development  into  a  usable  form  had  not  penetrated,  and  we  had  not 
thought  of  regarding  wing  drag  as  having  two  separable  '.’omponents,  form 
drag  and  induced  drag.  To  predict  the  aerodynamic  ch-.  racter’  of  the  wing 


USAFA-TR-84-5 


we  could  only  rely  on  tests  of  models  in  a  wind  tunnel  and  there  was 
much  controversy  about  the  magnitude  of  scale  effects.  The  only  method 
for  comparing  model  tests  with  full  scale  performance  was  to  compare  the 
overall  performance  of  the  whole  machine  with  that  predicted  by 
extrapolation  of  the  model  tests  and  in  this  comparison  there  were  so 
many  possible  sources  of  error  that  discrepancies  could  be  attributed 
to  a  wide  variety  of  causes  according  to  one's  fancy.  It  seemed  to  us 
that  direct  measurements  of  the  pressure  distribution  of  a  wing  in 
flight  could  be  made  and  that  a  comparison  could  then  be  made  between 
these  and  similar  measurements  over  a  model  in  a  wind  tunnel.  I  do  not 
think  that  any  measurements  of  that  type  had  been  made  in  1915  when  I 
started,  so  I  had  to  design  the  necessary  apparatus.  I  got  the 
workshops  to  make  a  metal  strip  with  20  pressure  holes  and  insert  it  in 
the  wing  of  a  BE2g  machine.  This  was  a  really  stable  machine.  I  could 
not  have  made  the  measurements  I  did  if  it  had  not  been  stable.  I  made 
a  multitube  manometer  with  photographic  paper  which  could  be  wound  over 
tha  bank  of  tubes,  each  of  which  was  connected  with  one  of  the  pressure 
holes,  and  arranged  a  light  which  I  could  switch  on  for  an  instant  when 
I  had  the  machine  flying  level  and  steadily  on  course  on  a  calm  day. 

Fig.  3  shows  the  results.  In  the  upper  section  the  positions  of  the 
holes  are  shown  and  in  the  lower  the  pressure  coefficient  for  each  hole 
over  a  range  of  speeds  from  50  to  97  mph.  Since  I  only  made 
measurements  over  one  section  I  could  not  achieve  a  comparison  between 
the  known  weight  of  the  aeroplane  and  the  integrated  upward  lift 
produced  by  the  observed  pressure,  but  I  did  get  quite  good  lift  and 
drag  curves  for  the  section,  and  since  the  object  of  the  exercise  was 
to  obtain  a  comparison  between  something  which  could  be  measured 
directly  on  a  full  scale  wing  and  on  a  geometrically  similar  model  wing 
in  a  wind  tunnel,  I  had  done  what  I  set  out  to  do.  When  the  comparable 
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Figure  3-  (Top)  The  position  of  the  pressure  holes 

(Bottom)  The  pressure  coefficients  for  each  hole 
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Before  leaving  this  phase  of  my  aeronautical  x  p  e  r  i  <;■  n  e  I  r  iv’,!. 
mention  one  of  the  less  pleasant  machines  I  used.  It  was  t!a> 
experimental  machine  BE2b  which  had  warped  wings  instead  of  ailcrcris  a  r.  d 
was  longitudinally  very  unstable.  If  it  got  up  to  90  mph  while  diving, 
the  pull  on  the  stick  was  enough  to  require  most  of  my  strength  t(.)  pull 
it  out  of  its  dive.  I  must  confess  that  when  this  happened  I  was  rather 
frightened.  Since  there  must  be  rather  few  people  about  who  have  flown 
warped  wing  machines  I  might-  describe  how  they  operated.  They  had  no 
ailerons.  Wings  in  those  days  were  so  flexible  that  unless  they  were 
stiffened  by  numerous  bracing  wires  they  could  be  twisted  without 
applying  much  force,  until  the  angle  of  incidence  at  the  wing  tip  was 
changed  by  a  degree  or  so.  In  warped  wing  machines  some  of  the  bracing 
wires  were  led  to  the  control  stick  so  that  a  movement  of  the  stick  to 
the  right,  say,  would  decrease  the  angle  of  incidence  of  the  right  wing 
tip  and  increase  that  of  the  left.  They  were  very  uncomfortable  in 
gusty  weather  because  if  a  gust  struck  one  wing  with  greater  force  than 
the  other  the  control  wires  moved  the  control  stick  with  much  greater 
force  than  the  pilot  could  apply.  Thus  although  the  machine  would  bank 
on  application  of  a  steady  force  on  the  stick,  the  stick  itself  would  be 
moving  violently  and  uncontrollably.  I  found  this  most  uncomfortable  in 
the  BE2b  and  also  in  a  Caudron  warped  wing  machine  which  I  used. 

I  left  Farnborough,  I  think,  in  1917,  and  became  a  meteorological 
adviser  to  the  Royal  Flying  Corps.  One  of  the  jobs  I  did  then  was  to 
help  organise  night  flying  experiments  and  in  particular  to  help  set 
up  In  France  a  method  devised  by  Bertram  Hopkinson  for  night  bombing. 
This  was  quite  simple.  Two  searchlights  were  set  up  about  a  mile  apart 
on  the  line  to  the  target.  Their  vertical  beams  could  be  seen  on  most 
nights  owing  to  a  slight  haze  In  the  atmosphere.  The  pilot  got  them  in 


line  and  set  a  compass  course  which  would  take  him  along  the  line 
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Joining  them.  He  started  a  stop-watch  when  passing  over  the  first  and 
stopped  it  at  the  second.  By  multiplying  the  time  thus  recorded  by  the 
ratio  of  the  distance  of  target  to  distance  between  the  searchlights  the 
pilot  estimated  his  time  of  arrival,  assuming  the  wind  remained 
constant.  At  the  estimated  time  he  dropped  his  bombs  and  went  home  — 
at  least  that  was  the  theory,  but  it  was  never  put  into  practice.  We 
had  set  up  the  equipment  in  France  near  the  south  end  of  the  line  and 
were  expecting  to  bomb  Coblenz  the  next  night,  but  happily  peace  was 
declared  during  the  day  and  put  a  stop  to  the  operation. 

Methods  like  those  I  have  described  must  seem  very  amateurish  to 
modern  aviators,  but  it  must  be  remembered  that  the  idea  that  one  could 
navigate  at  night  in  an  aircraft  that  was  heavier  than  air  was  quite  new 
and  very  few  people  could  do  it.  I  think  a  stimulus  which  made  it 
important  to  fly  at  night  was  the  advent  of  night  bombing  raids  by  rigid 
airships.  The  first  time  I  saw  one  of  these,  it  was  coming  down  in 
flames.  I  was  on  leave  one  night  at  my  parents'  home  in  St  Johns  Wood 
when  I  was  wakened  by  anti-aircraft  fire  and  suddenly  I  saw  the  whole 
sky  lit  up.  I  realised  that  something  spectacular  must  be  happening  and 

I  went  to  the  back  of  our  house  and  saw  an  airship  falling  vertically  in 
flames.  I  watched  it  till  it  disappeared  behind  a  house  whose  garden 
backed  onto  ours.  I  noticed  that  it  appeared  just  the  same  length  as  a 
chimney  behind  which  it  disappeared.  I  knew  that  most  airships  were 
then  about  600  feet  long  so  as  soon  as  it  was  light  enough  to  see 
anything  I  got  out  my  sextant  and  measured  the  angle  subtended  by  the 
chimney.  I  also  took  its  compass  bearing  and  so  could  pin-point  the 
position  on  the  map.  This  point  lay  in  a  NNE  direction  and  was  about 

II  miles  away  near  a  place  called  Cuffley.  There  was  a  railway  station 
near  the  spot  on  the  Great  Northern  line.  I  took  the  earliest  train 


there  In  the  morning.  I  was  one  of  the  first  there  and  saw  a  grisly 
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Sight  with  the  bodies  of  the  crew  still  lying  about.  The  whole  ruined 
structure  of  the  airship,  which  turned  out  to  be  a  Schutte^^Lanz ,  was 
confined  to  an  area  about  equal  to  its  cross  section. 

Although  my  only  experience  as  an  aircraft  pilot  was  with  powered 
aeroplanes  and  a  parachute  I  have  been  in  almost  all  other  types  except 
gliders,  as  a  passenger.  Perhaps  the  passage  I  enjoyed  most  was  a  trip 
in  a  free  balloon.  It  was  so  beautifully  quiet.  During  the  1914  War 
balloons  were  sometimes  used  at  night  as  targets  for  training 
searchlight  crews  and  one  night  I  was  one  of  three  passengers  and  a 
pilot  in  the  basket.  We  went  up  from  Kennington  Oval  and  after  a  short 
time  we  went  into  the  clouds  at  about  4000  feet.  Since  we  were  then  of 
no  use  to  the  searchlight  crews  the  pilot  pulled  a  cord  and  released 
some  gas.  We  began  to  go  down  and  very  soon  broke  through  the  clouds. 
Putting  my  hand  out  of  the  basket  I  could  feel  a  strong  upward  breeze 
so  we  were  evidently  descending  pretty  fast.  After  a  time  the  pilot 
thought  it  was  time  to  let  out  some  sand  which  we  carried  as  ballast, 
for  otherwise  we  would  have  hit  the  ground.  The  sand  rapidly  dispersed 
and  fell  more  slowly  than  the  balloon  so,  relative  to  us,  it  went 
upwards.  Then  as  our  downward  velocity  stopped  and  reversed  the  sand 
caught  up  and  filled  the  air  all  round  us.  It  was  like  a  sandstorm. 

The  landing  was  also  interesting.  When  he  wanted  to  land  the  pilot  let 
out  a  long  and  rather  heavy  rope  which  I  should  have  thought  would  have 
broken  any  telephone  wires  it  met.  As  the  rope  dragged  along  the  ground 
it  kept  us  at  a  more  or  less  constant  height  of  about  100  feet  or  less 
and  slowed  up  our  drift.  The  pilot  then  directed  an  electric  torch  beam 
downwards  and  as  the  light  fell  on  trees  or  houses  or  gardens  the 
Illuminated  patch  Jumped  up  and  down.  After  a  time  it  stayed  down  and 
the  pilot  said  we  must  be  over  a  field,  so  he  pulled  the  gas  release 
cord  and  we  descended  with  a  moderate  bump  and  bounced  I  guess  some  30 
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f  eet . 

When  the  war  ended  I  was  expecting  to  return  to  Cambridge  in 
October  1919  as  lecturer  in  mathematics.  I  had  half  a  year  to  spend 
the  gratuity  which  the  Government  handed  out  to  demobilised  officers, 
so  I  was  glad  when  Handley  Page  offered  me  the  Job  of  watching  the 
weather  and  teaching  celestial  navigation  to  the  navigator  of  his  entry 
for  the  Daily  Mail  prize  for  the  first  direct  crossing  of  the  Atlantic. 
The  crew  consisted  of  H.  G.  Brackley,  pilot,  Tryggve  Gran,  a  Norwegian, 
as  navigator,  and  Admiral  Mark  Kerr,  whose  function  I  never  really 
discovered.  The  engineer  who  did  most  of  the  organising  was  Stedraan, 
who  afterwards  became  head  of  the  technical  side  of  the  Canadian  Air 
Force.  There  were,  if  my  memory  is  right,  four  competing  parties  who 
all  assembled  in  St  John's,  Newfoundland,  some  time  in  May  1919. 

Although  four  firms,  Vickers,  Martinsyde,  Sopwith  and  Handley  Page 
entered,  the  affair  was  essentially  a  sporting  event  rather  than  an 
organised  commercial  competition.  There  were  no  aerodromes  then  in 
Newfoundland  and  each  party  set  to  work  to  make  itself  its  own,  erect  a 
hangar  if  needed  (in  fact  I  do  not  think  any  were  completed)  and  then 
build  its  machine.  Fig.  4  is  a  photograph  of  the  Handley  Page  under 
construction  in  the  open.  At  first  we  all  lived  in  the  only  reasonably 
comfortable  hotel  in  St  John's.  Every  morning  each  party  rushed  off  to 
its  own  site  and  tried  to  get  ready  before  the  others.  At  the  beginning 
of  a  month  there  were  cheerful  parties.  This  timing  was  a  side  effect 
of  prohibition.  When  prohibition  was  introduced  in  Newfoundland  doctors 
were  allowed  to  prescribe  alcohol  in  extreme  cases  for  those  patients 
who  needed  it  medicinally.  They  did  this  at  a  dollar  a  time.  One  of 
them,  I  was  told,  prescribed  about  30  000  a  month.  A  rule  was  then 
introduced  which  permitted  each  doctor  to  prescribe  only  a  limited 
number  of  "scrips"  per  month.  At  the  beginning  of  the  month,  therefore, 
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The  first  outfit  to  get  ready  was  the  Martinsyde.  They  crashed  at 
the  take-off  and  were  never  airbourne.  The  navigator,  a  man  called 
Morgan,  was  slightly  injured  I  believe  and  when  the  papers  gave  out  that 
the  machine  was  to  be  repaired  but  that  Morgan  might  not  be  able  to  act 
as  navigator,  people  came  rushing  up  from  all  over  America,  even  from 
the  Argentine,  to  offer  themselves  as  substitutes.  The  hotel  became  so 
crowded  one  could  hardly  move.  The  next  to  start  was  the  Sopwith  with 
Hawker  as  pilot  and  McKenzie  Grieve  as  navigator.  I  saw  them  take  off, 
which  they  did  apparently  without  difficulty.  They  were  not  heard  of 
again  for  over  a  week,  when  they  turned  up  in  a  Norwegian  freighter. 

They  had  been  picked  up  at  sea  and  everyone  marvelled  at  the  ways  of  a 
beneflclent  providence  which  supplied  a  ship  in  the  unfrequented  area 
where  they  came  down  at  exactly  the  time  and  place  necessary  to  save 
them.  The  fact  was  that  when  things  began  to  go  wrong  —  I  think  an 
oilplpe  fractured  —  they  searched  the  sea  for  a  ship  and  ditched 
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themselves  alongside  when  they  found  It. 

The  next  to  take  off  were  Brown  and  Alcock  in  the  Vickers.  That 
was  a  really  spectacular  affair.  Their  aerodrome  was  a  rough  stretch 
on  some  high  ground  behind  St  John's.  They  ran  the  full  length  without 
appearing  to  take  off  and  disappeared  from  our  view  apparently  below  the 
level  of  the  ground  we  were  standing  on.  People  near  where  I  was 
standing  thought  they  must  have  crashed,  but  in  a  few  minutes  they  came 
back  at  a  level  which  seemed  to  be  about  100  feet  and  shot  straight  out 
to  sea:  Brown,  the  navigator,  was  not  very  knowledgeable  about  celestial 
navigation  but  he  pretended  to  be  less  than  he  was,  and  I  remember 
talking  to  him  about  how  he  would  know  where  he  had  got  to  when  he  saw 
land  below  him  on  the  other  side.  He  replied  "Oh  well,  if  I  see  ladies 
in  mantillas  I  will  know  it's  Spain.  If  I  see  people  eating  frogs  I 
will  know  it's  France  and  if  I  see  people  hitting  one  another  over  the 
head  it  will  be  Ireland."  This  is  not  such  a  very  exaggerated  caricature 
of  the  method  of  navigation  actually  employed. 

When  the  Vickers  got  across,  the  Handley  Page,  being  a  lot  bigger 
than  the  others,  had  not  completed  her  flight  tests,  so  they  wisely 
decided  not  to  attempt  the  Atlantic  crossing  but  to  go  down  to  New  York, 
where  they  arrived  after  some  delay  caused  by  a  forced  landing  in  a 
strawberry  field  in  Nova  Scotia.  One  of  the  more  absurd  things  that 
sticks  in  my  memory  was  the  reaction  of  Tryggve  Gran  when  he  heard  the 
news  that  Brown  and  Alcock  had  been  knighted.  He  was  a  very  courageous 
fellow  not  very  good  at  English  and  not  very  clever.  He  had  recently 

married  an  English  actress  and  sent  her  a  telegram  every  day  to  tell  her 

of  his  progress.  When  he  heard  that  Brown  and  Alcock  had  been  knighted 
he  thought  in  his  slow  way,  "If  the  Handley  Page  had  got  across  first  I 

would  have  been  knighted,  and  then  my  wife  would  have  been  Lady  Gran"  so 

his  telegram  to  her  that  day  was,  "Sorry  you  aren't  a  lady".  He 
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couldn't  think  what  was  wrong  with  it. 

In  the  years  1914-19  aeroplanes  had  increased  greatly  in  numbers, 
power  and  reliability  but  the  basic  knowledge  of  fluid  mechanics  which 
was  needed  for  a  proper  understanding  of  how  they  work  had  not  advanced 
so  rapidly.  The  advantage  of  reducing  resistance  by  streamlining  struts 
and  wires  was  appreciated,  but  it  was  not  until  Melvill  Jones  pointed 
out  that  the  part  of  the  drag  due  to  skin  friction  and  turbulence  could 
theoretically  be  reduced  to  near  that  of  a  flat  plate  of  the  same  area 
as  the  wing  that  designers  had  an  ideal  against  which  they  could  measure 
their  performance.  The  conception  that  drag  is  separable  into  two  parts 
is  a  natural  consequence  of  the  Lanchester-Prandtl  theory  of  lift.  The 
combination  of  this  theory  with  the  Jou  "■wsky  theory  of  how  circulation 
is  related  to  aerofoil  section  and  the  boundary  layer  theory  of  surface 
friction,  provided  the  main  themes  of  aerodynamic  research  until  speed 
increased  so  much  that  effects  of  compressibility  became  Important. 

Although  these  ideas  provided  a  basis  for  design  there  was  a  gap 
between  them  and  the  classical  theory  of  hydrodynamics,  so  well 
expounded  in  the  great  treatise  of  Horace  Lamb.  Although  some  of  the 
conclusions  of  classical  hydrodynamics  had  been  verified  experimentally, 
those  on  waves  and  on  slow  viscous  flow  are  instances,  those  on  flow 
past  solid  bodies  had  not.  In  particular  the  sphere,  geometrically  the 
simplest,  is  aerodynamically  one  of  the  most  complex  of  bodies,  and  the 
inability  of  theory  to  describe  flow  round  what  seemed  to  be  the 
simplest  possible  shape  told  heavily  against  it  in  the  minds  of 
practical  engineers.  When  I  returned  to  Cambridge  in  1919  I  aimed  to 
bridge  the  gap  between  Lamb  and  Prandtl.  In  following  that  course  I  was 
fortunate  enough  to  interest  Rutherford  who  gave  me  a  room  near  his  own 
in  the  Cavendish  Laboratory  to  work  in.  The  reason  for  the  discrepancy 


between  classical  theory  of  bodies  moving  in  a  fluid  and  what  is 
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Observed  was  understood  at  that  time.  The  classical  theory  permitted 
perfect  slipping  at  solid  surfaces  and  the  existence  of  viscosity  makes 
this  surface  condition  unrealistic.  I  therefore  looked  for  effects  in 
the  motion  of  fluids  which  do  not  depend  on  surface  conditions.  One  is 
that  of  rotation  which  affects  fluid  flow  through  the  body  of  the  fluid 
rather  than  at  the  solid  surface.  Here  I  found  that  classical  theory 
could  predict  exactly  in  some  cases  how  motions  which  themselves  would 

not  be  described  by  that  theory  would  be  modified  by  rotation. 

Encouraged  by  this  success  I  looked  for  and  found  a  number  of  other 
things  that  could  be  predicted  correctly  by  classical  theory.  Since 
then,  I  have  witnessed  an  increasing  willingness  on  the  part  of 
designers  to  use  the  work  of  theoretical  hy dr odyn am i c i s t s . 

At  this  stage  I  think  it  is  time  to  stop  before  I  expose  ray 
ignorance  of  modern  aeroplanes.  I  am  glad  to  have  lived  during  the  time 
they  were  simple  enough  for  a  reasonably  intelligent  amateur  to 
appreciate  them  as  mechanisms.  They  have  developed  in  the  kind  of  way, 
but  much  faster  than  we  thought  possible  *10  years  ago.  They  have  also 
had  some  of  the  regrettable  effects  on  our  lives  that  I  foresaw  and 
expressed  in  my  Wilbur  Wright  lecture  to  this  Society  years  ago.  I 
then  pointed  out  that  every  improvement  in  transport  facilities  makes  it 
less  worth  while  to  travel  on  a  globe  of  limited  size  because  the  place 
of  arrival  becomes  more  and  more  like  the  place  of  departure.  Although 
aeroplanes  have  not  had  so  detestable  an  effect  on  our  country  as 
automobiles,  the  activities  of  travel  agents  are  forcing  them  along  the 
same  destructive  path.  You  may  think  this  an  old  man's  view,  but  it  is 
what  I  thought  when  I  was  young  and  what  I  still  think. 

By  way  of  contrast  with  the  primitive  state  of  aeronautics  sketched 
in  this  lecture  I  may  mention  that  since  writing  it  I  travelled  in 
comfort  from  Boston  to  London  in  a  VC10  —  neither  ground  nor  sea  were 
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visible  at  any  time  during  our  flight.  Shortly  after  leaving  the 
American  coast  our  Captain  told  us  that  were  were  under  the  guidance  of 
an  automatic  pilot  controlled  from  London. 


1  10 
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UNCERTAINTY  ANALYSIS:  A  NEW  COMPUTER  SUBROUTINE 
Christopher  A.  Boedicker* 

Editor's  Note 

An  earlier  version  of  this  paper  won  third  place  in  the  under¬ 
graduate  division  of  the  1983  American  Institute  of  Aeronautics  and 
Astronautics  Region  V  student  paper  competition. 

Abstract 

At  the  USAF  Academy,  I  have  developed  and  implemented  a  computer 
subroutine  (UNCERT)  that  greatly  simplifies  the  immense  problems 
associated  with  uncertainty  analyses.  This  paper  describes  the  theory 
behind  the  uncertainty  analysis  and  expl  ins  how  this  theory  can  be 
implemented  in  a  computer  program. 

I .  Introduction 

Uncertainty  analyses  are  vital  for  the  test  engineer.  Consider 
this  example.  A  firm  wishes  to  purchase  a  turbine  engine  compressor, 
and  the  only  criterion  of  Importance  is  polytropic  efficiency.  Company 
A  boasts  an  efficiency  of  90  percent,  and  their  number  is  guaranteed  to 
within  *  0.1  percent.  Company  B  has  a  compressor  with  a  91  percent 
efficiency,  but  this  number  can  only  be  guaranteed  to  within  +  2 
percent.  Which  compressor  should  the  firm  buy?  The  importance  of 
uncertainty  analysis  may  be  somewhat  clearer  now. 

Analyzing  the  uncertainty  before  a  test  gives  the  engineer  an 
indication  of  the  expected  accuracy  of  the  test  results.  In  the  absence 
of  universally  accepted  standards,  it  is  usually  the  engineer  who 
determines  the  uncertainty  values  that  are  acceptable  for  a  particular 
test.  For  example,  in  a  test  for  which  only  mass  flow  and  Mach  number 
are  calculated,  an  engineer  may  specify  that  he  needs  to  know  mass  flow 
within  t  0.5  percent,  but  a  Mach  number  within  t  3  percent  is 
acceptable.  If  pretest  uncertainties  are  too  large,  corrective  action 
can  often  be  taken.  Specific  methods  for  reducing  uncertainty  will  be 
discussed  later. 

I  have  developed  a  computer  subroutine  called  UNCERT  to  facilitate 


*Cadet,  USAF  Academy 


1 S AFA-  I  R- ^4  -  4 


uncertainty  analysis.  One  can  attach  UNCERT  to  a  computer  program  tfiat 
contains  performance  calculations  in  a  separate  subroutine.  Pretest 
analyses  show  the  test  engineer  where  problems  exist  so  corrective 
action  can  be  taken,  if  necessary.  Once  the  test  has  been  completed, 
UNCERT  provides  the  user  with  a  reportable  measure  of  the  accuracy  of 
the  results.  This  paper  discusses  the  theory  behind  the  uncertainty 
analysis,  explains  how  to  use  the  subroutine,  and  presents  the  results 
of  an  actual  application  of  the  subroutine. 

1 1 ,  Theory 

Uncertainty  analysis  concerns  the  mathematical  reliability  of  an 
experiment.  A  detailed  discussion  of  the  theory  behind  uncertainty 
analysis  appears  in  Ref.  3:  the  paragraphs  below  summarize  the  main 
ideas  from  that  article. 

In  looking  at  the  results  of  an  experiment,  one  must  distinguish 
between  error  and  uncertainty.  Error  is  a  certain  fixed  value  and 
refers  to  the  difference  between  a  measurement  and  the  true  value  of  a 
quantity.  Uncertainty,  on  the  other  hand,  is  a  measure  of  accuracy  and 
represents  a  range  of  values  within  which  the  result  is  expected  to 
fall . 

Two  types  of  error  are  associated  with  any  measurement:  random 
(precision)  error  and  fixed  (bias)  error.  Random  errors  occur  because 
small,  Independent  Influences  prevent  the  measurement  system  from 
delivering  the  same  reading  when  supplied  with  the  same  input  value. 
Deviation  from  the  mean  approximates  a  normal  distribution  (Figure  1). 
Standard  deviation  is  used  as  d  measure  r,r  tiie  precision  error,  and  the 
statistic  S,  called  the  ;'-e,:'<i  index,  is  used  to  estimate  the 


standard  deviation. 


.  s  ^  i  V  ••  r 
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where  N  is  the  number  of  measurements,  and  X  is  the  average  value  of 
individual  measurements  X^.  By  making  Z  independent  measurements,  the 
standard  deviation  can  be  reduced  by  Jz ,  as  shewn  in  Eqn.  (2). 

S 

s  _ _ 


Bias  errors  will  cause  repeated  readings  to  be  incorrect  by  roughly 
the  same  amount  but  for  unknown  reasons.  Bias  is  determined  by 
comparing  a  measurement  with  a  known  standard,  an  often  impossible  task 
in  a  laboratory  situation.  Judgment  of  the  most  knowledgeable 
Instrumentation  expert  is  often  relied  upon  to  obtain  bias  limits. 

Figure  2  illustrates  how  the  bias  error  combines  with  the  precision 
error  to  form  measurement  error. 

The  bias  and  precision  error  estimates  are  combined  to  yield  an 
uncertainty  interval.  The  National  Bureau  of  Standards  (NBS)  and  much 
of  industry  recognize  a  standard  represented  by  the  bias  limit  plus 
a  multiple  of  the  precision  index.  The  uncertainty  is  centered  about 
the  measurement  value  and  is  evaluated  using 

U  =  iCB  tg^S)  (3) 

where  B  is  the  bias  limit,  S  is  the  precision  index,  and  t 95  is  the  95th 
percentile  for  the  2-tailed  Student's  "t"  distribution.  tg^  is  based  on 
the  sample  size.  When  t  =  2  (the  value  usually  selected),  Eqn.  (3) 
becomes 

U  =  1  (  B  +  2S)  ( 14  ) 

This  uncertainty  value  does  not  specify  a  statistical  confidence 
interval,  since  the  bias  limit  is  based  on  Judgment.  Uncertainty,  as 
stated  earlier,  is  the  largest  error  expected  in  the  parameter  of 
interest.  The  uncertainty  Interval  is  shown  in  Figure  3. 
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Measurement 


Figure  3.  Measurement  Uncertainty  (Ref.  1) 


Performance  parameters  are  rarely  measured  directly  in  a  test. 
Basic  quantities  like  temperature,  pressure,  force,  and  position  are 
determined,  and  performance  parameters  are  then  evaluated  as  a  function 
of  these  more  basic  measurements.  Errors  in  the  measurements  are 
propagated  to  the  performance  parameters  through  the  particular 
equation  for  the  parameter,  and  this  propagation  may  be  approximated  by 
a  first-order  Taylor  series  method. 

Consider  a  test  in  which  there  are  L  performance  parameters  or 
calculations  based  on  n  measurements.  The  bias  limit  for  the  Jth 
calculation  may  be  determined  using 


where  c,  is  the  jth  calculation,  is  a  measurement,  and  B is  the 
bias  limit  of  the  kth  measurement.  The  j.arge  M  subscript  specifies  an 
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estimate  of  the  error  for  a  measurement:  the  large  C  subscript  specific 


an  estimate  of  the  error  for  a  calculation.  A  similar  method  is  used 


to  determine  the  precision  index  for  each  of  the  calculations. 

rz — - ~2 


n  .3C. 


k  =  iV\ 


where  Sv,  represents  the  precision  index  associated  with  the  kth 


measurement.  Once  these  numbers  have  been  obtained,  they  may  be  added 


using 


•Jj  =  i  > 


to  Obtain  the  uncertainty  in  the  jth  calculat ion , Ined  parameter  is  the 


III.  Reasons  for  Computer  Appl Ication 


Although  methods  to  determine  uncertainties  are  available,  some 


tests  are  so  complex  that  analytic  solutions  to  the  problem 


cannot  be  obtained  easily,  if  at  all.  The  propagation  of  errors 


(through  numerous  measurements)  may  also  prevent  an  analytic  solution. 


Numerical  solutions  must  therefore  be  used.  The  partial  derivatives  in 


Eqns.  (5)  and  (6)  must  be  evaluated  using 


aC.  AC. 
J  J 


3Mk  AM^ 


where  C  is  a  calculated  value,  and  M  is  a  measurement.  To  facilitate 


calculations,  a  non-dimensional  number  should  be  used: 


r  ^  * 


^1^  \  Cj/  nominal 


This  is  the  sensitivity  coefficient,  expressed  in  percent  per  percent. 


At  this  point,  we  can  use  a  time-saving  computer  subroutine  to 


calculate  these  sensitivity  coefficients  and  then  combine  them  with  all 
the  other  proper  computations  to  determine  uncertainties. 


IV .  General  Working  of  the  Subroutine 

The  performance  of  a  test  article  is  usually  determined  by  a 
computer  program.  The  main  part  of  the  program  is  used  for  inputs,  and 
the  performance  parameters  are  calculated  in  a  separate  subroutine.  If 
such  a  program  already  exists,  the  UNCERT  subroutine  can  be  attached  in 
order  to  yield  the  uncertainties  in  the  calculations. 

Bias  limits  (Bw  's)  and  precision  Indices  (Sm  's)  for  each  of  the 

k 

measurements  must  be  input  into  DATA  statements  within  the  UNCERT 
subroutine.  The  values  must  be , determined  from  manufacturer's  data  and 
expert  opinions. 

To  determine  the  uncertainties,  we  must  numerically  evaluate 
partial  derivatives.  First,  a  1  percent  change  in  the  measured 
parameters  is  made.  Then  the  subroutine  that  calculates  the  performance 
parameters  (TAPERF)  is  called  (once  for  each  measured  value)  to  obtain 
a  new  value  for  the  performance  parameter  to  correspond  to  the  new  value 
of  the  measured  parameter.  The  partial  derivatives  are  then  computed  as 


(C  )  -  C  .)  /  Cj 

(M  ■  - 


(10) 


Where 

M',  =  1  .01  M,  (11) 

K  k 

and  C'j  is  the  calculated  value  (performance  parame'ter)  evaluated  using 
M ,  Once  the  measurement  bias  limits,  precision  indices,  and  none 
dimensional  partial  derivatives  are  available,  the  bias  limits  and 


precision  Indices  for  the  performance  parameters  are  found  using 
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and 


These  can  be  evaluated  through  matrix  multiplication. 


(12) 


(13) 


ft  X  n]  fn  X  l]  =  [l  X  l] 

where  L  and  n  are  the  maximum  number  of  calculations  and  measurements, 

2 

respectively.  The  values  are  calculated  in  the  same  manner.  A 

j 

2  2 

general  matrix  product  subroutine  evaluates  Bq  and  S values,  square 

j  j 

roots  are  taken,  and  then  uncertainties  in  the  performance  are 
determined  using 

Uj  =  t(Bc^  *  t95  Sc,)  (7) 

where  tgj  is  generally  set  equal  to  2.0.  The  subroutine  outputs  the 
final  values  of  and  Sq  ,  the  equation  used  for  the  uncertainty 

j  j 

calculations,  and  . 


V .  Application  of  Subroutine 

UNCERT  has  been  created  in  order  to  perform  uncertainty  analyses. 
A  listing  of  the  subroutine,  as  used  within  a  specific  example,  and  an 


« 

•  i 

3 


■j 


'i 
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explanation  of  how  the  subroutine  can  be  used  appear  in  Appendix  A. 
Because  of  the  length  of  the  program,  only  those  portions  necessary  to 
understand  UNCERT  will  be  examined. 


)lication  Example 


The  UNCERT  subroutine  has  been  used  to  calculate  uncertainties 


for  a  centrifugal  compressor  test  at  the  Wright-Patterson  AFB  Aero 
Propulsion  Lab.  The  main  program,  used  to  input  the  proper  measure¬ 
ments,  was  called  RUNTA  (Run  Test  Article),  and  the  performance  was 
computed  in  a  subroutine  called  TAPERF  (Test  Article  Performance). 

The  relevant  parts  of  this  program  are  listed  in  Appendix  A. 

The  compressor  to  be  evaluated  was  designed  and  manufactured  by  the 
Garrett  Turbine  Engine  Company.  Figure  i)  shows  a  cross  section  of  this 
compressor,  with  measurement  station  designations.  This  compressor  has 
variable  geometry  with  a  flow  rate  of  about  25  Ibm/sec  of  air  at  design 
speed  (23i000  rpm),  and  it  develops  an  8:1  pressure  ratio  in  the  single 
stage.  RUNTA,  using  TAPERF,  computes  64  different  performance/aero¬ 
dynamic  parameters  based  on  20  different  measurements.  The  continuity 
equation,  velocity  diagrams,  and  the  Euler  Turboraachlnery  equation  are 
used  with  1  ^--dimensional ,  steady-flow  assumptions  to  evaluate  the 
compressor's  performance.  We  also  employ  iterative  schemes  in  some 
cases  to  find  certain  performance  parameters.  The  final  calculations 
Involve  the  use  of  simple  definitions,  such  as  pressure  recoveries  or 
ef f 1 c lenc ies . 

The  results  of  the  RUNTA  program  appear  in  Appendix  B.  The 
nominal  performance  values  are  listed  first.  Then,  in  accordance  with 
Dr.  Abernethy’s  recommendations  (Refs.  1  and  2),  the  equation  for  the 
uncertainty,  the  bias  limits,  the  precision  indices,  and  the  final 
uncertainty  values  (in  +  percent)  are  given  as  output. 
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2.9  3.0 


Figure  *) .  Compressor  Flowpath  and  Station  Designation 

The  results  provide  some  interesting  insights  into  the  uncertainty 
analysis.  Some  of  the  uncertainties  found  are  in  excess  of  +10  percent. 
An  examination  of  the  data  shows  that  several  different  factors  may  have 
led  to  these  high  uncertainty  values.  First,  the  method  we  used  to 
calculate  the  performance  parameter  may  have  resulted  in  high  values  for 
the  Influence  coefficients.  At  the  test  article's  operating  point,  the 
performance/aerodynamic  parameter  may  have  been  highly  sensitive  to 
small  changes  in  a  measured  value.  Second,  the  instruments  used  to 
determine  the  value  of  the  measurement  may  have  had  high  bias  limits 
and/or  high  precision  indices.  Finally,  the  uncertainty  may  have 
resulted  from  the  point  of  operation  of  the  test  article.  For  example, 
transducers  are  most  accurate  when  they  are  operating  near  100  percent 
of  their  range.  At  lower  values,  the  uncertainty  associated  with  a 
transducer  can  Increase  significantly. 
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In  looking  at  the  results,  one  immediately  questions  the 
uncertainty  of  +0  percent  listed  for  C(2).  C(2)  represents  the  pressure 
loss  coefficient  associated  with  the  inlet  guide  vane  (IGV)  of  the 
compressor  (shown  in  Figure  ^).  This  coefficient  is  evaluated  using  a 
look-up  function  in  the  performance  subroutine,  and  it  depends  on  the 
angle  given  to  the  inlet  guide  vane.  The  loss  coefficient  does  not 
change  its  value  until  the  IGV  angle  approaches  30  degrees.  For  this 
example,  the  IGV  angle  is  .01  degrees.  This  small  nominal  value  is  used 
since  a  value  of  0  degrees  causes  DELTAC,  (C'-C)/C,  to  be  undefined. 
Small  changes  in  the  IGV  angle,  for  its  small  nominal  value,  do  not 
change  the  value  of  the  loss  coefficient.  The  influence  coefficients 
for  this  example  are  all  0;  thus,  no  uncertainty  is  found  in  the  loss 
coefficient. 

An  in-depth  analysis  of  one  of  the  uncertainty  values  may  shed  more 
light  on  the  subject.  The  38th  calculation  in  the  program  involves  the 
absolute  Mach  number  at  station  2.9  (Figure  4),  the  diffuser  exit.  The 
uncertainty  in  this  calculated  parameter  [C(38)]  is  +5.4  percent,  as 
listed  in  Appendix  B.  Table  I  below  shows  the  bias  limits,  precision 
indices,  and  uncertainties  for  the  5  absolute  Mach  numbers  calculated 
in  this  example. 


Table  I 

BIAS  LIMITS  (B^  ),  PRECISION  INDICES  (Sc  ),  AND  UNCERTAINTIES  (U,) 

j  j 

FOR  5  MACH  NUMBER  CALCULATIONS 


Parameter 

J 

Uj  ( 

t%) 

AMI  9 

1 

.6867  E-02 

.3265  E-02 

.  1  3 

E  +  01 

5 

.7249  E-02 

.3447  E-02 

.  1  4 

E  +  01 

20 

.5162  E-02 

.5064  E-03 

.62 

E  +  00 

AM29 

38 

.1935  E-01 

.1732  E-01 

.54 

E  +  01 

AM30 

42 

.2126  E-01 

.1458  E-OI 

.50 

E  +  01 
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The  uncertainties  in  the  absolute  Mach  numbers  at  stations  2.9  and  3.0 
are  both  over  5  percent.  Looking  at  Table  I,  one  can  see  that  the  bias 
limits  and  precision  indices  for  these  calculations  are  significantly 
higher  than  those  for  the  other  3  Mach  numbers.  Looking  back  at  Eqns. 
(12)  and  (13)i  one  can  see  that  high  values  for  the  partial  derivatives 
or  the  measurement  bias  limits  or  precision  indices  would  cause  high 
values  in  8^  and  S  ,  and  thus  large  uncertainties.  If  a  problem 

j  j 

exists,  there  are  several  ways  to  correct  it.  Lower  values  for  B., 
should  not  be  input  in  order  to  lower  the  overall  uncertainty.  The 
first  estimate  was  based  on  expert  opinion,  and  changing  this  value 
until  a  desired  result  is  achieved  defeats  the  purpose  of  the 
uncertainty  analysis.  The  S|.(^  values,  however,  could  be  lowered  in 
order  to  alleviate  high  uncertainties.  Either  more  or  better 
transducers  could  be  used  to  lower  the  S;.]^  values.  More  could  be  used 
since 


S 


(1 


where  Z  is  the  number  of  transducers  whose  outputs  are  averaged  together 
to  give  the  measured  value.  However,  cost  and  space  available  for 
transducers  may  be  prohibitive.  A  final  method  to  lower  the  uncertainty 
would  be  to  change  the  equation  used  to  calculate  the  parameter  of 
Interest.  A  different  equation,  if  one  exists,  might  favorably  alter 
the  value  of  the  partial  derivative. 

We  can  use  the  Mach  number  example  to  Illustrate  some  of  the 
problems  that  might  be  encountered  using  a  given  calculation  method. 
Since  the  bias  limit  and  precision  index  for  any  given  calculation  are 
based  on  partial  derivatives,  the  equation  becomes  very  important.  If 
static  and  total  pressure  are  known  at  a  given  point,  Mach  number  may  be 
calculated  using 


123 
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The  influence  coefficients  would  be  obtained  by  multiplying  the  above  2 

I  I 

equations  by  —I  and  —  I  ,  respectively.  Note  that 
M  I  nom  M  I  nom 

3M  Pg  3M 

M  3P^  M  3P^  ^ 

The  change  in  the  influence  coefficients  in  the  above  equation  with 
Pg/P^  (or  Mach  number)  is  shown  in  Figure  5. 

Figure  5  shows  that  a  Mach  number  calculation  using  Eqn.  (16)  is 
quite  sensitive  to  errors  in  P^  or  Pg ,  especially  for  low,  subsonic  Mach 
numbers.  Table  II  is  a  list  of  some  of  the  partial  derivatives  of  the 


Mach  numbers. 


AMI  9 

PS25 

0 

PS29 

0 

PS30 

0 

PT29 

0 

PT30 

0 

AM20 

0 

0 

0 

0 

0 

(C.)  AM25 

AM29 

.1261)  E  +  00 

0 

0 

0 

0 

.  1)306  E-03 

.8898  E+01 

0 

.8191  E+01 

0 

AM30 

.1)311)  E-03 

0 

.9815  E-01 

0 

.8958  E+01 

The  partial  derivative  values  shown  above  actually  represent  the 
percentage  change  in  a  calculated  value  that  can  be  expected  for  a  1 
percent  change  in  a  given  measured  value.  Thus,  one  can  expect  an  8.898 
percent  increase  in  the  absolute  Mach  number  at  station  2.9  for  a  1 
percent  Increase  in  the  static  pressure  at  that  station.  A  negative 
value  for  a  sensitivity  coefficient  means  that  one  can  expect  a  decrease 
in  the  calculated  value  when  the  measurement  is  Increased.  The  partial 
derivatives  for  the  Mach  numbers  at  stations  2.9  and  3.0  (AM29  and  AM30) 


I  S-'  l  A  -  rR-H4-  => 


are  very  large  when  taken  with  respect  to  static  and  total  pressures 
near  the  compressor  exit.  At  stations  1.9  and  2.0,  the  air  flow  angle 
is  known,  and  the  continuity  equation  is  applied  in  a  rather  simple 
iterative  scheme  to  find  the  Mach  numbers.  At  the  impeller  exit 
(station  2.5),  the  air  flow  angle  is  not  known.  Thus,  we  use  the  Euler 
Turbomachinery  equation  with  the  continuity  equation  in  a  double 
iteration  to  determine  Mach  number.  At  stations  2.9  and  3.0,  we  use 
static  to  total  pressure  ratios  to  find  Mach  number.  But  the  possible 
problems  with  this  last  method  have  just  been  illustrated.  The  test 
engineer  may  wish  to  change  the  calculation  procedure  if  the  situation 
allows.  An  iterative  scheme  could  be  employed  to  find  the  Mach  numbers 
at  stations  2.9  and  3.0.  The  desired  uncertainty  levels  and  the 
feasibility  of  the  change  will  determine  whether  or  not  the  calculation 
procedure  should  be  changed. 


UNCERTAINTY 
%  READING 


Figure  6.  Uncertainty  Variation  with  Instrument  Operating  Range 

The  operating  point  of  the  test  article  is  very  important.  The 
uncertainty  results  apply  to  only  one  operating  point.  It  has  already 
been  shown  that  the  operating  point  affects  the  influence  coefficient; 


it  may 

also  affect 

the  bias 

;  1  m  i  t 

a  r.  p  r*  iM  "  i  v'  ri 

i  n  'j  .< 

i :)  .‘j  o 

given 

transducer  . 

Figure  6 

should 

clarify  thl.3 

point  . 

I  f 

article  Is  not  operating  so  that  the  transducers  and  instrumentation  are 
near  100  percent  of  their  prescribed  range,  then  u n oer t a i n t  i  es  can  be 
higher  than  those  obtained  at  a  "design  point."  If  the  test  must  be 
performed  far  from  the  design  point,  then  better  or  different 
transducers  might  be  needed,  meaning  large  expenditures  of  time  and 
money.  Unfortunately,  in  many  cases  nothing  can  be  done.  Either  the 
uncertainties  must  be  accepted,  or  the  test  must  be  cancelled.  If  the 
uncertainty  analysis  is  performed,  however,  it  still  enables  the  test 
engineer  to  tell  how  accurate  the  answers  from  that  test  will  be. 

In  summary,  high  uncertainties  can  result  from  many  different 
sources.  Calculation  methods  may  lead  to  large  sensitivity 
coefficients.  If  this  is  true,  new  calculation  procedures  should  be 
sought.  Also,  the  bias  limits  and  precision  indices  may  be  too  large. 
Better,  more  expensive  transducers  can  be  used  to  alleviate  this 
problem.  The  measurement's  precision  index  may  be  reduced  by  measuring 
the  same  variable  with  more  transducers.  Finally,  the  point  of 
operation  of  the  test  article  may  cause  high  uncertainties.  Better 
transducers  could  be  used  in  this  situation,  or  the  engineer  might 
choose  to  change  transducers,  depending  on  the  specific  point  of 
operation.  The  feasibility  of  adopting  many  of  these  options  depends 
on  the  time  and  money  available  to  the  test  director. 

VII.  Summary 

The  UNCERT  subroutine  has  been  developed  to  perform  uncertainty 
calculations  for  various  experimental  test  articles  with  existing 


performance  programs.  The  work  required  to  attach  the  subroutine  and 
make  the  required  additional  statements  is  not  excessive,  especially  if 
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one  considers  that  at  present  most  uncertainty  analyses  are  either  very 
time-consuming  or  are  never  accomplished.  Using  the  UNCERT  subroutine 
before  a  test  is  conducted  can  identify  potential  problem  areas  so  that 
the  test  engineer  can  decide  whether  to  continue  the  test.  If  the 
alterations  lead  to  excessive  costs,  the  engineer  might  cancel  the  test 
before  wasting  large  amounts  of  money.  If  the  test  is  conducted, 
results  of  the  final  uncertainty  analysis  can  be  compared  to  the  pretest 
analyses.  The  UNCERT  subroutine  can  provide  a  valuable  resource  for  any 
test  facility. 
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Symbols 

bias  limit  associated  with  a  performance 
parameter 

bias  limit  of  an  individual  measurement 
elemental  bias  limit 
the  Jth  calculation 

total  number  of  calculated  parameters 

the  kth  measurement 

number  of  samples 

total  number  of  measurements 

precision  index  associated  with  a  performance 
parameter 

precision  index  of  an  individual  measurement 
elemental  precision  index 
estimate  of  the  standard  deviation,  or 
precision  index 
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Greek  Symbols 

A  A 


(T 


Superscripts 


Subscr 1 Dts 

J 
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ANSWER 


BCJ 

BMK 

C 

CPRM 


DC 


uncertainty;  an  estimate  of  the  error 
an  observed  value:  an  individual  measurement 
sample  average  of  measurements 
number  of  transducers 

sensitivity  (influence)  coefficient 
degrees  of  freedom  involved  in  a  sample 
true  standard  deviation  of  repeated  values  of 
measurement 

represents  a  new  measurement  or  calculation 
value  when  a  specified  change  has  been  made 
in  a  measured  value 
denotes  a  non-dimensional  value 

specifies  a  number  from  1  to  L 
specifies  a  number  from  1  to  n 
Computer  Variable  Names 
final  uncertainty  value  determined  by  U  = 

B  +  2S 
B(-  values 

j 

B  ^  values 
calculated  value 

(C')  represents  the  value  of  a  calculated 
parameter  after  it  has  been  evaluated  using  a 
value  of  M' 

dummy  array  containing  the  performance 
parameters  calculated  by  the  performance 
equations 
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Appendix  A 
Computer  Subroutine 


To  use  the  UNCERT  subroutine,  one  requirement  must  be  met:  the 
performance  of  the  test  article  must  be  evaluated  in  a  separate 
subroutine.  I  will  now  illustrate  what  must  be  done  to  implement  the 
UNCERT  subroutine. 

The  computer  program,  which  has  been  named  RUNTA  (Run  Test 
Article),  is  listed  on  lines  100  through  14000.  The  UNCERT  subroutine 
is  shown  on  lines  14100  through  28400,  and  the  general  matrix  product 
subroutine  (GMPRD)  is  listed  on  lines  28600  through  303OO,  The 
performance  calculations,  contained  in  TAPERF  (Test  Article 
Performance),  are  given  in  lines  30400  through  87400. 

The  main  part  of  the  program,  along  with  the  performance 
subroutine,  must  be  modified  in  order  to  evaluate  the  required  partial 
derivatives.  A  1  percent  change  is  made  in  the  first  measurement, 
resulting  in  M'j*  The  performance  subroutine  is  then  called,  and  C'@ 


values  are  determined  based  on  The  partial  derivative  is 


(c!  -  c.)  /c. 

'  -1  j  ) 


K  -  V 


/M, 


The  nominal  values  of  the  performance  parameters  and  the  measurements 


(Cj  and  M  )  are  stored  and  used  to  calculate  the  influence  coefficients 


(ICj/lM^). 


1 


is  then  converted  back  to  its  nominal  value,  M 


The 


second  measurement  is  changed,  performance  parameters  are  evaluated,  and 


the  influence  coefficients  associated  with  are  determined.  This 


process  continues  until  all  of  the  influence  coefficients  have  been 
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i 

computed . 

For  the  UNCERT  subroutine  to  evaluate  sensitivity  coefficients 
correctly,  certain  assignment  statements  must  be  added  to  the  main 
program  and  the  performance  subroutine.  In  the  main  program,  we  must 
place  the  various  measurements  into  the  M  array.  For  example,  if  there 
are  five  different  measurements,  then  the  required  assignment  statements 
would  resemble 

100  M(1)  =  Measurement  Variable  Name  #1 

200  M(2)  =  Measurement  Variable  Name  #2 

300  M(3)  =  Measurement  Variable  Name  #3 

1100  M(i4)  =  Measurement  Variable  Name  #4 

500  M(5)  =  Measurement  Variable  Name  #5 

The  variable  names  would  correspond  to  the  measured  parameters  on  which 

the  performance  calculations  depend.  This  can  be  seen  on  lines  8800 

through  10700  in  this  Appendix. 

In  the  performance  subroutine  (TAPERF),  we  place  all  of  the 

performance  calculations  into  a  dummy  array  at  the  end  of  the 

subroutine.  This  is  seen  on  lines  79800  through  86100.  The  values 

stored  in  the  dummy  array  are  then  placed  into  either  a  Cj  or  C'j  array. 

When  we  use  the  nominal  values  of  the  measurements  to  evaluate  the 

performance  parameters,  we  can  then  place  and  store  the  values  of  these 

calculations  in  the  Cj  array.  When  a  1  percent  change  in  a  measurement 

(producing  M'j^)  is  made,  the  resulting  performance  calculations  are 

placed  in  the  C'j  array.  Once  the  C'j  values  have  been  stored,  we  can 

calculate  the  Influence  coefficients  using  Eqn.  (A1).  The  user  must 

employ  a  "flag"  so  that  the  program  can  determine  whether  to  place  the 

performance  calculations  in  the  Cj  and  C'j  array.  The  flag,  IFLAG,  must 

be  placed  into  COMMON  blocks  in  both  the  main  program  and  the 

performance  subroutine  (lines  2000  and  37800,  respectively).  Before 

calling  the  performance  subroutine  in  the  main  program,  the  flag  is  set 

to  0  (line  10900).  In  the  performance  subroutine,  the  code  appears  on 
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lines  861400  through  87100.  The  DC  array  is  a  dummy  array.  For  nominal 

performance  parameters,  we  set  the  flag  at  0  and  transfer  the  D  (j  array 

J 

values  into  the  C  .  array.  Then  the  flag  becomes  1  (line  86500)  after  we 

have  made  the  initial  assignments.  The  C .  values  are  stored  and  can 

J  J 

be  recalled  when  required.  When  the  flag  is  1  (whenever  one  of  the 
measured  values  has  been  altered  by  1  percent),  we  place  the  dummy 
values  into  the  C'j  array  so  we  can  determine  the  influence 
coefficients. 

UNCERT  requires  certain  inputs.  The  array,  containing  bias 

'  ^ 

limits  for  each  of  the  measurements  and  the  S„  array,  containing  the 

'  \ 

precision  indices  for  each  of  the  measurements,  must  be  input  (as 
fractions)  within  the  UNCERT  subroutine  (lines  20100  to  20600).  The  B„ 

A 

values  come  from  expert  opinions,  manufacturers'  data,  or  prior  testing 

of  the  instrumentation  against  a  known  standard.  Sv,  values  can  come 

k 

from  manufacturers'  data  or  from  prior  testing  of  the  Instrumentation. 

If  Z  transducers  are  used  to  determine  the  same  measurement,  then 

1  transducer  ) 


The  resulting  value  should  be  input  in  the  S  array. 

Finally,  proper  array  sizes  must  be  used  in  the  main  program,  the 
performance  subroutine,  and  the  UNCERT  subroutine.  When  FORTRAN  IV  is 
used,  as  is  the  case  in  this  example,  array  sizing  must  be  done  in 
COMMON  statements.  Although  COMMON  statements  are  not  a  requirement, 
they  do  facilitate  the  programming. 

We  need  two  array  sizes  in  the  program,  L  and  N;  L  is  the  total 
number  of  performance  parameters,  and  N  is  the  total  number  of 
measurements  on  which  the  calculations  depend.  In  the  UNCERT 
subroutine,  we  must  size  the  BCJ,  SCJ,  and  ANSWER  arrays  using  L.  DCODM 
is  a  two-dimensional  array  holding  the  influence  coefficients  and  must 
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'.‘  *■•  '> 'A  ■'•  '>  ■.■•'•T*.* ' 


be  sized  DCODMfL,N).  The  M '  j.  array,  MP,  must  be  dimensioned  in  a  REAL 
declaration  using  N.  We  size  the  M  array  in  the  COMMON  statement  using 
N.  In  other  COMMON  statements  the  C  and  CPRM  arrays  are  sized  using  L. 
The  DIMENSION  statements  and  proper  sizing  statements  are  seen  in  lines 
19300  through  19900  in  the  subroutine. 

In  the  main  program,  we  must  size  the  M,  C,  and  CPRM  arrays  in  the 
COMMON  statements  as  was  done  in  the  UNCERT  subroutine.  In  the 
performance  subroutine,  the  C  and  CPRM  arrays  are  again  dimensioned  in 
the  COMMON  statements.  We  must  write  the  COMMON  block  for  the  M  array 
(in  the  performance  subroutine)  as  shown  in  lines  36200  through  36^00, 
since  the  measurement  values  will  be  input  into  different  variables. 

This  special  COMMON  block  could  appear  in  either  the  main  program  or  the 
performance  subroutine,  but  it  should  not  appear  in  the  UNCERT 
subroutine.  The  COMMON  block  name,  DSPAR  here,  must  be  consistent 
throughout  the  main  program  and  the  two  subroutines.  This  is  true  for 
all  COMMON  block  names. 

In  the  call  to  the  UNCERT  subroutine  (line  13700),  two  integer 
values  must  be  passed:  L  and  N,  After  all  requirements  have  been  met 
and  modifications  completed,  the  program  can  be  run  to  yield  uncertainty 
results.  The  program  outputs  the  nominal  performance  of  the  test 
article,  the  equation  for  uncertainty,  bias  limits,  precision  indices, 
and  final  uncertainty  values  for  the  performance  parameters.  Other 
output  statements  may  be  added  by  the  user. 


u 


r.-i 


134 


s'. 
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In  summary,  the  application  of  this  subroutine  depends  on  several 
activities.  First,  a  computer  program  must  be  used  to  find  the 
performance  of  the  test  article.  Then  the  performance  calculations  must 
be  contained  in  a  separate  subroutine.  Finally,  certain  assignments 
must  be  made.  These  are  listed  below. 


1.  The  various  measurements  must  be  assigned  to  the  M  array  in  the  main 
program. 

2.  Performance  calculations  must  be  placed  in  dummy  array  DC  at  the  end 
of  the  performance  subroutine.  Proper  statements  must  be  added  so 
that  these  values  are  placed  into  the  C  or  CPRM  array,  depending  on 
a  flag  value. 

3.  The  M,  C,  and  CPRM  arrays  must  be  placed  in  COMMON  statements  in  the 
main  program  and  the  two  subroutines.  The  COMMON  block  for  the 
measurements  contained  in  the  performance  subroutine  must  be  input 
in  a  special  manner  (see  previous  page). 

1)  .  All  required  inputs  must  be  made.  The  BMK  and  SMK  arrays  must  be 
input  as  fractions  in  the  UNCERT  subroutine. 

5.  All  array  sizing  must  be  accomplished.  Some  of  the  sizing  must  be 
accomplished  in  the  COMMON  statements. 

6.  The  UNCERT  subroutine  must  be  called  with  the  two  integer  values  that 
represent  the  number  of  performance  parameters  and  the  number  of 
measurements,  respectively. 

After  running  the  program,  the  user  can  program  additional  output  as 


desired. 
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UNCIHTAINU  CALCOLAUQMS 

STSIEM- 

aUhHUOaHS  6900 
PNCOHAMMIN- 
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PURP05F- 
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functional  area- 


usage- 


CALLING  SEQUENCE- 


CCHHGN  STORAGE-  ' 

HCAURAY  ECLOING  NQHINAL  RLAStKEO  9ALUES) 

CCAHRAY  HGLUING  NQHINAL  CALCULATED  VALUES) 
CFRH(AFfAV  HCLOING  VALUES  CF  THE  CALCULATED 

PARARETERS  AFTER  A  CHANGE  IN  THE  HEASUREO 
VALUE  HAS  BEEN  HADE) 

EXTERNAL  RLFLPENCES- 

GKPRUCGENERAL  PATHIX  PRODUCT) 
lAPERFCTESf  ARTICLE  PERFCRHANCE) 

FILES  USEU- 
INPUT/OUTPUT- 

ThE  bias  LIHITLUHK  iN  x>  AND  THE 

PKEL'ISIUR  lAOEXCSHK  IN  X)  FIR  EACH  OF  INC  PEASUREHENTS 
HLST  be  INPll. 

BIAS  ANU  PHECISICN  ERRORS^  AtCRC  MllH  THE  aNSHER* 

AHE  CLTPLI  fCR  EACH  Cr  the  CalCllATEO  pArAhETErS. 

these  are  given  in  percent. 


SUBRUui INe  UnCERT(L>N) 

DIHENSION  UCar;H(64*20  ).I)C.(b4)#SCJ(b4) 

UXNENSiUN  UELTAH(20)»UtL  T AC (b 4 > ^ AKSHE  R C 64 ) . BRK C 2o )«  SMK( 20) 

REAL  M*MPC20) 

CONHUN  /U3PAH/  M(20) 

CaPMUN  /NQPVS/  C(e4) 

COHNGN  /PHIHS/  CPNM(b4) 

UAIA  BRK/.COie».CC2n».OOS«0.r  .0021S*  .OC  |€  ..00  lb»9  *.00  16»2*.  00  3» 

A  .002».002S/ 

OAIA  SPN/.aOO47#.COOl^».CO2Sr.O0CCI6#.CQ01C*.O03S4»*00M8>.O0ir7* 
A  . GO 12S» .000  66, . 0 0 086 , . 002 S ^ . 0 Q2S, . 0 028 9 , . 002 69  *  .00  1/7, 

B  .0004, .0004, .000^4, .OOI2S/ 

DO  lU  N=1,N 

lF(H(N).tQ.Q.)Pt)S)-.Ot/l,Cl 

NP(K)3l.01*H(M 
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00  ^=1.N 

^1400 

UUM=M(K) 

ZlbOO 

2lbua 

C 

21/00 

CALL  t*PLH( 

21dU0 

C 

21^00 

H(K)=0(JM 

22000 

c 

22100 

ir(H(K).La.o.)kHiU(t:.n) 

22200 

15 

rOKMATCDl.'PRCGRAM  STCPPCC  StNCC  A  HLASLRED  PARAMCIER  IS  2LNQ') 

22500 

IFEHClEI.EQ.O.IsraP  nil 

22400 

c 

22500 

00  1/  J=1.L 

22600 

ir(C(j).Ea.a.)iiRirc(6.t6> 

22/GQ 

16 

FChHArElA>'PKCGRAti  SIEPFEO  SIACC  A  CALCULAIEO  PARANEIEfi  IS  ZERO*) 

22600 

IFECE  JKEfl.O.  ISICP  1112 

22900 

1/ 

C  ONI IRUL 

25000 

C 

25100 

0£LTAMEK)  =  £HP(IO-H(K)1/H(R) 

25200 

00  10  J=1>L 

25300 

D£LrACtJJ=ECf«F(J)-C<JJ>/CEJ> 

25400 

OCaOM(J.K)=OCL1ACE^)/CtL1AN(Kt 

25500 

16 

CONTINUE 

256UQ 

C 

25/UO 

20 

coNriNoe 

25800 

C 

259G0 

DO  70  JsI.L 

24000 

00  70  K:1.N 

24100 

OCOOM(J.N)=(OCOOMCJ«K)><*2 

24200 

/o 

CCMINLE 

24500 

C 

24400 

00  75  K=I,N 

24500 

ei'K(K}=apK(K)>*2 

24600 

SHK(K)=5RNE«)**2 

24/00 

/5 

CONTINUE 

24600 

C 

24900 

CALL  GNPRO(DCOON.0MK.BCw>L»N.l) 

25000 

CALL  GNPRO(OCODN>SHK>SCO»L>N.  1) 

25100 

c 

25200 

00  77  J=1.L 

25500 

tfCJEJI-SURECUCJE J>) 

25400 

SCjE  JI-<ORICSCJ(  jn 

25500 

ANSWLH<J)=bCJ(J)«2.0«SCJEjl 

2  5600 

rf 

COM  INLE 

25/00 

c 

256U0 

OU  7b  J-l.L 

25900 

ANSMLf<EJ)  =  ANS»llRCJ}*lCC. 

260u0 

/a 

CONIINIE 

26100 

C 

26200 

HNllLEC.bQJ 

26  5uO 

RK  1  TL{t,.b7  » 

26400 

C 

26500 

NHlUCb^bOJ 

26ouO 

NKlILEb.bS) 

26/00 

00  85  J=1.L 

268G0 

MNITL(b.(;61J>UCJ(Jl>SCJ(J>>ANSIiEKEJ) 

26900 

65 

CONT INLE 

2/000 

C 

2/100 

C 

2/200 

bO 

EUNNAKIHI  ) 

2/500 

61 

FORHAK  1  X>  10(  E  1 1  .  A,2X)  > 

2/400 

65 

FCFNAIEIX.SX.'FASAPEIEF  •*.7X.*B1AS  L I P  I  T *.  1 1  X.'PRl C  I  S ICN  INDEX’. 

2/500 

A  bX.’LNCLNIAINIf’./E 

2/600 

bb 

EaRMAI(lX.I0X.I2.9X.ln.A.|2X.EIl.A.lnX>L9.?) 

2//00 

b/ 

F ChM AIE 1 X. luX .’THE  ECIAIICN  USED  IC  OEIEbMlNE  iHE  L N C 1 H I  A  1 N II E S * 

2/ioo 

A*  IS-.//.  1  lX.’L  =  H  ♦  2  •  S’.//.11X.*B  IS  DEEIMD  AS  TliE  BIAS’ 

2/900 

B’  LINIT.  AND  S  IS  OEFINtU  AS  IHE  PRECISICN  I  NOE  X’ . / /. 1  1 X. 

26000 

C’THE  VALLE  CF  195(2)  CORES  FROR  IHl  TBD-IAILEO  SIUDLMS  I’ 

26100 

0’  OISTNIBUI ION’) 

26200 

C 

26500 

RE  lUKN 

26400 

END 
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zaiuo  c»«*««**»»*«**««*»«****»**»*»***»» 

ZdbuO  C 

^a^UO  SUUKQUIINE  GHPNO(«.0>b.b,M,L) 

EdaOQ  UIMhSlUN  J.h(l) 

^ayuo  iR=a 

2*000  ia=-M 

29100  oc  10  a^i>L 

29200  IK-IK*M 

29ioa  oc  10  j=i>a 

29b00  1R-IU*1 

29b00  JI=J-N 

29600  lU^IK 

29700  RtlR)-0 

29aou  DC  10  i=i>a 

299U0  JI=JI*N 

iOQOO  IR-ie*l 

Juloo  10  H(  IR  1  =  K( 1R)*«( Jl l•U( la  ) 

10200  RCIUNN 

iOlOO  LRD 

J0600  c***************************  . 

loboa  c 

iObOO  C  NAHL  - 

10700  C  lAPERf  5/2*i/e2 

loaoo  C  OAhRCTI  CCklKirUGAL  CCAFRESSOR  PLRl CRHAACt 

10900  C 

llOOO  C  &7$rLN  - 

11100  c  laa  170/1^6  (aaia  corfcuri 

11200  C 

lllOO  C  PRUGRARMtRS- 

llbOO  C  HARK  S.  CRR  AfkAL/POIl  FURCUC  UMvCfiSlir  CC'OP  SIIOCAI 

11^00  C  CHRISIOPHCR  A.  BOCOICKtR  OSAF'  ACA0£Nr/CS*20 

llbOO  C 

31700  C  PURPOSt  - 

iiaoo  C  FHIS  RCQLLL  OGCS  COHPFCSbCR  PCRFORMANCC  CALCULAUONS 

11900  C  SPlClFltO  UY  GARRtll 

12000  C 

12100  C  rUNCriUNAL  AREA  - 

12200  C  RAIR  CGRPUUH  CCFPRtSSCR  PLRECRMARCE 

12100  C 

12A00  C  USAGE  - 

12S00  C  CALELO  dV  IFPERF  FLR  EACH  CAFA  ACCUISIIICA  POINT  THAI  IS  |C 

12b0a  C  BE  PROCESSED. 

12700  C 

12800  C  CAlLIAG  SEAUEACE  - 

12900  C  (HI  EU  -  MEASURED  PARAMETERS  ENGINEERING  UNITS 

11000  C  (H)  PC  -  PERFORMANCE  CONSIANIS 

11100  C  (R)  mode  -  OAIA  ACOUISIIICA  R COE ( I =NCN> P^ST AT IC. 1  =  IR AN. A- ABOfc  I  > 

11200  C  (R)  BARON  -  RARCMETNIC  PRESSURE 

llluo  C  (Ml  KKK  -  BUFFER  CONIAINIAG  COMPRESSOR  CONFIGURATION  CODE 
llAOO  C 

IISOO  C  COMMON  STORAGE  - 

11600  C  NONE 

11700  C 

liaOO  C  EXTERNAL  REFERENCES  • 

11900  C  ACGS  .  AIAN  .  COS  .  CPTBL  •  ITHAT2  .AlCG  .  SIN  .  SORT  ,  TAR  .  UVAr 

lAOUO  C 

lAlOO  C  FILES  USED  - 

14200  C  DAIA  BASE  RTCORD  NUMBER  20  •  ENGINEERING  UNITS*  M 

14100  C  UAIA  BASE  KECCRO  NUMBER  2S  -  PERFORMANCE  CCNSTANIS*  R 

14400  C 

14S00  C  IhPUI/OUIPUT  - 

14600  C  NCNE 

14700  C 

14800  C*«*«**«*  .  •••••*.••••••  . 

14900  SUBRUUIINE  lAFtNF 

ISOOO  C  DIMENSION  ALPH (1 8 1 . C(  6 4 1 .CPRM (64  1 

ISIOO  C  LINE  ABOVE  REPLACED  MIIH  LINE  BCLOH  TO  UEEETE  C  (DIN  AI  11/0).  AnO 

1S2U0  C  CPRM  (DIM  AT  >380). 

ISIOO  DIMENSION  A L P H ( 1 B ) . DC ( 6 4  ) 

1S400  REAL  J*N*MU2S.NCCR*N(20) 

ISSOO  COMMON  /KONSI/  AI9  .C019  *420  .C02C  .R20ME  .H2aSH0*G> 

ISbOO  A  J  .RG  .rsTO  ,PSIC  .N2S  .A2S  .BEIAbL. 

IS700  B  CU2S  >CDlHN  .IhRCAI.ARU  «ANl  .COEE  .RLE  . 

isaoo  C  ALP9LE 
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jbyoo 
iboua 
5bl  UO 
ibZOO 
ibiOO 
30400 
ib;>oo 

SbbOO 
30700 
50800 
30900 
3/000 
3/100 
3/200 
5/300 
3/400 
3/500 
5/000 
3//00 
3/800 
5/900 
38000 
58100 
30200 
18300 
58400 
30500 
18600 
58/00 
38000 
18900 
19000 
39100 
39200 
19500 
39400 
59500 
39600 
39/00 
59800 
39900 
40000 
4C100 
40200 
40300 
40400 
40500 
40600 
40/00 
40800 
4O900 
41000 
41100 
41200 
4  1  300 
41400 
41500 
4  1600 
41/00 
4  1800 
4  1900 
42000 
42100 
42200 
423O0 
42400 
42500 
42600 
42/00 
42800 
42900 


C 

C 

C 

COHNQN  /USP4H/  P5t9  ,1120  ,N*V  .6  ,1150  ,PS2ShD,PS25 

4  PS29  ,fS10  ,PI29  ,PljO  ,p5IHKU,PSlhHL,P5£XO 

a  PStXl  ,PSLt  ,«OIff  ,a1PHIG,SH,225 

c 

COMMON  /FARAM/ 

A 

a 

C 
0 
e 

f 
0 
H 
1 
J 

c 

COMMON  /FLAGO/ 

COMMON  /NOMVS/ 

COMMON  /PRIMS/ 

C 

OAIA  ALPH/U. 0,00.0,0. 3, 14. 1,0.0, .02, 5.0, .02, 10.0, .02, 20.0,. 02, 
A50.0,. 02 05,40.0, .026,60.0,. 039/ 

500  rcNMAi(5x,(:i2.5) 

GC=32.2 
PI-3.14159 
f  A=0.0 
HC=0.0 

C***  CALCULAIL  HSIU  bAStO  ON  1510 

CALL  CPl6L(ISI0,H51D,FFSIC,FA,HC,SH,l,FC,CAMMA,CP,kA,CX> 

C»*»  CALCULAU  GAM2,H120.  ANC  RP20  AT  1120 

CALL  CPIHLl  lT20.Hl20,FP20,f A,HC,SH.l,fiG,GAF2,CF,NN,CX) 

C 

C 

C***  IIERAIE  fCR  AM19 

U1-.5<(GAM2*1.C)/(GAm2*I.O) 

G2=.5*(GAM2-1 .0) 

NDIMl^O.O 
N0IM2  >0.0 
AM19=.5 
AMNA^l.O 

HDIM-UA/'SORI (TT20)/(PT19*C019«A19) 

5  CONTINUE 

C 

UUMl -GAM2*CC/NG 

IFLOUMl.LT.O.O)  MRIIL(6,500>  DUNl 
IFLOLMl.LT.O.O)  STOP  1001 
C 

i(0IM3  =  5aRHGAN2,GC/RG)«AM19»I<  1.0*G2«(  AM19«*2))«*(-Gl  )> 
ll(A8S(MClMl-hUrM).L1..000l )  GCTC  / 

CALL  IIRAT2CNOTHl,AMNl,kOTM2,AMN2,HOTM3,AM19,MOTN,AMN4) 

GOIO  5 

/  CONTINUE 

C««*  TAIILE  LOOKUP  FOR  ALU520 

ALL520>U9AH(ALFH,ALPHIG1 

C 

C 

C**«  CALCULATE  PSPI19 

PSPTI9x(i.0»g2«(AM19..2)).*I-(0AM2/ICAM2-1.0)») 

OPIICV=PII9*ALL520*(1,0-PSPT19) 

PT20=PI19-0P1IC9 

01GV=0.0 

ALFH20>ALPH1G*FI/180. 

C***  ITERATE  FOR  AM20 
NOIM1=:0.0 
N0IM2=C.C 
AN20-.5 
AMN4=1.0 

M0IR-NA9>S0RT(  II20)/CF  I2C>C020aA20<C05(  ALPH201 


AM19 

VU20 
V2CSH0 
U25 
PRSI 1 
9U?5A 
ALUS29 
ETAU  30 
CPU 

PRECLT 

PNT  r 


.ALGS20, 
.9M20  , 

.81 IA2S. 
.AM25  , 

.DEVI  , 

.sr 

.AM29 
.FRSI 
■  CPL 
.  AMLE 
•  IFIT 


PI20 
U20M 
RM2SHD 
ALPH25 

E  TAT 
OLTGIA 
EIA029 
EIAST 
EFFU 
ALPHLL 
ETAII 


, ALPF20 
,B£ IA20 
.Of  LF 
.a£ IA25 
n2|i1 
,DL  I  CT 
,PR£C50 
,£  TaFCL 
,£f  f  1 
, AINCO 
,£TAAO 


,  AN20 

,RP20 

,OELHC 

,RN25 

,hRlSE 

,FI25 

,AL0S3O 

.PSIHR 

.PRECU 

.CAPRAI 


,920 

,AN2ShD. 
.CLICIC, 
.PRIII  , 
,PAVBF  , 
,PMEC29. 
,AM30  , 
,PIIRR  , 
,PR£CL  . 
,HPBF  , 


IFLAG 
C(64  ) 
CPRM(64  5 


US  AKA 


k  iOUO 
4i  100 
4  i^UO 
4  iiOO 
4i40U 
4  i:>00 
4ib00 
4  J^UO 
4  idUO 
4iV00 
44UU0 
44100 
44^00 
44  500 
444U0 
44^00 
44bU0 
44  ?G0 
44  600 
44900 
4;3UuO 
4t>l00 
4'>^U0 
4tJ  iOO 
4‘j4U0 
4'aSOO 
4bb00 
4‘J/OO 
4^800 
4!)9Q0 
44000 
4olU0 
4b<;00 
4bi00 
4b400 
4b^00 
4bb00 
46  ?C0 
46600 
4b900 

4  ;aoo 

4/100 
4/200 
4/500 
4/400 
4/'>00 
4/600 
4  //OO 
4/dGO 
4/9J0 
48000 
48  100 
48200 
48500 
48400 
48‘>U0 
48600 
48/00 
488U0 
48900 
49000 
4 '/I  00 
492U0 
49500 
49400 
49'>00 
49600 
49/00 
49800 
4  9900 
*>0000 
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IS  CONTlNLi: 

c 

c 

WCrM5-SCf.TCuilH2*Gf/HGl*AH?C*((  1  -  0  ♦G?  •  C  4M2  0  ••  2  )  “GI  )  ) 

ir ( Ac^CMCT H 5-MUT ^ >-L 1  - -CCCl )  GlIL  1/ 

CALL  irRAr2CUCTMl*AHNl#ti0TH?#AMN2*UOlHl»AN2O»MUTH.AMM4) 

Gora  IS 

1/  CGHTlALt 

C 

C 

DLH2  =  1  .Q«G2*( AH2C  ft«2 ) 

IKQUM2.Lr.O.O)  «RIT£t6^S00l  0UH2 
IF<UUM?,Lf .0.0)  STOP  l002 

C 

92 AT20=A«20/SflfiTCl.O*G2*( AH20*«2)) 

C 

OLH  5=  u AH2 ‘GC •HG* I  1 ?0 
IMLUX5.L1.C.0)  hh  ITf  C  b.SOO  )  UUhS 
If  ([JUM5.LI.O.O)  STOP  1005 
L 

AT2C<SwhUCAI>^«GC*ftG*1l2C) 

/20-92AT20«AT?0 

VU20=V20*<SIS{ALPH20)) 

Vr20-V20*(CLSCALPH20)J 

U20Mi?.0*PI*N*F20Mt/(bU.»l2.) 

HErA20=4rAN(( U20M-VU20 )/VM2Q ) 
HP2OsA82C*(CCS(ALPH2O)/CCS(0ilA2O)) 

C 

DUM4  =  C((PT20/PS2SHCi*«<G?*2,/GAM2))-1.0>/G2 
If  (CUH4.L1.o.O\  hf  ITE(b»SOU)  0U|t4 
If (0UH4.LT.0.01  SrCP  1004 
C 
C 
C 

4>»2SH0*S4H1  ((  ((PT20/PS2Sh0)**C62*2./GAP2  ))-1.0)/G2) 

C 

DUMSi I .0*G2*(AH2SH0**2) 

IMGUPS.L  1.0.0)  hflKIb^SOC)  OUMS 
U COUMS.LT.O.O)  STOP  1005 

0 

9AI2SU2A^2Sh0/SCHT(UC4G2*(  AN2SH0**?)) 

V?OS>tD-yAT2SU*AT20 
VU2SH0-y20Sh0«SlN( ALPH20 ) 

9H2SHO=V20SHO«COSr ALPh2C ) 

U2SHU=020M*h205HC/H20Mt 

Ut  f A2S-4I AN(<  U2bHD-VU2SHD )/VM2SM0) 

hP2SH0=AI'25hO*CCS(ALPf2O)/CCS<HE  lApS) 

i<2O5H0-KN2SH0«S0RT(  CAH2*CC*RG*TT20/CU0^G2«C  AH2SH0»*2>)) 
C***  calculate  hT2C  a^ 

CALL  CPTbLC  ir20»HT20.PP20#f A.HC»SH#|»RCrGM2#CP2»NNrCX) 

C**«  CALCULATE  HT50  AT  1150 

CALL  CPTbHlT30»Hl50#fiF3C»fA»HC»SH»l»RGrCA^5FCP3»AA-CX> 

C 

C 

OELH=HT3C-HT2C 
0ELHC=DELH*  TSTC/TT20 
HT30C^HSrU*DtLMC 
C***  CALCULATl  TTiOC  AT  HTiOC 

CALL  CfT»L(TT5UC»HI50C#RP30C»fA»HC»SH»2#RG#GP5C»CP3C»44»CX) 
OTTC^TT  50C-TSTU 
DL  TC  IC=DTTC/TSTO 

C»**  CALCULATE  CP  AT  ( T T iOC ♦ I T20 )/2 .0 
T  JCSU2  =  CITi0C»IT20)/2.0 

CALL  CPrOL(T3CSn2»L#NP#fA*HC»SHrURG«CAM5CS«CPlA/#NK»CX) 
U?5=U20M*H2S/R20HE 
C»**  ITLRAIE  rcR  TS2S 
TT3CHl=0.0 
rT5CM2=0.0 
IS25=TI5C 
TSN4  =1130-10. 

C 


2 
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10100 

‘>o^oo 

10300 

lOOOO 

10100 

10600 

10/00 

loaoo 

10S/00 
IIUJO 
1 1 100 
11^00 
:>1  300 
110  00 
11100 
lluoo 
11/00 
iiaoo 
11/00 
i<;ooo 
i<!iaQ 
i^<;oo 
1<!300 
12000 
12100 
12600 
12/00 
12000 
12900 
13000 
13100 
13200 
1  3  500 
13000 
1  3100 
1  3600 
13/uO 
13800 
1  3900 
lOOUO 
16100 
1*200 
1*  loo 
1**00 
1*100 
1*600 
1*/00 

I  *  a  0  n 

1*9o0 
110  00 
1110  0 
112  00 
11100 
11*00 
1110  0 
11600 
11 '00 

I I  a  0  0 

11900 
16000 
16100 
16200 
16  5  0  0 
16*  00 
16100 
16600 
16/  00 
16800 
16900 
1/000 
1/100 
1/2C0 
1/300 


VU21i(CPl«V«QI  IC*OC«J*L2CP*VU20)/O2‘j 
PCA=PS21*C021. «?1 


CQMIMJP 

CALL  CPTBL(1521.t.HP.f*,hC,SH.1.6i,GAp21.CP.*N.CA) 

If (u*H2l.ia.0.0)  IIlP  1121 

801=1902  1*«.  J/(C.*H21*GC*0o.ir5OC) 

03=.1«(0AP2l-l.C) 

AA=&3'(G5«*2  1*0(11 

(302  =  1  .  C“(0  AM21-1 ,0  )  *0  0  1 

l}LM6=i(G/(GAM?1*Gr  1 

I  (  1 oUHb . L 1 . C . U 1  Hhllllb.lOO)  OUOb 
It  ILUft.L  f .0.0)  Oltf  |C06 

CC=-(3Ul-((221.i(AV*iCKt(II30)/PCA)*SeKI(PG/(GAH21*bC))).*2 

ooM/:r(ia2*.2)-*.****cc 

IK  UUP/ .L  r.C.O)  MSIlflb.lOO)  UUH/ 

IKCUP/.L  I.O.U  1  5  Uf  100/ 

DOMt!  =  (-(3H2*SCKT(CLN/))/(2.***» 

IKGUMd  .1  I  .0 .0)  66111(6*100)  OUPB 

If tuuMa.Li.0.0)  Slop  1008 

AP21=1*KI(C-Bd2*SCHl((HB2**2)-*.«***tC))/(2.***)) 

rr3CP5=rS21*(1.0«C3*(AH21..2)) 

0  =  0*1. 

iru.OI.*0.)6Nin(b>993)tI3CNS.II5aC 
iriO.GI.AO  )GUIU  2/ 

fURM*r(H.*It3CM3='.fl1.$.5X.'tI30c=*,fl1.9,5*,*LOCP  EAIIH)') 

IK  AdSl  iriCHl- 11300.  I  I..001  )  GOTO  2/ 

IKC.Gt.cO.  in  3C  65  =  11)03*. 21 

CALL  llRAr2(II)C61.rSA)*II3CH2,|^N2,rT3CM3*TS21*n5  0C«IS**) 

GJIU  21 
COM  INoL 

U)A212=1.0*6)*(AP21*.2) 

rjLM9  =  ,jA621*GC/(HC*lllO) 

I  K  lyOH9.  L  I  .  0 .0  1  66111(6*100)  0UP9 

II  (L0P9.L  I.C.l)  Slip  1009 

1aC6RI=SCRIlCAP21*6C/((K*ll50)) 

l)CM  1  J  -  (  ()  I 

I  K  LKN  1  0 .1  T  ,  0 . 0  1  (tf  I  |(  (  6  *100  )  UURlO 
IKLOMC.LI.O.P)  SILP  1010 

AL  PM21=A  lAM  SCRKlibl  )*(6iA2 12  )*(P121. CC21.A2 1/(221* »aw))*SC6GR|) 

JOMll-l,  AP<'1*i,C*hi:*ll(uC 
1 K CUH 1  1 .L I .0 .0  )  ■61l((6.100)  nUPIl 

II  (LORI  I.L  I.C.O)  SILP  I6l  I 

GOP  12  I  .  0  *6  S*  (  AP2S.  *2  ) 

I  K  COP  I  2  .L  1  .0 .0  )  *1111(6.100)  nuM/ 

IK.IOPIC.LI.O.OI  Slop  1012 

V21  =  AP21*S6l(I(6*P21*CC*66*lI)OC)/Si.Rl(1.0*G5*(«M21**2)) 

VN21 =921*C0S( ALPH21) 

Bt  l*21=*IAi9(U21/*P21-I*M  *1  Ph21)) 

HP21  -  *P2S*CCS(  Al  pi’21  )  /(  (  n/.j  I 

PI2 1 -  PS2  1  *1(1. u*g5*(AP21**.'1)**'GAP.'1/(LAP2  1-1. 0)11 


;SA  l’A-  TR-S4-  ') 


5  /buO 
‘)//uO 
‘>/dOO 
5/900 
5OO00 
5tilOO 

soeoo 

5oiU0 

5d400 

5ti5oa 

50600 

50/00 

50000 

50900 

59000 

59100 

59200 

59S00 

59(i0g 

59500 

596U0 

59/00 

59000 

59900 

bOOOQ 

60100 

60200 

boioa 

60400 
60500 
60600 
bO/UO 
bUduO 
60900 
6  1000 
61100 
61200 
61  JOO 
6  1400 
6  1500 
falbOO 
61/00 
6  1000 
900 
62000 
62100 
62200 
b2i00 
62400 
62500 
62600 
62/00 
620UO 
62900 
6i000 
6il00 
6  i200 
6}SUQ 
6i4  00 
bi500 
6  J600 
oi/UO 
6  iaOQ 
6  i9u0 
64000 
64il60 
6%2U0 
64  iOO 
o^<iOO 
64  5U0 


I  T  i  -  H  r25  /  }’  1  19 
P:yST  L  -P52'./H1  19 
OE  V  !-■  Bl  I  92  5-(  ':l  r  AtJL*H  /  I  fcO  •  ) 

I’  5  I  =  9L25  /U2‘J 
PMi-9M25/U25 
H1N-HT20 

CALCULAIL  H  i  S(  N ( HP  I  5 f N  )  hH£kt  RP I S£ 6  =  HP 2 C ‘Pb F T  1 
HP15LN=HP2u*Pflni 
C 

DUM2iK--nPiSt  N 

IFCCUPPiH.LL.O.G)  bhnE(t>500)  DUI>2ih 
If CDUH2iR,LL.U.O )  STOP  1025 
C 

CAtL  CPTBL(T15l6*HlSfb*bPlS£N«rA*HC#SI-»3«bC#CANlSE»CPlS£*46»CX) 

LTAl-ChI5tN-Hl6)/(HtX-HlM 

i(2wl=VH25/lW2  05HU*CGS(btlA25)) 

hfi  ISl =MtX-Hlb 

PArfBf=PS25*Cl.C-(U25«*2)/(l6.*CC*«L*TTj0)> 

HU25  =  (2.25/£-0  )•( n  30**( 1-5)) /tri 30C* 19  0.6) 
HHC=PAVBf.l44./(225*RC»IIiO) 

C 

t 

fUN-HhU«U2  5*h2  5/C12.*Mt2  5*CC) 

CHBF-O.0il94/(HF.N**0.2) 

C 

C 

CM-hHa«(U25«*2)«tMbK«(H25**5l/(l/26.*2.«GC> 

HPbF=C0‘6»2 -• f 1/33000- 
OFOtL=55O.*HP0f/( J*MAV ) 

HAtHUl=HbI St-OFOtL 
HAtFU2=htX-0FDtL 

/U25 A s( GC»J*H At ROl *020 H*9U20)/U25 
C 
C 

Sf-VU?5A/(U25-VH25«IAK(BtTABL*Pl/l80.)) 

CALCULAIL  lAthtCIUf  6C2) 

CALL  CPTBLC  TACHa#HAEHC2»RPAeK0»FA»HC»$h»2»fiG»GANAtR^CPA£H,NK,Cx) 
OLIOTA=( rAtR0-rr20)/TT20 
ULrCr^(n5a-TT20)/TT20 
C*»*  CALCULAFE  GaMCFC TF 30^ F F2C  )/2  .) 

T3l2AV-(FT30»n2C)/2.C 

CALL  CPFtlLCF5T2A9»H52A9-hPi2A9»FA»HC*SH»l*FG»GAHC»CPi2A2»^N«CX) 
PF25  =  Pri9*((UL  For •£ T AI ♦! .0)*»(GAMC/(GAPC* I.O  )) ) 
PHtC29i(PS29-FS25)/(FT25-PS25) 

ALUS29=(Pr25-PT2g) /(Pl25-P4?5) 

€*••  CaLCUlAfI  GAPDC C FF  !0» F525  )/2.0  ) 

Tir25A=(lTJ0*TG25)/2-0 

CALL  CPTaLCl3l25A.H325AV^HP325ArFA#HC»SpFl*6G*6A^D-CF3254»fc6-CX) 
C 

DUM15-  (2,  /(  GAMD-|.))*(((PS29/PI29)««C-C6AH0-I.)/CA»<D))“1.) 
IF  (OUPl  S.LT.Q.O)  6MFt(6»500)  0UP13 
IF(UUH13.LT-0.0)  STOP  tOli 
C 
C 
C 

AM29  =  S0HT( (2. /CGAMU-1 .0 ) )•(£(? S?9/PTZ9)»*C-CGAMD-l-)/CAHD))-l-)) 
IS29  =  TI30/(I.C4((GAf'C  -  I  .  O  ) /?  .  >  •  C  A0Z9  •  *2  )  ) 

tTAD29=((CPS29/PS25)**((CAH0-|.>/GAMO))-i.)/((FS29/lS25)-l.) 

PHLC  J0=(PSiC-PS25)/CPF25-FS?5) 

ALUS jO=(Pt25-PI J0)/(PT25-PS25) 

C 

DLHI4  =  (Z-/(CAPn-I.))*l(CFS30/pT50).*{-(C4HO-U)/CAHO))'l-> 
IF  CO JHl 4.LT-0.0  )  hRiTL(b»5U0>  0UM14 
If  (LUm-LF-0-0)  STtP  1014 

C 

AHiU  =  S.-tm(2./(GAHlJ-l,J)«(((P55O/Pf5O)««{-(GAHU-l.>/GAP0))-i.*> 

I  !  J  -  F  I  3t /(  1  .  ♦  •5«  (GAFC-1  -  )  •(  46^0  ‘‘Z  )  ) 

LT4Ui0-(((PS3O/P525)««((CAMO-l.)/GAHO))-l->/C(TS30/FS25)-l.) 
PHFT^PF  ^G/PTl^ 

PFST -PS  5C/P  I  ig 

£•••  CALCULAU  UStMCHPISFl)  bhLHt  HP  I  SI  T  =  R  P20«  P  H  F  T 
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mV.S 


r  S  A  [•'  A  -  r  R  -  8  4 


b4b00 
64/00 
t>4d00 
64900 
bb  JUO 
b^lUQ 
0  ^  00 
b'jiUO 
b64  OO 

b  ‘j  6  0  0 

bbfaOO 
b^a/00 
bboOO 
b‘>9U0 
bbOUO 
bbl  00 
bb^UO 
bb  500 

bb4  qO 
6  b  6  0  0 
bbbUO 

bb/U3 

bbaoQ 

bb  9u0 
b/OOJ 
b  MOO 
b/^OJ 
b/  5UO 
b  /400 
fa/’igO 

b/bUO 
6  //OO 
b/  dUO 
b/VOO 
bOUUO 
bOl  UO 
bd^UO 
bo  300 

o  04  0  J 
bO'>00 
bObOO 
bO  /  uO 
bddOO 
bd9oa 
bVO  JO 
bVlOO 
bV^UO 
b9  ioO 
694  UO 
b9‘jU0 
bVbOO 
b  9/  uO 
bVdUO 
69900 

/  jouq 

/Ol  00 

1  OtT  0  J 

'0  iOJ 

'  U  4  «  0 
/  U  ■>  0  0 
/  0  b  U  0 
'0/00 
/  J  -10  0 
/  0900 
MOUQ 
/  1  I  JO 
M  ('OO 
/ 1  luo 
/  1400 

'  r^oo 

/ 1  b  j  J 

/ 1  Mo 


HPXSLT=RH;'0*PPTT 

c 

U0M.'4H^hP  ist.  r 

lMLU»>24K.LL.a.O)  hhllt(6#^00) 

:i  CDUM^4R.LL.0.0  )  SIGP  lGi?4 

c 

CALL  Ci'rbLCT!^tNT.-Ml5rM*RPIStl*FA,HC#Sh#i*RC^LA>'ISr^CflS£I#h4»Cx> 
C***  CALCULAIL  HICLKiCHPISLSl  bHLRt  RP I S t S =H p 20 •  P HS F 
HPiStS^HP2C*PRi)T 

C 

0LM26H=HPISES 

If  (CU^^^K.LL.O.O)  bhUt(6*SOO) 

If (OUM2bR.LL.O.O>  5IQP  l02^ 

C 

C 

C 

CALL  CPTUL(riSLNb#MlSfNS«HPlSLS#rA»HC#SH*5»RG#CAMlSS»CPiS£S#MK»CX) 
LrATTMMlSLM-Hr20)/ChI30-HT2  0) 

ETASI-(HIGtN3-HT20)/(Hl 30-HT20) 

PHlRLf ^0,^7639 
C 

0UM1‘5L  =hP20 

U  (CUPlbL.LL.0.0)  »ihIfC(6#*jOO>  OUPI5L 
If  tDUHl5L.LL.0.0)  Slop  lOl'J 
C 

PHir20i(hG*ALCG(hP2C)/J)4PhlHfF 

C 

OUMlbL=HPJO 

If  (CUM6L.Lb«0.a)  liPirC(6»^00)  DUPlbL 
If CUUHlbL-LL-0.0)  STOP  lOlfa 
L 

PHI  r50  =  (hG*  ALCCChP  30)  /  j)»PMlHf.f 
C***  CALCULATE  PH  1 1 S t ( HP  I  S£ T  ) 

C 

0tMl7L=hPlStT 

U  (0UHl/L-Lt.U,0)  b}ilU(b»‘>00>  0UM171 
If (UUH1/L,L£  .0.0)  SrCP  lOi7 

c 

c 

c 

PHII  6£ -CHG* ALLGChP I  St T>/J)*PHlb£f 
LF  APUL  =  (PHnSt-PHn20  )/(Phlli0*PHlT20  ) 

L 

C 

PS  rH;l=(PSTHHO*PSrHKL)/2.0 
AlhH-  lHHCAl«tJOirf  *21. 

G4=2./('SAMi*l.O) 

GS=.^«IGAH3*U0)/(GA^5-I.0) 

Gb=.‘)*(GAMi-l.O) 

i^OTHI-0.0 

hOFH2=G,C 

AH  FHH 

AHN4=1.0 

NUFM=MAV*SQRT(rT5QJ/(C01HR*ATHH*PSTHR) 
iS  CONFINUt 
L 
C 

OUMlrt^GAM  i«GC/KG 

U  (  GUf' 1  a  .L  F  .0.0)  kLIIL(b.SOO)  OUkl^ 
n cjUMid.LT.o.o)  siop  lOid 
c 

OLH  19  :  I .C ♦Gb •( AP  FPP**2  ) 

If (UUHI9.LT.0.0)  kRITL(b.SOO>  OUMI9 
If COUMI9.LT.O.O)  STOP 
C 

kJFMS-SQHF(GAMi«GC/RG)»AHFhR*SOHI(l.O^Cfc*(AHTHH**2)) 

1 f ( AbS( M Cl H  i-WUl P  ) -L r  .  ,CCO  I  )  CCTC  3f 

CALL  FTKAF2(kOFHl^AHM^IiCFH2»  ANN2  »M0  F  M  i  *  AH  T  HR  »  kO  F  N  »  AHk4) 

GIJFU  5S 
3/  CCbFIkLt 

PFFHK:PSIHH*((1.0*Gb*(AMlHN»*2))**(GAHj/(GAH3-1.0)>) 

LPUMPSIAU-PSlhPD/CfnPk-fSlhbU) 


i 


I'SAl'A-  I  K-S.'*-  ) 


/ 1  ouu 
IVUU 
JOO 
r^liJO 

'  <;  iuo 
r^k  uj 

'ilbOO 

fdtiOO 
'i?9uO 
/■  iouo 
?'ilUO 
'  i^OO 
/iiUJ 
/  i400 
r  i-aUQ 
t'ibuu 
nroQ 

f  itiUQ 
/  i9ua 
/4oaa 
?4lUO 
/<l^UO 
fk  iOO 
f  <i<*uO 
/4lOU 
/46U0 

/<«dUO 

/4900 

.»^ooo 

/‘H  00 
'b?UO 

/‘j4  00 

't>doa 

/■i/UO 

/•jooo 
'b9uO 
^bOOO 
/blOO 
^fa^UO 
^b  JUO 
/b«UO 
'  b  j  0  0 

/bbU  J 
^b  ^UO 
/  tduo 
i'bVOO 
'  /GdO 
f/\  00 
f  f 

n  luo 

'  f  4G0 

//'>dO 

//too 
nr^Q 
//ioo 
//^oo 
'0000 
/bl  00 
^0^00 
ioo 

/  04  00 

r«'>ao 
'  a  b  c  0 
/  d/  uo 
/ddOO 
^d900 


CI*L^(Pi>LXL-PjfHHL)/(Pnn«-PSrKKL) 
AAirHH^(b4»(Cl,U*Gb«(*HTHfi*  •-?))••&*>  ))/il»‘lH(i 
AAbf XU-AAGIHK«AHL 
AAGlXL^AASIHN^AHl 
AA^l =0.0 
AAb<e=G.Q 
AMtXU  =  -‘i 
AMN4=l-0 
4b  CCMINLt 
C 
C 

AAGi  =  Clj4«(Cl.4Ub*CAHt.XL»»^))»»bb))/AHtXU 
If  (AliG(AA5i-AA5LXU).LT..0UO|  )C0  TO  4/ 

CALL  ITKAT^(AASUAKM»AAS^#AHb<*#AASi>A^£XU#AAS£XU#AHN4) 

tiC  TL  4b 
4^  CUNTINUL 
C 
c 

AAbl =0.0 
AAb«?  =  0.0 
A^tXL  =.b 
AHN4=1.U 
bb  CONTINUE 

C 

c 

AAi>i  =  (u4*((1.0*06*(AHfXL**<^))«*li'5))/AH£XL 
If  I  ALbC  A  Ab  i- A  ASLXL  ).L  !  .  .COOl  >  OtlTC  bf 

CALL  lIKAT<'CAA5l#AHNWAA52#AHKi?#AA5^»A«£Xl»AAStXL*AHN4T 

oUTO  bb 

b'  CChllNUE 

C 

c 

HhLXU  =  (i.O*Gb*(AhEXU**i?))»*(-(C»Ar5/(OA>3-l.C)>) 

PbLXUl  :PnEXU*PTTHK 

PHlXL  =  (1.0*Lb*CAHfXL«»i^))*«(-(GAMS/(GA«3-l.0))> 
PbLXLl-Pbt  XL'PMMh 

If  »  U=CPblXC-PSrHhU)/CPbt  XL J-PbTHHU) 
tffL=(PSLXL-PSTHHL)/(HS£XLl-PbTHHL) 
PHLCb-<PbtXL-FiTHHU)/(PTIHh-PSlHhU) 

P.UCL=(P$f XL-PSTHfiLl/(PTTHH-PSTHKL) 
Prt£CLr=(PSrHH-Pb<?b)/(PI?b-PS^'b) 

PfiLE  =HbLE /P  T2b 
C 
C 
C 

0UM2d  =  ((PHLL*«(-(GAMi-U0)/CAH3)  1-1.0  >/Gb 
If (0UM20-Lf-0.O)  WRITCCb.bOOl  0Uk20 
If  (i;Uh20,l  T.O.O)  STOP  1020 

c 

AHLE=b4RKl(PHlf«*£-(tA^3-l.O)/OA»'3)  )-  I  -  J)/Gb) 

C 

0UM2l=r,  AP5/(Hu*UC) 

If  (LUP2  l.L  r-O.G)  ^SiltCb.bOO)  tUf'21 
If  (1)0/121. LT. 0-0)  blOP  1021 
C 

DCM22- 1 .04Gb«( AHlt  ‘*2  ) 

If  (OUf'22 -L  T.O  .0)  Kf- I  rtCb^bOO)  0UP22 
If (OUM22.LT.O.O)  STOP  1022 
c 

c 

c 

CH  Ui>aitHGAMi  )  )•  AP.1.1  •S9Km  .0«Gb«(*fl.l**2)) 

*fFF=M*lC5ilHI(ITlc)/tCHI«PSlt«GC) 
C0=AFFf/(<'.0»PI*HLF*HDIFF1 
ALPHLF  =  AKCOS(CC/COLF  ) 

AIACU=ALPHLF-AI.P*ll«H/l(lO. 
tAPKAI-AIIIh/AFFF 
PRTI iPI  iO/P  I  19 
r h I T  =  I r  JC/I I2C 

t[AAI)  =  (PHII**((GANC-l.l/(,APC»-l.l/(IHI|-U) 

L 

AL  FHPO-AIPHJO  •  leO  ./P I 
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/vuou 

8£I*<;U-iltT*^a*ldO./PI 

/VlUO 

b£  f A^s-a£ I A?s* lao./Pl 

rtzQO 

ALPH,?^-ALPH?i»10O./PI 

B£  -d£  I  A^b«  IdU./P  1 

^9400 

U£VI^0LVI«1«0./PI 

^9^00 

ALPHL£;ALPHL£»iaC./Pl 

^9t>U0 

AiNco-AiNco*iaa./Pi 

79700 

C 

79800 

OCCl J=AM19 

7990U 

DC(<?)=ALUS20 

30000 

OC(J  )=PIPQ 

30100 

DC(A)=ALPH^O 

30^00 

OCC^  )=AM<,'J 

30SOO 

DCCb  )-\<Z0 

80^00 

DC(£)=VUJO 

30*300 

DC{a)  =  rfMi'0 

aoboo 

UC  C9 )  =  UPOM 

80700 

octici  =  a£rA?(i 

80800 

octii 

80900 

act )  =  AM^ShU 

aiooQ 

UC113)=V20SMO 

81100 

OCtlA  )=3£TA?S 

^1^00 

DCtl5)=KH2SHa 

dliuo 

DCC  lh)iU£LH 

A  UOO 

QC(1/I=0£LHC 

81^00 

DC(18)=ULTUfC 

81600 

DC(19)=U2^ 

81700 

UC(£0  )=AM?5 

81800 

UCtZU^ALPH^b 

819U0 

DCt22)  =  fUIA2i 

8<;ooo 

UC(211=KP2‘j 

8^100 

UC(2A)=PHriI 

8^200 

0C(Z5  )  =  Ph5I I 

d^iOO 

DC(2b)=DLVI 

8^400 

UC(2£>=£IAt 

82500 

0C(2a)-M£Ht 

8^600 

0C(2»>-HHt3£ 

82700 

aC(50)=HAV0P 

82800 

OCC31)=Vl2bA 

82900 

0Ct32)=SP 

diOOO 

oc(  i3)=0LrbiA 

8  ilUO 

DC(3A)=aiT0T 

31200 

UC(  5")  l  =  PI2‘j 

diiOO 

&C(3b)-PH£C?9 

814  00 

0C(3£)-AL0529 

81500 

JClSa  )  =  AM2V 

81600 

OCC 39)=tIAb29 

81700 

OCCAO  )=PK£C  50 

81800 

DCIAIJ^ALOSSO 

81900 

OC(A2)=APiO 

84000 

UCCA  3  )  =  £TA030 

84100 

UCtAA)=PfiSr 

34200 

0ClAb)=£IA5I 

84  100 

DC(Ab)=£lAPaL 

34400 

0C(A/)=P5IHR 

34500 

UC(A8)-P  rlHR 

34600 

DCCA  9  >  =CPU 

84700 

BC(^0 J=CPL 

84600 

DC(5l )=£FFU 

84900 

0Ctb2)=EfFL 

85000 

UCta3)=PR£CU 

85100 

UCt^A )=Ph£CL 

85200 

OCCbbl=PR£CLT 

85100 

UC(‘>b  )  =  AMLE 

85400 

DC<5£)=ALPhL£ 

85500 

DC(^a)=AINCD 

85600 

QCC99)=CAPRAI 

85700 

UC(bO)-HPUf 

85800 

DCCbl  )^PNTl 

85900 

0C(b2 )=THTr 

L'SAi-'A-TR-84-5 


Appendix  B 

Sample  Equations 


■OmMl  9*  C4LCUI.«ltO 


Mtt 

•LOSPO 

ruo 

«LPM20 

AII20 

VU20 

9PP0 

U20M 

•CIAPO 

•«S}?C*00 

.tooai'Oi 

•  1’  Scl«02 

.tQCOC*Ot 

•4444C«0C 

.9C12C-01 

.SUSE40S 

•99SSC*D1 

^  .4PSIC*0P 

IPtSNC 

•  IU2S 

•tiasHo 

DElii 

OtLPC 

OLTDTC 

UPS 

AMPS 

•  I0iu«0l 

.S41PI400 

•)«4SC«0J 

.6«ooc«e? 

•1536C«OI 

• 1?P9E*05 

•l2t4C40S 

•954SC«00 

•t0IPC*04 

•ll4AC*OI 

011*21 

ill2S 

putfl 

PUSH 

0C«1 

CTAT 

II2MI 

NS|S£ 

PAvOf 

.Mtii<a2 

«S29iC*0Q 

•«|42C*OI 

•4cs7r«ci 

»S?0bE«0l 

•904PC«Oe 

•|P2«C*01 

.404SC40P 

msA 

or 

otror4 

oiior 

PIPS 

4L0SPV 

»»2* 

CHOP* 

PBCCIO 

.•MJ[*00 

«i4flie«oo 

•1SS5C«03 

.fc9|2t*00 

•P1ISC»00 

.ZtlOt.OO 

•r4pBr*oo 

•AvaPC*eo 

not 

MSO 

CIAOSO 

fUJI 

C14ST 

tTIPCL 

PStHR 

PIIMO 

CPU 

Cpl 

.22  0*1  *00 

•763iC*oa 

•  TH5C*Ot 

•6C«1C*00 

.b;s(e«oo 

•nuc»os 

.It2tt*0> 

•  SCUC400 

•tSCiAC^OO 

crru 

irri 

pnccu 

PfiECL 

PPtClf 

ixif 

4LPHLC 

*iaco 

CAPPAT 

Mpor 

•ro}*c«oo 

.«II2C*0« 

«?C  t<( «00 

•4J44C«00 

•46P4C400 

.  JU4IU0I 

•  2J5t£  *0^ 

.iiPArtoi 

»I04«C»0P 

pail  rail  tiaii  iiaao 

.•0141*01  .ifi«c*oi  .ajeo^aco  .aiiae«oo 


THE  EQUATION  USED  TO  DETERMINE  THE  UNCERTAINTIES  IS 
U  =  B  +  2  +  S 

B  IS  DEFINED  AS  THE  BIAS  LIMIT,  AND  S  IS  DEFINED  AS  THE  PRECISION 
INDEX  THE  VALUE  OF  T95(2)  COMES  FROM  THE  TWO-TAILED  STUDENTS 
T  DISTRIBUTION 
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ONE -DIMENSIONAL  UNSTEADY  ISENTROPIC  GAS  FLOW  (Part  II) 

Daniel  H.  Daley* 

Ed  i  tor ' s  Note 

This  article  is  the  second  ii  a  two-par*  tutorial  series  presenting 
new  insights  into  familiar  topics.  Part  I  .1  the  last  issue  of  the 
Digest)  examined  the  wave  equation  and  its  solutions.  Part  II  applies 
the  wave  equation  to  unsteady  isentropic  flow  of  a  perfect  gas. 

I  .  Nonlinear,  One - D  i  mens i onal ,  Unsteady  Flow 

The  following  equations  govern  nonlinear,  one-dimensional,  unsteady 

flow: 

Continuity, 


Euler's  equat ion , 


30 

3P 

3u 

u 

-*■  P 

=  0 

(  1  ) 

3t 

3x 

3x 

3u 

3u 

1 

3p 

*  U 

— 

=  0 

(2) 

3t 

3x 

P 

3x 

and  the  Isentropic  relations. 


•  7  •“ 
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In  th,  tMi  n  n  t  i n ,  :■  p  r  n  ;3  a  r- e  ,  u  =  velocity  in  the  x  direction,  o  = 

density,  T  =  t  ert  p  ^  r  1 1  u  r  e ,  c  =  the  speed  of  -sound,  and  k  =  ratio  of 
s  p  e  0  1  f  i  0  i'.  p  a  t  s  .  The  subscript  1  designates  some  known  reference  state. 
This  set  of  five  equations  relates  the  five  dependent  variables  (o,  u, 
p,  T,  and  c)  to  the  two  independent  variables  of  position  and  time 
(x,t).  Since  p  is  a  function  of  p  only  as  given  by  Eqn.  (3)t  we  can 
write  Dp/'X  in  Eqn.  (2)  as 


3p 

dp 

ap 

3x 

3x 

3x 

where  c“  =  dp/dp.  When  we  substitute  the  above  expression  for  3p/3x  in 
Eqn.  (?),  the  result  is  three  equations  in  the  three  unknowns  (p,  u, 
and  c)  —  Eqns.  (1),  (6)  [modified  from  Eqn.  (2)1,  and  (5)  —  as 
follows: 

3p  3p  3u  ^ 

3t  3x  3x 


2 

3u  3u  c  3p 

3t 


c  = 

These  equations  are  exact  as  long  as  the  motion  under  consideration  is 
isentroplc.  Under  that  condition,  a  solution  of  the  equations  would 
completely  describe  p,  u,  and  c  In  terms  of  x  and  t  for  a  system  of  gas 
contained  in  a  constant  area  tube  (Figure  1). 


We  can  then  find  the  pressure  and  temperature  through  Eqns.  (3)  and  ( 4  )  . 
Since  terms  such  as  u  (Ou/lx)  make  the  equations  nonlinear,  however, 
they  are  not  easily  solved.  There  are  two  methods  for  arriving  at  a 
solution.  One  is  based  on  sm'^11  perturbations  in  p  and  u  and  leads  to 
the  classical  wave  equation.  The  other,  the  method  of  c h a r ac t er  i  s t  i  cs  , 
originated  by  Riemann  and  Earnshaw  in  the  nineteenth  century,  leads  to 
a  stepwise  solution  of  the  nonlinear  equations,  a  procedure  that  is  more 
exact  than  the  linear  solution.  In  the  next  section,  we  will  linearize 
the  equations  to  obtain  the  so-called  acoustic  equations. 


1 1 .  Linearized  One - D i mens  i  ona 1  Unsteady  Flow 

Consider  small  perturbations  in  density  and  velocity  from  an 
initial  uniform  condition  with  property  values  of  o^,  p^,  T^,  and  c^  and 
zero  stream  velocity.  The  small  perturbation  in  density  is  :  =  u  -  , 

and  the  small  perturbation  in  speed  is  '5u  =  u  -  0  =  u.  Here  we 
introduce  a  non-dimensional  perturbation  density,  .7  ,  defined  by 

P-Pj^  5p 

Pi  Pi 


so 

I 


that 


P  =  Pi  (  1  +  P  ) 
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Substituting  as  given  by  this  expression  into  Eqns.  M),  (b),  ani  (b) 

and  placing  Eqn.  (5)  into  Eqn.  (6),  the  equations  hecorr,'?' 


3p 

3u 

3P 

+  1  i 

3u 

+  p 

3t 

3x 

3x 

3x 

Su  3u  „  _  u  o 

-  +  u  -  +  C,  (  1  +  P  -  =  0 

3t  3x  3x 


We  can  now  linearize.  This  means  neglecting  terms  that  are  small 
compared  to  ot!,er  terras  in  the  equation.  Thus,  we  drop  terms  that 
contain  products  of  small  perturbation  quantities.  In  addition,  since 
p  is  small  compared  to  1 ,  we  drop  Z  relative  to  unity  in  the  last  term 
of  the  second  equation.  These  operations  give  the  following  linear 
equat ions : 


3p 

3u 

3t 

3x 

3u 

2 

+  r  ^ 

3t 

^  3X 

2  2 

Since  3  u/axlt  =  3  u/3t3x,  we  can  differentiate  the  first  equation  with 

respect  to  t  and  the  second  with  respect  to  x  and  equate  them  to  get 


c 


2 

1 


a'p- 


( 


1  65 


■  ’  - 


Ia  ■JV.O  -v*  O.  ' 
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Similarly,  because 


r/3x'it  =  '  ,, /'if'X,  we  obtain 


3^u 


( 1 


-1 


3x^  3t^ 

which  was  discussed  in  Part  I  of  this  series.  Eqns.  (9)  and  (10), 
which  are  in  the  form  of  the  wave  equation,  show  that  any  disturbances 
in  density  and  speed  propagate  through  the  gas  without  change  at  the 
speed  of  sound  in  the  undisturbed  gas,  c^  . 

If  we  write  the  general  solution  to  Eqn.  (10)  as 


u  =  6u  =  f(x-c^t)  +  g(x  +  Cjt) 


( 1 


then,  since  p  is  related  to  u  through  Eqns.  (7)  and  (8),  the  solution  to 
Eqn .  ( 9 )  must  be 


6p 


1 


P  = 


f  (  X  -  t  )  =  —  g  (  X  t  t  ) 


( I 


One  can  verify  that  Eqns.  (7)  through  (10)  are  all  satisfied  by  Eqns. 
(11)  and  (12).  If  we  consider  only  simple  rightward-  and  leftward- 
propagating  waves  (for  which  g  =  0  and  f  =  0,  respectively)  we  can  write 
Eqns.  (11)  and  (12) 

6p  =  6u 


P 


1 


=  + 


(  1 


where  the  +  and  -  signs  refer  to  rightward-  and  1 ef t war d -pr opaga t  i  ng 
simple  waves. 

Here  is  the  interpretation  of  Eqn.  (13).  If  the  piston  (Figure  1) 
imparts  a  positive  incremental  change  in  velocity,  u,  to  the  gas 
particles  adjacent  to  the  piston  face,  this  disturbance  will  prop  .a  gate 
to  the  right  through  the  gas  at  the  speed  Cj.  This  disturbance  in 
particle  velocity  is  accompanied  by  a  positive  increment  in  gas  density 
given  by  Eqn.  (13).  If  we  move  the  piston  to  the  left,  imparting  a 


negative  increment  in  iu,  then  this  rightward-propagating  disturbance 
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in  particle  speed  is  accompanied  by  a  decrease  in  gas  density.  .Similar 
remarks  apply  to  the  density  effect  produced  by  leftward-  [negative  sign 
in  Eqn.  (13)1  propagating  disturbances  in  particle  speed  generated  by 
a  piston  placed  at  the  right-hand  end  of  the  tube  (Figure  1). 

For  example,  suppose  the  piston  of  Figure  1  imparts  a  sinusoidal 
velocity  disturbance  to  the  air  particles  adjacent  to  the  piston  face, 
so  that  at  t  =  0  the  gas  velocity  distribution  in  the  tube  is 

ft  X 

6u  =  10  -  sin  2  it  — — 

sec  4  ft 

0  £  X  £  2  ft. 

The  perturbation  density  distribution  in  the  gas  is  then  given  by  Eqn. 
(13)  with  the  plus  sign.  Thus,  if  c^  =  1000  ft/sec  and  =  0.076 

•3 

Ibm/ft  ,  corresponding  to  standard  air  conditions,  then 

6u 

6p  =  - 


Ibm 

=  0.076  — 
ft-’ 


ir,ft  2  It  X 

10 - sin  -I  -Ti.*— 

sec  4  ft. 


1000  ft/sec 


Ibm  2  It  X 

=  0.00076  ■  sin  ~  (3  t  =  0 

ft-’  4  ft. 


These  disturbances  in  7u  and  5p  propagate  through  the  gas  along  the 


rightward  characteristics  (Figure  2). 


:■  we  must  also  find  the  pe  r  t  u  r  ba  1 1  o  n  s  in  p  r  es  s  u  r  ••  and  t  em  po  r  a  t  u  r  o 
a  0 'a :  m  p  a  n  y  i  ri  p  the  density  and  pa  r  t  i  o  1  e  -  s  pe  e  d  d  i  s.  t  u  r  t  a  n  c  g  s  .  These  ar'" 
piven  t.  y  the  i  sen  tropic  relations,  Ect.  s.  (3)  and  ( '!  )  ,  which  tak"  the 
1  1  0  w  i  n  g  f  o  r  m  when  written  in  terms  of  trie  perturbation  density,  , 

/  Pi  \  k  Pi  r  r  -  1 

p  =  I  — ^  IP  =  — ^  Pl  (1  p)  “Pi  ••• 

\^Pi  y  Pi  L  j  L  j 

-r,  "'^cppi  r.g  terms  of  higher  order  than  q  in  the  above  binomial 


e  X  a  r  3  i  o  r.  , 


p  =  p,  (  1  +  k  p  ) 


w  ‘  can  write  the  equation,  using  'S  p  =  p  -  p  ^  ,  a  £ 

6p  5p 


i  rr,  i  :  ■  I  i  y  ,  E  q  n  .  (  i)  )  gives 


=  (  k  -  1 


Jsinv  these  results,  W'-^  can  determine  the  pressure  and  temperature 
listurhances  associated  with  the  velocity  and  density  disturbances 
illustrated  by  Figure  2.  According  to  Eqn.  (Ill),  when  p^  =  1680 
It  f/ft"  at  time  t  =  0, 


6p  =  kp 


6p  =  1,4(1680  psfa)  10”^  sin  2ti  — — 

L  J  4  FT. 


6p  =  23.5  sin  2Tt  — — 

ft^  4  ft. 

id  a  c  r;  r  d  i  n  g  to  Eqn.  f  1  S  ;  ,  if  T  |  =  1)1  5  degrees  R  ,  then  at  t 

"■  m  '1  r  a  t  u  r  d  i  s  t  u  r  b  a  n  c  e  i  s 


=  0  the 
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6T  =  (0.4)  (415°R)  10  sin  2ir 


6T  =  1.66“R  sin  2ir 


We  can  locate  these  disturbances  graphically  at  any  time  t  by  reference 
to  the  characteristics  of  Figure  2.  But  in  this  particular  example,  it 
is  Just  as  easy  to  write  the  equation  for  each  property  disturbance  in 
the  gas.  Since  we  are  dealing  only  with  rightward-traveling  waves,  eact 
disturbance  must  be  written  using  the  argument  (x  -  Cjt).  In  terms  of 
this  variable  the  property  disturbances  are 


6u  =  10  ft/sec  sin  (x  -  c,  t) 

4  ft  ^ 


6p  =  0.00076  Ibm/ft^  sin  -il-  (x  -  c,  t) 

4  ft  ^ 


6p  =  2.35  Ibf/ft^  sin  -22-  (x  -  c,  t) 

4  ft 


6T  =  1.66‘’R  sin  (x  -  c,  t) 

4  ft 


This  example  describes  rightward-propagating  simple  waves.  But  we 
can  obtain  similar  results  for  leftward -propagating  simple  waves  by 
letting  a  minus  sign  represent  the  relation  between  density  and  particle 
velocity  in  Eqn.  (13). 


I V .  Finite  Isentropic  Waves 

The  essence  of  the  linearized  theory  presented  in  the  preceding 
section  lies  in  the  concept  (resulting  from  the  linearization  of  the 
governing  equations)  that  small  per tur iia t i ons  in  gas  properties  travel 


through  the  gas  at  a  constant  wave  speed. 


■  I  • 


In  relation  to  wave 


speed,  we  igir^re  the  effects  of  variation  in  the  speed  of  sound  and 


gas-particle  velocity.  But  if  the  gas  disturbances  are  not  sufficiently 
small,  we  must  allow  for  the  change  in  speed  of  the  wave  from  point  to 
point. 

Consider  an  isentropic  but  finite  r i ght war d - 1 rave  1 i ng  compression 
wave  in  a  gas  (Figure  3).  We  can  approximate  the  continuous  curves 
forming  the  waves  of  u,  u  ,  and  c  versus  position  in  the  manner  of 
calculus  by  the  infinitesimal  step  "curves"  shown.  In  this  procedure, 
the  disturbances  propagate  through  the  gas  as  an  infinite  number  of 
infinitesimal  pulses.  Because  each  pulse  produces  infinitesimal  changes 
(du,  dp,  dc,  etc.)  in  the  gas  properties,  it  travels  into  the  gas  in 
front  of  the  pulse  at  the  speed  of  sound  c  (always  plus)  in  that  gas. 

The  gas  itself  has  a  velocity  u  (vector  quantity).  Relative  to  the 
stationary  wall  of  the  tube  containing  the  gas,  the  pulse  speed  or  local 
wave  speed  w  (always  positive)  is 

w  =  C  +  u  (16 

where  the  plus  sign  applies  to  r i ghtwar d -tr a ve 1 i ng  waves  (Figure  3)  and 
the  minus  sign  applies  to  leftward-moving  waves.  Since  c  and  u  vary 
throughout  the  wave,  the  wave  speed  also  varies  from  point  to  point  on 
the  wave. 

We  can  find  the  wave  speed  at  any  point  on  an  isentropic  finite 
wave  in  terms  of  the  gas  velocity.  First,  the  local  speed  of  sound  is 


k-1 


Second  , 


du  -  ^  C  (  1  8 

where  the  plus  sign  refers  to  rightward-  and  the  minus  sign  to 
1 e f t war d -mo V  i  ng  waves.  We  obtain  this  result  by  application  of  the  mass 


conservation  principle  to  flow  through  a  control  surface  surrounding  and 
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Eqn.  18  is  similar  to  Eqn.  (13),  which  we  developed  by  applying  linear 
theory  (see  Section  II).  The  two  equations  are  identical  if  Eqn.  (13) 
is  applied  locally  so  that  pj  =  p  and  c^  =  c. 

Since  c  is  known  in  terms  of  o  (from  the  isentropic  relation  given 
by  Eqn.  (17)1,  Eqn.  (18)  can  be  integrated  to  obtain  u. 


1-^  i 


k-1 

1  T 

(p^t 


k-i 

T  1 


k-1  (  k-1 

T  'T 


u  ..  ^  (C-Cj) 

If  we  solve  this  equation  for  c  and  then  substitute  it  into  Eqn.  (i6), 
the  wave  speed  at  any  point  on  the  wave  (in  terms  of  the  speed  of  sound 


in  the  undisturbed  gas,  Cj,  and  local  gas  velocity,  u)  becomes 

w  =  ^  “7”  ^  ^ * 

In  the  compression  portion  of  a  wave,  u  increases  according  to  Eqn 

(18),  since  the  pressure  and  density  increase  as  shown  by  Eqn.  (3). 
Because  of  this,  wave  speed  Increases  as  the  wave  moves  through  the  gas 
and  thus  trailing  parts  of  the  wave  will  travel  faster  than  leading 
parts.  This  causes  the  wave  to  steepen  as  the  trailing  parts  overtake 
the  leading  parts.  We  hear  this  steepening  as  a  sonic  boom  when  an 
aircraft  traveling  faster  than  the  speed  of  sound  flies  overhead. 
Conversely,  the  opposite  effects  arise  in  the  rarefaction  portion  of  a 


This  article  demonstrates  how  the  equations  governing  non-linear, 
one -d i mens  1  on  a  1 ,  unsteady,  Isentropic  flow  can  be  reduced,  or 
linearized,  when  variations  in  properties  of  the  gas  are  small.  The 
governing  equations  have  the  same  form  as  the  classical  wave  equation. 
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Thus,  disturbances  in  velocity,  density,  pressure,  and  temperature  will 
propagate  through  a  gas  as  waves  with  the  properties  described  in  Part 


I  of  this  series. 
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THE  ENGINEER’S  BOOKSHELF: 

THE  MYTH  OF  OBJECTIVITY  IN  THE  WRITING  OF  ENGINEERS 
Robert  M.  Hogge* 

Picture  in  your  mind  a  wind  tunnel  in  an  aeronautics  laboratory.  A 
scientist,  dressed  in  a  white  lab  jacket  and  wearing  horn-rimmed 
glasses,  enters  the  observation  area.  With  a  professional  air,  he  scans 
the  equipment,  makes  sure  everything  is  ready  for  the  test.  Jots  down  a 
few  notes  on  a  clipboard  checklist,  then  flips  a  switch  to  begin  the 
experiment.  During  the  test,  he  looks  at  gauges  and  dials,  observing 
exactly  how  his  model  (a  forward-swept-wing  aircraft)  behaves  in  0.8 
Mach  coinditions.  Then  when  the  test  is  over,  he  reviews  his  findings 
and  writes  up  his  report.  He'd  like  to  write  his  report  completely  in 
the  language  of  math,  but  he  knows  he'll  also  have  to  use  the  English 
language  to  communicate  his  findings  to  some  members  of  his  target 
audience.  So  he  picks  his  words  carefully,  striving  for  absolute 
objectivity  .  avoiding  personal  pronouns,  using  passive  verbs,  and 
writing  in  what  he  appropriately  feels  is  the  "scientific  style."  He 
then  reads  the  paper  to  himself.  It's  based  on  empirical  evidence;  his 
reasoning  is  sound;  and  other  scientists  will  be  able  to  replicate  his 
experiment.  So  his  writing  must  be  objective.  Or  is  it? 

Of  course  it  is,  you  say.  Scientists  don't  write  about  their 
feelings,  opinions,  and  emotions.  These  "seekers  after  truth"  focus  on 
objective  reality  (the  empirical  evidence)  and  record  simply  and 
uniquely  what  they  detect  through  the  five  senses.  That's  perfect 
objectivity  .  isn't  it?  But  how  well  does  the  English  language  serve  as 
a  vehicle  for  that  objectivity?  Two  scientists,  Cox  and  Roland,  have 
some  grave  concerns  about  the  written  language  and  what  kinds  of  words 
scientists  use  to  present  their  findings:  "Since  scientists  agree  that 
their  observations  and  conclusions  should  be  presented  as  objectively  as 
’Major,  USAF,  Associate  Professor  of  English,  USAF  Academy 
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possible,  rhetoric  should  be  avoided  assiduously  in  scientific  writing." 
(Ref.  1  ,  page  14  0)  Here  Cox  and  Roland  imply  that  subjectivity  evir.ts 
in  scientific  writing,  and  one  of  the  great  culprits  is  "rhetoric," 
language  designed  to  persuade  or  impress  an  audience.  So,  like 
Puritans,  scientists  attempt  to  strip  rhetoric  from  their  specialized 
use  of  language,  usually  ending  up  with  a  cold,  lifeless,  dreary,  and 
boring  set  of  verbal  alternatives. 

But  where  did  we  get  this  scientific  penchant  for  objectivity  in 
writing?  Probably  it  originated  in  the  Greek  society  with  Plato  who  had 
no  room  in  his  Republic  for  men  of  literature  (those  who  used  words 
imaginatively  and  emotionally).  Or  maybe  it  came  as  a  reaction  against 
the  Sophists  who,  by  clever  and  intricate  arguing,  could  make  any  set  of 
facts  seem  true.  Or  perhaps  it  was  from  Aristotle,  the  patron  saint  of 
scientists,  who  argued  that  the  mind  is  simply  a  passive  receptor  and 
recorder  of  external  stimuli.  Or  maybe  it  originated  later  in  France 
with  Descartes  (the  Cartesian  method  of  analysis)  or  Ramus  (who 
separated  rhetoric  from  dialectic  and  argued  for  an  objective  language 
free  from  emotion).  Whatever  its  origins,  "objective"  writing  became 
thoroughly  entrenched  in  the  minds  of  scientists,  almost  a  sacred  ritual 
they  passed  on  from  generation  to  generation.  Today  many  scientists  and 
engineers  still  firmly  hold  to  this  "objective"  theory  of  writing. 

But  there  is  another  (and  more  expansive)  way  to  look  at  writing. 

To  show  you  what  I  mean,  let’s  look  at  a  major  ph i 1 osoph i ca 1  debate 
between  David  Hume  (1711-1766)  and  Immanuel  Kant  (1724-1804).  Hume,  an 
empiricist,  argued  that  all  thought  comes  to  us  through  the  five 
sen3es--that  the  mind  is  a  passive  observer  and  receptor.  This  is  the 
philosophical  position  many  engineers  adhere  to  when  they  write.  Kant 
too  believes  that  all  knowledge  begins  with  sensory  experience  but  that 
a  separate  phenomenon  (what  he  calls  a  priori  knowledge)  exists  already 
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in  the  mature  mind  and  does  not  come  from  the  senses.  As  examples  of  ^ 
priori  knowledge,  Kant  specifies  space  and  time,  two  aspects  of  reality 
the  mind  understands  intuitively.  Because  of  our  awareness  of  space  and 
time,  the  mind  is  able  to  interpret  and  arrange  sensory  data  into 
meaningful  patterns. 

If  Kant  is  correct  in  his  theorizing,  what  does  it  mean  for 
scientists  and  engineers  today?  Robert  Pirsig  provides  an  interesting 
explanation; 

If  you  presume  that  a  priori  concepts  in  our  heads  are  independent 
of  what  we  see  and  actually  screen  what  we  see,  this  means  that  you  are 
taking  the  old  Aristotelian  concept  of  scientific  man  as  a  passive 
observer,  a  'blank  tablet,'  and  truly  turning  this  concept  inside  out 
(Ref.  2,  page  119). 

Does  this  mean  that  the  whole  tradition  of  "objective"  writing 
might  be  in  error?  Throughout  the  years,  have  scientists  and  engineers 
limited  themselves  unnecessarily  by  clinging  to  a  seemingly  passive, 
impartial,  disinterested,  or  purely  objective  mode  of  writing?  Can  the 
human  mind  do  more  than  passively  record  external  stimuli?  Yes!  The 
mind  not  only  receives,  but  it  also  creates. 

To  show  how  the  mind  creates,  Robert  Pirsig  describes  an  analytical 
process  you're  already  familiar  with— The  Scientific  Method.  First  he 
lists  the  six  steps: 

1.  statement  of  the  problem, 

2.  hypotheses  as  to  the  cause  of  the  problem, 

3.  experiments  designed  to  test  each  hypothesis, 

A.  predicted  results  of  the  experiments, 

5.  observed  results  of  the  experiments, 

6.  conclusions  from  the  results  of  the  experiments. 

(Ref.  2,  page  93) 

Using  the  Scientific  Method  is  a  logical  and  empirical  way  to  arrive  at 
truth.  Or  is  it?  How  does  the  scientist  formulate  hypotheses?  Where 
do  they  originate?  Albert  Einstein  developed  a  law  to  answer  these 
questions:  "The  number  of  rational  hypotheses  that  can  explain  any 

given  phenomenon  is  infinite."  (Ref.  2,  page  100)  If  Einstein's  theory 
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is  true,  then  do  we  ever  arrive  at  one  truth  or  one  r(',i!ity?  hoes  the 
human  mind  have  the  awesome  responsibility  to  respond  to  these  infinite 
hypotheses,  selecting  those  that  are  the  most  credible  and  prioritizing 
them?  Yes,  the  mind  is  so  much  more  than  a  passive  receptor.  Thomas  S 
Kuhn  describes  the  almost  mystical  and  often  erratic  workings  of  the 
mind  as  it  develops  hypotheses  to  solve  scientific  puzzles: 

Scientists  then  often  speak  of  the  'scales  falling  from  the  eyes' 
or  of  the  'lightning  flash'  that  'inundates'  a  previously  obscure 
puzzle,  enabling  its  components  to  be  seen  in  a  new  way  that  for  the 
first  time  permits  its  solution  (Ref.  3,  page  122). 

So  science  does  not  always  move  forward  "objectively,"  by  induction  or 
deduction.  Occasionally  (and  in  the  best  of  situations)  the  mind 
provides  insights  intuitively. 

The  mind,  in  other  words,  creates  or  shapes  reality.  To  explain 
how  the  mind  works,  Kuhn  discusses  his  concept  of  "paradigms."  For  him 
a  paradigm  is  a  perception  or  a  model  for  what's  out  there — reality. 

For  example,  he  shows  how  astronomy  first  of  all  matured  within  a 
Ptolemaic  paradigm,  the  philosophy  that  the  universe  is  geocentric 
(earth  centered).  Ptolemy  not  only  perceived  it  that  way,  but  he  also 
developed  an  extensive  mathematical  construct  to  support  that 
perception.  Throughout  the  years,  different  scientists  suggested 
alternative  world  views.  But  it  wasn't  until  scientists  (with  advanced 
mathematics  and  instrumentation)  revealed  a  major  flaw  in  Ptolemy's 
paradigm  that  the  scientific  world  looked  for  another  way  to  interpret 
reality.  Thus,  the  Copernican  Revolution.  Then,  Kepler,  Newton,  and 
Einstein — each  with  a  different  perception  of  what's  out  there. 

How  the  mind  perceives  reality  makes  all  the  difference.  And  we 
don't  always  perceive  stimuli  In  the  same  way.  According  to  Kuhn, 
researchers  in  modern  psychology  have  shown  that  "two  men  with  the  same 
retinal  impressions  can  see  different  things."  (Ref.  3,  126-127)  Also, 

in  1938  Albert  Einstein  said:  "Physical  concept'’  are  free  creations  of 
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the  human  mind,  and  are  not,  however  it  may  seem,  uniquely  determined  by 
the  external  world."  (Ref.  page  8)  In  other  words,  Isaac  Newton's 

mind  created  the  Law  of  Gravity  as  it  relates  to  interpla.etary  motion; 
that  specific  law  or  concept  didn't  exist  before  Newton.  And  it's  the 
same  with  all  other  laws,  hierarchies,  and  concepts  in  science  —  all  are 
inventions  of  the  human  mind. 

How  "objective"  is  the  human  mind?  When  we  move  especially  from 
the  ma cr omo 1 ec u 1 ar  (Newtonian)  universe  to  the  m i cr omo 1 e c u 1  a r  (quantum 
mechanics),  objectivity  does  not  exist.  John  Wheeler,  a  physicist  at 
Princeton,  explains 

'Participator'  is  the  incontrovertible  new  concept  given  by  quantum 
mechanics.  It  strikes  down  the  term  'observer'  of  classical  theory,  the 
man  who  stands  safely  behind  the  thick  glass  wall  and  watches  what  goes 
on  without  taking  part  (Ref.  4,  page  29). 

But  I  submit  to  you  that  the  "participator"  mode  is  not  Just 
reserved  for  those  studying  quantum  mechanics.  It  also  applies  to 
engineers  working  in  aeronautics  labs,  with  computers,  on  construction 
projects.  As  engineers,  we  are  much  more  than  passive  observers  and 
recorders  of  objective  data.  Our  minds  shape,  mold,  organize,  and  even 
create  what's  out  there. 

It's  time  now  for  a  Copernican  Revolution  in  scientific  writing. 
Let's  leave  behind  once  and  for  all  the  Ptolemaic  notion  of  "objective" 
writing.  Instead  let's  use  the  English  language  fully,  uniting  rhetoric 
and  dialectic,  to  create  even  more  effective  "paradigms"  in 
engineering  —  new  discoveries,  new  ways  to  see  the  world. 
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